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ABSTRACT 


This  thesis  investigates  the  estimation  of  Lanchester 
attrition- rate  coefficients  in  a  deterainistic  aggregated 
combat  sod el.  Numerical  values  for  the  attrition-rate  coef¬ 
ficients  in  the  Lanchester-t ype  model  are  taken  to  be  the 
same  as  those  for  a  corresponding  continuous-time  Harkov- 
chain  model,  and  maximum  likelihood  estimators  are  developed 
for  these  Harkov-chain  coefficients.  A  discussion  of  the 
historical  background  of  Lanchester’s  equations  preceeds  the 
theoretical  development  of  components  of  the  attrition  rate 
equations.  Using  the  functional  forms  and  procedures  devel¬ 
oped  by  Gordon  N.  Clark  in  his  doctoral  thesis,  a  simple 
computer  simulation  program  is  developed  for  a  short  dura¬ 
tion  battle  consisting  of  homogeneous  forces.  Recommenda¬ 
tions  are  made  to  araify  this  program  to  model  heterogeneous 
forces  in  a  similar  battle.  A  discussion  of  the  theoretical 
background  supports  these  recommendations. 
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I.  INTRO DOCT ION 


The  purpose  of  this  thesis  is  to  investigate  methods  for 
estimating  Lanchester  attrition-rate  coefficients,  used  in  a 
deterministic  aggregated  combat  model,  from  data  generated 
by  the  detailed,  high  resolution  Monte  Carlo  combat  simula¬ 
tion  model  called  STAR  (Simulation  of  Tactical  Alternative 
Responses) .  An  additional  purpose  is  to  construct  a  simple 
Lanchester-type  combat  model  capable  of  using  the  estimated 
attrition  rate  coefficients  to  predict  the  outcome  of  STAR 
simulations. 

STAR  is  a  model  which  has  been  developed  primarily  by 
O.S.  Army  students  attending  the  Naval  Postgraduate  School. 
It  has  recently  gained  acceptance  in  the  Army  Community  as  a 
viable  combat  model,  useful  for  studies  such  as  the  United 
States  Army  Mortar  Study  to  be  conducted  in  1982.  The  cur¬ 
rent  version  of  STAR  is  written  in  SIMSCRIPT  and  requires 
about  2  megabytes  of  core  storage.  The  run  true  for  one 
replication  of  STAR  involving  a  company  size  force  defending 
against  a  battalion  size  attaching  force  is  approximately 
four  CPU  minutes  on  the  IBM  3033  computer. 

The  basic  conceptual  difficulties  found  in  this  thesis 
are  summarized  in  Figure  1.  The  numerical  values  for  the 


attrition  rate  coefficients  in  the  deterministic  Lanchester- 

t 

type  aodel  are  assumed  to  be  the  sane  as  those  for  a  corre¬ 
sponding  continuous-time  Harkov  chain  model.  With  this 
assumption,  maximum  likelihood  estimators  (HLE's)  can  be 
developed  for  the  Markov  chain  model  coefficients.  Specific 
data  needed  for  the  HLE  computations  is  extracted  from  the 
high  resolution  Monte  carlo  simulation  and  used  to  compute 
"estimated"  coefficient  values.  These  values  for  the 
continuous  time  Markov  chain  aodel  are  then  substituted  into 
the  "analogous"  deter ainistic  Lanchester-type  aodel. 

Constructing  a  Lanchester-type  attrition  model  sue  s 
that  proposed  in  this  thesis  is  of  significant  value.  <  a 
respect  to  the  hierarchy  of  models  concept  [fief.  1],  it 
could  possibly  fill  the  requirement  as  a  'hook*  or  link 
between  STAR  and  low  resolution  models  that  exist  at  higher 
levels  in  the  pyramid  cf  models.  However,  possibly  the  most 
important  contribution  this  aodel  can  make  is  as  a  filter 
aodel  for  STAR  to  reduce  costs  and  time  to  run  the  many  var¬ 
iations  of  a  scenario  in  STAR.  To  explain  further,  consider 
the  upcoming  mortar  study.  To  conduct  such  a  study  many 
runs/replications  of  a  battle  using  the  STAR  model  will  have 
to  be  made.  These  numerous  runs  are  required  to  compare 
various  types  of  mortars,  various  types  of  force  mixes. 
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Figure  1:  Relationships  of  Models 
weapon  sixes,  and  scenarios.  Since  STAR  is  a  Sonts  Carlo 
siaulation,  Tariance-reducing  replications  are  required. 


This  requires  significant  computer  time  at  a  significant 
cost.  By  using  a  siapler  analytical  model 'which  produces 
comparable  output  to  a  STAB  run,  many  force  mixes/weapon 
aixes/scenarios  can  be  eliainated  from  further  consideration 
at  a  minimal  cost.  STAB  can  then  be  used  to  concentrate 
only  on  the  more  promising  options  for  force/weapons  mixes. 

This  thesis  records  the  progress  made  in  developing  the 
small,  simple,  Lanchester-type  combat  model  discussed  above. 
It  also  lays  the  foundation  for  future  work  on  this  model. 
Chapter  2  of  this  thesis  provides  an  introduction  to  Lan- 
chester-type  combat  models.  Chapter  3  outlines  the  develop¬ 
ment  of  the  Deterministic  Aggregated  aodel  of  STAB  on  the 
Apple  computer  (DAHSTAC)  ,  which  represents  the  first  step  in 
building  an  analytical  model  comparable  to  STAB.  Chapter  4 
explores  various  Lanchester  attrition  rate  functional  forms 
which  may  be  used  in  the  DAHSTAC  or  in  future  generation 
models,  chapter  5  details  some  desired  enhancements  of  the 
model,  while  Chapter  6  discusses  the  work  required  to  tran¬ 
sition  from  the  current  homogeneous  force  DAHSTAC  model  to  a 
more  realistic,  albeit  more  complicated  heterogeneous  model. 
Chapter  7  contains  the  concluding  remarks. 

Appendix  A  lists  the  variables  used  in  the  computer  pro¬ 
gram  for  the  DAHSTAC  model.  Appendix  B  lists  some  notes  of 


interest  concerning  the  prograa  and  model,  which  aay  not  be 
apparent  to  the  reader.  Appendix  C  discusses  the  utility 
prograas  developed  to  support  the  aodel.  Appendix  0  is  a 
discussion  of  additional  work  done  by  the  authors  in  the 
fora  of  a  coaputer  aodel  that  explicitly  plays  the  terrain 
effect.  Appendix  E  consists  of  work  done  by  two  other  stu¬ 
dents  on  a  postprocessor  which  supports  the  coaputational 
requirements  of  the  aodel  and  acts  as  an  interface  between 
STAR  and  the  model.  Finally,  Appendix  F  contains  the  compu¬ 
ter  listings  of  all  prograas  written  for  this  thesis. 
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II.  HISTORICAL  BACKGROUND 

Lanchester-type  combat  models  are  analytical  models  in 
which  differential  equations  are  used  to  express  the  rate  of 
change  of  force  levels  to  represent  the  effects  of  the 
attrition  process.  The  concept  of  representing  attrition 
effects  in  land  warfare  with  such  mathematical  equations  was 
originated  by  F.  ».  Lanchester  in  1914  in  order  to  quantita¬ 
tively  justify  concentration  of  forces. 

Lanchester  hypothesized  that  under  conditions  of  "modern 
warfare"  attrition  between  two  homogeneous  forces,  each 
using  "aimed"  fire,  could  be  modelled  by: 

dx/dt  =  -ay  (t)  with  x(0)  =  x0 

dy/dt  =  -bx(t)  with  y(0)=  y# 

where  x  (t)  and  y  (t)  are  the  number  of  X  and  I  firers  alive 
at  time  t,  and  a  and  b  are  positive  constants  that  are 
called  Lanchester  attrition-rate  coefficients.  The  constant 
'a'  represents  the  effectiveness  of  an  x  firer  in  killing  y 
targets  per  unit  of  time  and  1 b'  represents  the  equivalent 
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effectiveness  of  the  y  fixers.  The  equations  written  above 
are  usually  called  "Lanchester • s  equations  for  aodern  war¬ 
fare".  Clearly,  these  functional  forms  for  attrition  rates 
that  Lanchester  developed  are  directly  proportional  to  the 
attrition  rate  coefficients  and  the  number  of  firers  fron 
the  opposing  force.  The  assumptions  that  aust  be  made  in 
order  to  use  Lanchester* s  aiaed  fire  equations  to  aodel  can- 
bat  attrition  are  very  restrictive.  These  assumptions  can 
be  stated  as  follows  £8ef.  2]: 

1.  Two  homogeneous  forces  are  engaged  in  coabat; 

2.  Each  unit  on  either  side  is  within  weapon  range  of  all 
units  on  the  other  side; 

3.  The  effects  of  successive  rounds  in  the  target  area 
are  independent; 

4.  Each  unit  is  sufficiently  well  aware  of  the  location 
and  condition  of  all  enemy  units  so  that  it  engages  only 
live  eneay  units  (one  at  a  time)  and  kills  them  at  a 
constant  rate,  which  does  not  depend  on  the  enemy  force 
level,  when  an  eneay  target  is  killed,  search  begins  for 
a  new  target,  with  the  rate  of  acquiring  a  new  eneay  tar¬ 
get  being  independent  of  the  eneay  force  level; 

5.  Pire  is  uniformly  distributed  over  surviving  enemy 
units. 
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Lanchester  also  hypothesized  that  shea  both  sides  used 
'•area"  fire,  the  attrition  of  forces  could  be  modelled  by: 

dx/dt  =  -axy 
dy/dt  =  -bxy 

The  first  three  conditions  under  which  Lanchester "s 
equations  for  area  fire  apply  are  identical  to  those  for  the 
equations  for  aodern  warfare,  but  assumptions  4  and  5  are 
replaced  by  [Ref.  3]: 

4.  Each  firing  unit  is  aware  only  of  the  general  area  in 
which  eneay  forces  are  located  and  fires  into  this  area 
without  feedback  about  the  consequences  of  its  fire; 

5.  Pire  froa  surviving  units  is  uniforaly  distributed 
over  the  area  in  which  eneay  forces  are  located; 

6.  Each  unit  presents  the  saae  vulnerable  area  to  eneay 
fire.  This  vulnerable  area  is  such  larger  than  the  effec¬ 
tive  area  of  a  single  round  of  eneay  fire.  Additionally, 
the  number  of  hits  required  for  a  kill  obeys  a  geometric 
probability  law. 

The  state  equation  that  relates  the  force  levels  and 
does  not  explicitly  contain  time  for  the  aiaed-fire  assump¬ 
tions  is  known  as  Lanchester's  square  Law, 
b(x/  -  x*)  =  a(i;2  -  y*)  . 
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The  state  equation  for  the  area  fire  model  yields 
Lanchester *s  Linear  Law, 

b  (Xo  -  x)  *  a  (y#  -  y) . 

Proa  these  state  equations  many  important  results  can  be 
deduced,  for  example,  battle  outcome  prediction  conditions. 
Further  analysis  of  Lanchester* s  classical  equations  and 
their  implications  will  not  be  done  here.  Suffice  it  to  say 
that  Lanchester* s  original  work  forms  the  basis  for  the  ana* 
lytical  combat  model,  with  his  simple  model  serving  as  a 
point  of  departure  for  the  more  complicated  and  elaborate 
analytical  models  discussed  later  in  this  thesis.  Readers 
of  this  thesis  and  students  who  intend  to  continue  with  work 
initiated  h9re  should  have  a  sound  understanding  of  Lanches¬ 
ter*  s  basic  models  before  continuing  further. 

Lanchester* s  equations  for  modern  warfare  and  his  equa¬ 
tions  for  area  fire  have,  over  the  years,  been  thoroughly 
dissected  and  analyzed  by  students  of  differential  equations 
and  students  of  analytical  combat  models.  One  such  student 
was  Gordon  Clark,  a  doctoral  student  at  Ohio  State  tfniver- 
sity  in  the  late  1960*s.  He  was  well  aware  of  the  power 
that  simulation  models  such  as  Carmonette,  a  high  resolution 
combat  model,  had  for  describing  the  complex  process  of  com¬ 
bat.  However,  the  major  drawback  to  simulation,  as  he 


perceived  it,  was  that  the  aethod  was  slow  and  expensive. 
Clark  proposed  that  the  use  of  an  additional  aodel  to  inter¬ 
pret  siaulation  results  would  greatly  benefit  Aray  studies 
which  required  nuaerous  runs  with  the  siaulation  aodel.  He 
went  on  to  develop  this  analytical  aodel  which  he  called  the 
Coabat  Analysis  Model  (COMAN) .  This  such  siapler  and  faster 
running  aodel  would  greatly  reduce  coaputer  execution  tiae 
and  costs,  and  could  serve  as  a  filter  aodel  for  Caraonette. 
A  siailar  approach  to  Clark's  methodology  in  developing 
COHAN  will  be  used  here  to  develop  the  analytical  aodel  for 
STAB.  The  next  chapter  details  the  first  steps  taken  to 
build  the  analytical  coabat  aodel  for  STAB. 
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III.  DAflS£A£  asait 

The  aggregated  force-on- force  analytical  model  is  an 
integral  part  of  combat  modelling  and  serves  several 
significant  purposes.  Bith  the  high  cost  of  computer  time 
and  the  limited  computer  core  storage  available  to  most  mil¬ 
itary  users,  the  analytical  combat  model  offers  a  viable 
alternative  to  detailed  high  resolution  simulations. 

Analytical  combat  models  can  be  used  in  a  stand  alone 
configuration  or  in  conjunction  with  a  high  resolution  simu¬ 
lation.  In  the  stand  alone  mode,  the  models  are  usually 
"transparent",  or,  in  other  words,  easy  to  understand.  They 
are  also  small,  quick  but  most  importantly,  they  increase 
the  user's  understanding  of  the  process  of  combat,  in  con¬ 
junction  with  high  resolution  combat  simulations,  the  ana¬ 
lytical  model  offers,  as  Gordon  Clark  states  in  his  doctoral 
thesis,  the  opportunity  to  "extrapolate  the  results  of  the 
combat  simulation  to  predict  the  outcome  from  force  mixes 
that  were  not  explicitly  simulated."  In  this  way  it  can  act 
as  a  filter  to  eliminate  certain  force  sixes  from  further 
investigation.  Also,  Clark  identifies  another  use  for  it  as 
"an  integral  component  of  a  large  unit  model  to  predict  the 
outcome  from  individual  unit  actions."  [Ref.  u] 
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&  aajor  objective  of  this  thesis  vas  to  build  a  siaple 
aggregated  coabat  model  >u  the  Apple  II  aicr ocoaputer.  The 
input  paraaeters  for  the  aodel  were  estimated  using  the  out¬ 
put  data  froa  a  STAR  run.  This  concept  of  estinating  the 
attrition  rate  coefficients  involved  additional  wori.  in 
developing  a  data  postprocessor  compatible  with  STAR ,  and 
significant  effort  in  the  statistical  analysis  of  the  infor- 
aation  provided  by  the  postprocessor. 

There  are  several  reasons  why  the  Apple  II  computer  was 
chosen  for  the  development  of  the  aodel  rather  than  a  main 
frame  computer.  First,  the  need  exists  for  a  'portable* 
aodel  that  is  easily  accessible  to  the  Aray  analyst.  A  desk: 
top  type  systea  that  can  run  a  model  quickly  and  with  mini¬ 
mum  setup  time  will  be  a  valuable  asset  to  the  analyst, 
flith  the  Apple  II  computer  gaining  wider  acceptance  in  the 
Aray  for  various  uses,  the  accessibility  of  a  model  on  such 
a  systea  is  greatly  increased. 

Second,  the  nature  of  the  model  is  such  that  it  lends 
itself  to  a  small  systea  foraat.  The  aggregated  form  of  the 
aodel  reduces  the  required  core  storage  froa  that  required 
by  a  high  resolution  aodel.  The  concept  of  DAMS I AC  is  to 
serve  as  a  filter  aodel  to  weed  out  unacceptable  force  sixes 
or  scenarios  before  turning  to  a  high  resolution  model.  The 


20 


idea  here  is  to  reduce  the  time  and  cost  of  running  numerous 
iterations  of  a  large,  high  resolution  model  on  a  aain  frame 
computer. 

Third,  the  cost  of  a  computer  system  to  support  the 
model  is  certainly  economical.  at  a  cost  of  under  S3000, 
the  apple  II  computer  system  with  a  model  such  as  DknSTkC 
provides  a  powerful  analytical  tool  for  a  minimal  cost,  cer¬ 
tainly  in  relation  to  the  cost  of  a  main  frame  computer. 

Pourth,  the  apple  II  has  the  capability  to  provide  addi¬ 
tional  analysis  functions  on  the  output  of  the  model.  Kith 
such  programs  as  apple  Plot  and  available  commercial  statis¬ 
tical  and  numerical  analysis  packages,  most  of  the  needed 
analytical  work  can  be  accomplished  on  the  apple  II  compu¬ 
ter.  This  can  result  in  additional  savings  in  time  by 
eliminating  the  requirement  to  transfer  any  data  from  one 
system  to  another. 

The  laurels  cast  on  the  apple  II  should  not  go  without 
some  discussion  of  its  limitations.  The  apple  II  is  limited 
in  capacity  to  64  kilobytes  of  Random  access  Memory  storage 
space.  So,  essentially  the  model  must  be  able  to  work 
within  the  amount  of  core  storage.  There  are  methods  of 
working  around  this  problem,  if  it  arises,  but  this  will 
usually  create  a  new  problem  or  limitation.  For  example. 
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soae  of  the  data  or  a  part  of  the  program  can  be  kept  on  the 
storage  diskette  and  put  into  computer  memory  only  when  it 
is  needed.  This  saves  space  in  the  computer's  memory  but 
increases  the  run  time  of  the  model  due  to  the  time  required 
to  access  the  diskette  to  retrieve  the  required  information. 

This  raises  another  limitation  of  the  ip  pie  II— its  run 
time.  Compared  to  a  main  frame  computer,  the  Apple  II  is 
very  slow.  This  is,  however,  relatively  speaking  since  com¬ 
puter  time  is  normally  measured  in  milliseconds  and  nanose¬ 
conds.  For  the  analyst  waiting  for  the  output  from  DAHSTAC, 
the  model  on  the  Apple  II  would  take  about  four  to  ten 
minutes  to  run.  This  is  probably  acceptable  for  most 
analysts. 

A.  SCENARIO 

The  model  portrays  a  battle  between  a  RED  tank  battalion 
attacking  a  BLOB  tank  company  in  the  10  z  10  KH  Fulda  Gap 
(Keyhole)  terrain.  See  Pigure  2  for  the  terrain  map. 

The  8L0E  defensive  positions  and  the  RED  positions  along 
with  their  attack  routes  (coordinating  points)  are  provided 
in  Table  1.  The  BED  units  attack  along  three  preplanned 
attack  routes  toward  the  BLUE  defensive  positions.  Grid 
coordinates  for  unit  locations  and  coordinating  points  were 
computed  as  the  centroid  of  the  elements  of  the  unit.  In 


22 


figure  2:  Terrain  3a? 
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Star's  data  base,  each  tank  has  its  own  individual  location. 
For  example,  BLUE  unit  1  has  four  different  locations  asso¬ 
ciated  with  the  it,  each  one  corresponding  to  a  tank  in  the 
unit.  The  centroid  of  the  unit  was  taken  to  be  the  average 
of  these  grid  coordinates.  Thus,  the  initial  location  for 
BLUE  unit  1  in  Table  1  is  the  average  of  these  coordinates: 

54150  100690  54150  100620 

54190  100520  54160  100790 

The  rate  of  RED  advance  is  a  constant  but  specified  by  the 
user.  BLUE  units  can  move  to  alternate  positions  at  a  rate 
designated  by  the  user.  However,  for  this  scenario,  BLUE 
units  remain  stationary  for  the  entire  battle. 

To  obtain  the  raw  output  data  that  would  correctly 
reflect  the  above  scenario,  it  was  necessary  to  "turn  off" 
many  of  the  submodels  of  STAB.  For  example,  the  Field 
Artillery  module  was  shut  off  by  scheduling  the  fir3t  fire 
mission  well  after  the  battle  ended,  weapon  systems  other 
than  tanks  within  the  fighting  units  were  rendered  useless 
by  setting  their  ammunition  counts  (basic  loads)  to  zero. 
Additionally,  the  communications  and  the  Air  Defense  modules 
were  not  played.  After  all  of  these  preliminary  actions 
were  taken,  the  result  was  a  battle  exclusively  between 
tanks.  This  effort  worked  well  to  reduce  a  complicated. 
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heterogeneous  force  battle  to  a  siaple  homogeneous  force 
battle.  One  minor  problem  arose,  however,  when  a  weapon 
system  other  than  a  tank  was  detected  and  fired  on  when,  by 
the  scenario,  it  did  not  even  exist.  The  effects  of  this 
target  acquisition  and  engagement  were  considered  insignifi¬ 
cant  and  were  disregarded. 


Table  1:  OMIT  LOCATIONS/COOHDINATING  POINTS 


BLOE 

X-COORD 

Y-COOBD 

SJZ  2 

REMARKS 

OMIT 

1 

54163 

100655 

4 

initial  location 

UNIT 

2 

54788 

98120 

4 

initial  location 

OMIT 

3 

57448 

97893 

4 

initial  location 

SI2 

X-COORD 

Y-COOBD 

si£S 

REMARKS 

ONIT 

1 

587  27 
52000 

101491 

96000 

10 

initial  location 
coordinating  pt 

UNIT 

2 

581  10 
52000 

101696 

96000 

10 

initial  location 
coordinating  pt 

OMIT 

3 

58727 

52000 

101491 

96000 

10 

initial  location 
coordinating  pt 

B.  ASSUMPTIONS 

The  attrition  rate  functional  form  chosen  for  the 
DAMSTAC  model,  which  will  be  discussed  later  in  this  thesis, 
requires  the  following  assumptions. 

1.  The  conditional  kill  rates, a  and  3  ,  are  constant  and 
independent  of  time  if  an  individual  firer  has  at  least 
one  acquired  target.  If  there  are  no  acquired  targets. 
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the  kill  rates  are  zero. 

2.  P  and  q,  the  probabilities  of  nondetection,  are 
assumed  to  be  constant  functions  of  tiae,  that  they  are 
independent  of  the  force  size  and  the  number  of  acquired 
targets,  and  that  the  probabilities  p  and  q  are  equal  for 
all  x-firer/y-target  and  y-f irer/x-target  combinations 
respectively. 

3.  Synergistic  effects  are  neglected,  therefore,  the 
actual  force  attrition  rate  is  equal  to  the  sum  of  the 
individual  kill  rates  for  the  opposing  forces. 

Additionally,  the  assumptions  listed  below  are  designed  to 
remove  the  distractions  of  a  complicated  scenario  and  to 
make  the  initial  model  easier  to  understand: 

1.  All  units  have  point  locations;  ie.  a  battalion  is 
represented  by  a  grid  coordinate  and  all  calculations  for 
the  battalion  are  made  from  that  location. 

2.  Units  are  homogeneous.  That  is,  they  are  composed  of 
identical  firing  elements/weapons  systems. 

3.  Fire  is  uniformly  distributed  over  surviving  enemy  tar¬ 
gets  within  range. 

h.  Sufficient  supplies  of  ammunition  for  both  sides  exist 
for  the  duration  of  the  battle. 

5.  Units  move  at  a  constant  rate  of  speed. 
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C.  SYNOPSIS  OP  THE  COMPOTES  MODEL 

The  DAMSTAC  computer  program  is  a  modular  type  program 
written  in  the  Applesoft  BASIC  language.  The  program 
accesses  external  data  files  as  a  source  of  input  of  initial 
data  for  the  two  forces.  These  data  files  are  produced  by 
the  user  with  the  help  of  the  Force  Maker  utility  program 
described  in  Appendix  C.  Additionally,  the  program  has  the 
capability  to  make  output  files  that  can  be  stored  on  a 
diskette  and  accessed  by  the  user  with  the  Results  Reader 
utility  program,  also  described  in  Appendix  C.  Results 
Reader  will  print  the  output  to  the  screen  or  to  a  printer. 

A  copy  of  the  program  for  the  model  is  included  as  List¬ 
ing  A  in  Appendix  F.  A  sample  run  of  the  model  including 
user  inputs  is  included  as  Figure  3.  Plow  charts  of  the 
main  program  and  subroutines  are  included  as  Figure  4,  Fig¬ 
ure  6,  and  Figure  7. 

1 .  In itialization  gf  Arrays  and  Variables 

Lines  30  thru  i70  dimension  the  arrays  and  matrices 
and  initialize  hardwired  data.  Dimensions  are  set  to  handle 
up  to  five  units  in  each  force  and  up  to  ten  coordinating 
points,  including  the  initial  location,  for  each  unit.  The 
user  can  modify  these  capacities  by  changing  the  dimensions 
of  the  appropriate  variables. 
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Figure  3: 


Sample  Bun  of  the  aodel 
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Lines  180  thru  350  read  BLOB  force  data  from  a  file 
to  include  the  number  of  units,  initial  locations,  any  coor¬ 
dinating  points,  and  rates  of  advance,  and  compute  total 
force  level  based  on  user  input  of  the  unit  force  levels. 

The  time  increment  (DT)  is  set  to  30  seconds  while  the  maxi- 
nua  time  of  the  battle  (HT)  is  set  to  1500  seconds.  These 
values  are,  of  course,  adjustable  by  changing  thea  in  line 
80.  For  both  forces,  the  breakpoint  criteria  is  set  at  0.3 
of  the  force  level  and  movement  criteria  is  set  at  0.5  of 
the  unit  strength,  even  though  the  movement  criteria  was  not 
used  for  any  runs  of  the  model. 

Lines  360  thru  520  perform  the  same  operations  for 
the  FED  forces  as  lines  180  thru  350  perform  for  the  BL0E 
forces. 

Lines  530  and  600  ask  if  the  printer  should  be 
turned  on  or  off  and  calls  the  printout  subroutine  to  print 
the  initial  data. 

2.  flaig,  program  (lines  61 0-780) 

The  main  program  performs  these  operations: 
increments  the  time, 

checks  if  the  time  exceeds  the  max  time, 

calls  the  movement  and  attrition  subroutines, 

recomputes  current  total  force  levels  after  attrition,  and 
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Figure  h:  Plow  Chart  of  the  Saia  Prograa 
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f 
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checks  force  levels  against  the  designated  breakpoints  for 
each  force,  to  determine  if  the  battle  has  ended. 

If  the  breakpoints  or  the  max  tine  are  exceeded,  the  battle 
terminates  and  the  program  prints  the  force  levels.  Addi¬ 
tionally  the  program  will  write  the  force  levels  for  each 
time  interval  onto  the  diskette,  (lines  1470  thru  1580) 

5.  aoveeent  Sa&ouiiaa  (Uasg  790-1J60) 

Both  BED  and  BLOB  forces  are  capable  of  moving. 
Movement  of  the  BIDE  force  is  nullified  by  setting  its  move¬ 
ment  rate  to  zero  in  the  initialization  part  of  the  program. 
The  movement  subroutine  uses  basic  geometry  to  compute  new 
coordinates  for  the  moving  unit  (see  Pigure  5) .  The  move¬ 
ment  process  is  based  cn  moving  a  unit  from  one  coordinating 
point  to  the  next  in  a  straight  line.  The  amount  of  move¬ 
ment  depends  on  the  rate  of  movement  (BB  (I)  or  BB  (I) )  and 


Figure  5:  Geometric  Representation  of  (Jnit  dovement 
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Figura  6:  Flo»  Chart  of  the  HOTeaent  Subroutine 


L 
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the  tine  interval  length  (DT) .  If  the  distance  computed  is 
farther  than  the  next  coordinating  point,  only  movement  to 
the  coordinating  point  is  accomplished.  Lines  800  to  970 
perform  computations  to  move  BED  forces  while  Lines  980  to 
1150  perform  the  same  computations  for  BLUE  forces. 

4.  Attrition  sjusiaaiiaa  UlQzli*. S) 

A  range  check  is  made  for  each  f irer/target  combina¬ 
tion  to  insure  the  units  are  within  firing  range.  If  so, 
attrition  caused  by  that  firer  on  the  target  is  coaputed. 
Sext  total  attrition  for  each  unit  is  computed  based  on  the 
partial  attritions  of  the  firers  as  a  fraction  of  their 
force  levels.  In  other  words,  when  a  firer  is  firing  at  3 
enemy  units,  each  target  unit  receives  1/3  of  the  firer *s 
fire  power. 

5.  Printout  Subrogtjpe  (lines  1490- 1530  &  lines 

1760-1940) 

There  are  two  different  output  formats  used  in  the 
program.  The  first  (lines  1490  to  1530)  prints  only  a  sum¬ 
mary  of  the  total  force  levels  for  both  forces  for  each  time 
interval.  This  is  the  same  data  that  is  saved  to  a  file  by 
the  final  portion  of  the  program.  The  second  format  (lines 
1760  to  1940)  prints  out  the  units  of  each  force  with  their 
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current  locations  and  strengths.  This  format  allows  the 
user  to  see  which  unit,  specifically,  is  being  attritted  and 
the  grid  coordinates  of  that  unit  at  the  time  of  attrition. 
6.  satgut  £lle  flaiiaa  Safesaalins  (lines  1540-17501 
This  subroutine  allows  the  user  to  sake  a  file  of 
the  output  from  the  current  run  of  the  program.  This 
routine  actually  produces  two  files,  one  consisting  of  the 
attrition  rate  coefficients  and  probabilities  of  non-acqui¬ 
sition,  and  the  second  containing  the  time  and  force  levels 
of  both  forces  at  that  time.  These  files  are  seperated  so 
that  the  second  file  can  be  used  with  the  APPLEPLOT  program 
to  produce  a  graph  of  the  two  force  levels  as  a  function  of 
time.  This  file  has  the  postscript,  " .BESOLTS",  added  to 
its  file  name,  while  the  first  file  with  the  attrition  rate 
coefficients  has  the  postscript  ".COEFS". 

D.  METHODOLOGY  FOR  COMPUTING  PARAMETER  ESTIMATES 

The  process  of  building  a  computer  simulation  combat 
model  requires  the  development  of  methodologies  for 
modelling  the  various  aspects  of  combat.  Constructing  a 
force-on-force  attrition  model  requires  methods  for  consoli¬ 
dating  combat  phenomena  into  simplified  algorithms  with 
coefficients  representing  one  or  more  facets  of  combat. 

This  section  discusses  one  significant  part  of  the 
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force-on-force  attrition  model:  the  methodology  for  estimat¬ 
ing  the  attrition  rate  coefficients  in  the  Lanchester  type 
attrition  rate  equations. 

In  deriving  the  equations  for  the  attrition  rates  (dx/dt 
and  dy/dt) ,  it  is  necessary  to  understand  the  processes  that 
must  occur  to  cause  a  firer  to  inflict  attrition.  These  key 
processes  can  be  summarized  as: 
the  environmental  process; 
the  target  acquisition  process; 
the  target  engagement  process. 

The  environmental  process  involves  the  terrain  and  weather 
factors  that  affect  firer/target  line  of  sight  (LOS)  .  The 
target  acquisition  process  represents  the  visible  capability 
of  the  firer  to  detect  a  target  on  the  battlefield,  while 
the  target  engagement  process  describes  the  firer’s  ability 
to  hit  and  kill  an  acquired  target.  These  three  processes 
will  all  be  represented  by  the  attrition  rate  coefficient, 

A,  in  the  attrition  rate  equation.  This  coefficient  will  be 
broken  down  into  its  three  components  which  describe  mathe¬ 
matically  the  phenomena  of  the  three  processes.  The  three 
components  are: 

Pv  *  the  probability  that  a  target  is  visible  to  the  firer 


(environmental  process) 


Pa  =  the  probability  that  a  visible  target  is  acquired 
(target  acquisition  process) 

<*  =  the  casualty  rate  at  which  a  firer  attrits  the 
acquired  target  (target  engagement  process) 

a  should  not  be  confused  with  k,  the  attrition  rate 
coefficient,  of  which  a  is  a  factor.  The  attrition  rate 
equation  may  be  represented  in  words  by: 

dx/dt  =  -(P(tgt  is  visible)][P(tgt  is  acquired)) (casualty  rate) 
The  second  probability  above  can  be  broken  down  further  to 
P(tgt  is  acquired  |  tgt  is  visible).  Note  that  this  is  a 
conditional  probability  since  it  is  assumed  that  the  target 
must  be  visible  to  be  acquired.  By  identifying  the 
mathematical  expressions  for  each  of  these  components  of  the 
attrition  rate  coefficient,  the  attrition  rate  equation  can 
be  obtained. 

1.  Tfre  5nv ironmental  Process 

There  are  generally  two  ways  to  represent  line  of 
sight  in  a  model: 

as  a  mathematical  simulation  of  actual  terrain; 
as  a  stochastic  process. 

The  mathematical  simulation  method  normally  involves  a 
digital  representation  of  terrain  elevation  over  a  segment 
of  the  terrain.  The  elevations  of  the  segments  between  the 
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firer  and  the  target  are  used  to  determine  if  line  of  sight 
exists.  Virtually  all  high  resolution  models  use  this 
method,  or  variations  of  it,  to  determine  line  of  sight. 

STAB  uses  a  different  form  of  mathematical  simulation  to 
represent  terrain.  This  is  called  functional  terrain,  in 
which  hills  and  forests  are  represented  by  elliptical  con¬ 
tours  of  elevation.  Hills  are  modelled  by  describing  them 
in  terms  of  bivariate  normal  equations  with  varying  parame¬ 
ters  for  each  hill.  Forests  or  built-up  areas  are  repre¬ 
sented  by  an  elliptical  pattern  on  the  ground  again  w.th 
varying  parameters  for  each  forest. 

For  a  low  resolution  or  highly  aggregated  model,  the 
mathematical  simulation  method  for  line  of  sight  is  usually 
not  practical  or  desireable  in  terms  of  core  space  and  com¬ 
putation  time  required.  Therefore,  the  second  method,  in 
which  LOS  is  represented  as  a  stochastic  process,  is  used. 
The  concept  of  this  process  is  as  follows.  The  effect  of 
terrain  on  the  process  of  line  of  sight  is  that  a  target 
falls  into  one  of  two  categories  or  states:  either  it  is 
visible  or  if  is  invisible. 

Assuming  a  steady-state  situation  for  the  terrain, 
probabilities  can  be  assigned  to  each  of  these  two  states. 
These  probabilities  represent  the  percentage  of  time  that 
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the  target  will  be  in  that  state  (risible  or  invisible)  for 
that  terrain  over  a  reasonable  aiotint  of  tiae.  Graphically 
the  states  of  the  target  can  be  represented  by  the  diagraa 
in  Pigure  8. 
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Figure  8:  Two-state  Stochastic  Graph  of  Inter  risibility 

Since  this  is  a  stochastic  process,  there  is  an 
exponential  distribution  associated  with  each  ox  the  two 
states.  This  distribution  represents  the  tiae  between  two 
occurrences  of  the  same  state.  On  the  graph  in  Pigure  8, 
the  time  between  the  end  ox  and  the  start  of  T*  is  expo¬ 

nentially  distributed  with  mean,  u .  Likewise,  the  tiae  bet 
ween  the  end  of  T*  and  the  start  of  T*  is  exponentially 
distributed  with  mean,  p  .  Graphically,  this  can  be  repre¬ 
sented  by  Pigure  9,  where  the  circles  represent  the  two 
states,  risibility  and  invisibility,  and  the  arrows 
represent  the  transition  from  one  state  to  the  other  at  the 


given  rate 


The  attrition  caused  by  a  firer  on  a  target  can 
occur  only  when  the  target  is  visible  (assuming  no  attrition 
from  wild,  blind  shots  or  indirect  fire)  .  Therefore,  the 
attrition  rate  is  a  function  of  the  probability  that  the 
target  is  visible  to  the  firer.  This  probability  will  be 
denoted  as  ,  and  can  be  expressed  as  the  proportion  of 
tiae  that  the  target  is  in  the  visible  state  in  a  steady 
state  situation,  p  *  -n  ^  . 

2.  Iks  tarask  agaaiaikkaa  ttasaag 

Assuaing  now  that  line  of  sight  exists  between  a 
firer  and  a  target,  the  process  of  target  acquisition  aust 
now  be  represented.  There  are  two  primary  modes  of  target 
acquisition: 

Serial  acquisition,  and 
Parallel  acquisition. 

Tn  the  serial  acquisition  process,  targets  are  acquired 
between  engageaents.  That  is,  in  concept,  the  firer  is 
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concantrat ing  on  the  engaged  target  and  does  not  or  is  not 
able  to  acquire  other  targets  on  the  battlefield.  Any  addi¬ 
tional  acquisitions  made  prior  to  engaging  the  current  tar¬ 
get  are  lost  and  must  be  re-acquired  after  the  current 
engagement  has  ended.  Exaaples  of  firers  who  use  serial 
acquisition  are: 

any  firer  who  is  "buttoned  up"  and  has  his  field  of  view 

narrowed  to  a  thin  band; 

wire  guided  missiles/anti-tank  systems. 

in  the  parallel  acquisition  mode,  the  firer  main¬ 
tains  his  list  of  acquired  targets  throughout  an  engagement 
(except  those  lost  through  loss  of  line  of  sight) .  Thus, 
when  an  engagement  has  ended,  by  the  target  being  killed  or 
loss  of  line  of  sight,  the  firer  may  immediately  engage  a 
target  that  was  previously  acquired  and  is  still  a  valid 
target.  All  systems  which  do  not  use  serial  acquisition  are 
assumed  to  use  the  parallel  mode. 

For  purposes  of  this  thesis,  all  firers  are  assumed 
to  use  parallel  acquisition.  Therefore,  serial  acquisition 
will  not  be  discussed  any  further.  The  reader  should  note 
that  the  serial  acquisition  process  does  have  validity  in 
modelling  action  on  the  battlefield,  and  should  be  consid¬ 
ered  in  future  enhancements  of  this  model. 


42 


To  explain  the  derivation  or  the  target  acquisition 
portion  of  the  attrition  rate  equation,  the  diagram  in  Fig¬ 
ure  9  must  be  aodified  to  that  in  Figure  10. 


INVISI8L£ 

Stats 


Pigure  10:  Three-state  Stochastic  Hodel  or  Intervisibility 


The  two-state  stochastic  process  has  been  changed  to 
a  three-state  process  by  dividing  the  visible  state  into  the 
visible-not  acquired  state  (VHA)  and  the  visible-acquired 
state  (Vi) .  ahat  has  occurred  is  that  a  target  can  aove 
from  the  visible-unacguired  state  to  the  visible  acquired 
state  at  a  rate  \  .  Since  both  of  these  states  coaprise  the 
visibility  3tate  in  total,  the  rate  of  moving  froa  either  of 
thee  to  the  invisibility  state  remains  at  |i  .  Notice  that 
there  is  no  direct  transition  froa  the  invisible  state  to 
the  visible- acquired  state.  This  is  only  a  matter  of 
representation,  since  it  is  assumed  that  no  more  than  one 
event,  ie.  moving  to  a  new  state,  will  occur  during  a  time 
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increment.  Thus,  the  situation  described  can  still  be 
represented  by  Figure  10. 

Another  point  of  interest  in  this  diagram  is  the 
fact  that  there  is  no  transition  from  risible-acquired  to 
visible-not  acquired.  This  implies  the  assumption  that  once 
a  target  is  acquired,  it  remains  acquired  unless  it  is 
killed  or  line  of  sight  is  lost,  in  which  case  it  transfers 
to  the  invisible  state. 

To  determine  the  effect  of  the  target  acquisition 
process  on  the  the  attrition  rate,  the  probability  that  a 
target  is  acquired  must  be  determined.  From  Figure  10,  and 
using  the  same  concept  of  a  stochastic  process  as  that  used 
in  section  III-D-1,  the  equation  can  be  written 
Pa  (tgt  acquired)  =  —  ^  *  -  -  .  (1) 

This  should  more  properly  be  stated  as  the  probability  that 
a  target  is  acquired  given  that  it  is  visible,  since  this  is 
actually  a  conditional  probability  with  the  state  of  acqui¬ 
sition  being  dependent  on  the  target  being  visible. 

Therefore,  the  probability  that  a  target  is  visible 
and  acquired  in  the  steady  state  is: 

P„a  =  P(tgt  vis  &  acq  ) 

Vd 

*  Pa  (tgt  acq  Jtgt  vis)  Pv(tgt  vis)  (2) 

“  (  \  +  n  )  (  V  +  n  ) 
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The  implications  of  this  equation  will  be  discussed  later  in 


Chapter  ?. 

3.  poraulatiso  the  AttritioQ-yate  Equation 

Using  equation  (2)  in  the  attrition- rate  equation, 
the  attrition  rate  of  a  single  firer  on  a  single  target  can 
be  summarized  using  expected  values  as: 

dx/dt  =  -Pva  a  (3) 


That  is,  the  attrition  rate  equals  the  rate  of  killing  an 
acquired  target  times  the  probability  that  a  target  is  visi¬ 
ble  and  acquired.  However,  since  there  are  likely  to  be 
more  than  one  firer  and  more  than  one  target  on  the  battle¬ 
field,  this  equation  must  be  modified  to  reflect  this. 

If  a  firer  is  going  to  cause  attrition  on  his  tar¬ 
gets,  he  must  be  able  to  acquire  one  or  more  of  these  tar¬ 
gets.  Por  firer  yl  with  x  targets,  the  probability  that  yl 
acquires  one  or  more  targets  is  just  one  minus  the  probabil¬ 
ity  that  he  acquires  no  targets.  Assuming  the  same 
probability  of  acquiring  any  of  the  targets,  the  probability 
of  acquiring  none  of  the  targets  is  (1  -  P„=)x  .  Therefore, 

v  a 


p 

a 


n 


(1  - 


?va^ 


0>> 
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The  attrition  rate  now  becomes 

dx/dt  =  -Pa  a  (5) 

for  a  single  firer  against  x  targets.  Since  there  are  more 
than  one  firer  on  the  battlefield,  each  assumed  to  be  acting 
independent  in  this  process,  the  attrition  rate  in  equation 
(5)  must  be  multiplied  by  the  number  of  firers(y)  and  the 
equation  is  now  complete: 

dx/dt  »  -P  a  y  .  (6) 

In  terms  of  the  probability  of  acquiring  a  target,  this 
equation  can  be  written  as: 

dx/dt  a  -(1  -  (1  -  Pva)X  )  ay  .  (7) 

This  is  the  same  equation  which  Clark  uses  in  his  COM&N 
model,  that  uses  results  from  a  high  resolution  model  to 
estimate  the  value  for  the  conditional  casualty  rate  and  the 
probability  that  a  target  is  not  acquired.  Henceforth,  (1  - 
P^)  will  be  referred  to  as  p  for  BED  firer/BLUE  target  com¬ 
binations  (or  q  for  a  BLUE  firer/BED  target  combination) . 

4.  £he  Target  Engagement  Process  (Bond§£  vs  Clark) 

The  target  engagement  process  can  be  modelled  in  its 
most  detailed  form,  as  in  the  high  resolution  model,  or  in 
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its  simplest  fora,  as  in  the  aggregated,  analytical  aodel. 
The  high  resolution  aodel,  such  as  STAB,  includes  explicit 
representation  of  every  factor  that  is  involved  in  the  pro¬ 
cess  of  aiaing  at,  firing  at,  hitting,  and  killing  a  target. 
These  factors  include: 

weapon  systea  characteristics; 
veapon/aaaunition  type  characteristics; 

target  characteristics  such  as  size,  shape,  vulnerability, 
etc. 

In  contrast,  the  low  resolution,  aggregated  aodel 
represents  the  engagement  process  by  one  value — the  attri¬ 
tion  rate  coefficient.  This  is  the  method  to  be  employed  in 
the  model  presented  here.  The  task  now  is  to  develop  a 
aethod  to  accurately  estimate  this  coefficient  for  use  in 
the  attrition  rate  equation.  Note  that  we  are  actually 
talking  about  two  attrition  rate  coefficients  and  two  equa¬ 
tions.  The  derivation  and  method  for  estimating  one  is  the 
same  as  the  other,  so  for  purposes  of  the  discussion,  refer¬ 
ences  to  an  attrition  rate  coefficient  will  be  in  the 
singular. 

There  are  presently  two  basic  approaches  to  the 
problem  of  estimating  attrition  rate  coefficients  for  paral¬ 
lel  acquisition  and  stochastic  line  of  sight.  The  first  is 
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the  use  of  the  independent  analytical  aodel,  involving 
weapon  system  characteristics  and  the  use  of  algorithas  to 
deteraine  the  expected  tiae  for  a  firer  to  kill  a  single 
target.  The  attrition  coefficients  are  then  the  reciprocal 
of  these  expected  tiaes.  This  aethod  will  be  referred  to  as 
Bonder's  aethod  after  Seth  Bonder  who  has  used  it  exten¬ 
sively  in  his  coabat  models  such  as  VECTOR  II. 

The  second  aethod  is  the  f it ted- parameter  analytical 
model,  involving  the  use  of  a  separate  high  resolution, 

Monte  Carlo  combat  simulation  to  produce  data  to  deteraine  a 
aaxiaum  likelihood  estiaate  of  the  attrition  coefficients. 
This  method  will  be  referred  to  as  Clark's  method,  since 
Gordon  Clark  popularized  it  in  his  doctoral  thesis  and  in 
the  development  of  the  Coabat  Analysis  Model  (COHAN) . 

Host  of  this  section  involves  a  discussion  of 
Clark's  aethod  and  the  practical  application  of  it.  How¬ 
ever,  in  order  to  provide  a  contrast  to  that  aethod,  a  short 
description  of  Bonder's  aethod  will  be  provided. 

In  general,  the  concept  of  Bonder's  aethod  is  to 
produce  an  algoritha  representing  the  expected  tiae  that  it 
takes  a  firer  to  kill  a  single  target.  The  eleaents  of  this 
algoritha  include  coaputations  for  line  of  sight,  target 
acquisition,  target  selection,  and  the  killing  process. 


Each  element  accounts  for  a  certain  anount  of  time  in  the 
expected  tiae  to  kill  a  target.  This  entire  process  is  seen 
by  Bonder  as  a  Harkov  renewal  process  where  the  outcoae  of  a 
round  fired  is  dependent  only  on  the  outcoae  of  the  previous 
round.  By  using  weapon  systems  characteristics  data  as  raw 
input,  tiaes  to  acquire,  select,  and  kill  a  target  can  be 
determined  and  used  in  the  algoritha  for  expected  tiae  to 
kill  a  target. 

Bonder *s  method  uses  an  independent,  analytical 
aodel  requiring  a  large  data  base  of  factors,  such  as  range 
dependent  weapon  system  characteristics,  quantifications  of 
line  of  sight  and  target  acquisition  parameters,  to  support 
the  computation  of  the  attrition  rate  coefficient. 

A  detailed  explanation  of  Bonder's  method  for  esti¬ 
mating  the  attrition  rate  coefficients  can  be  found  in  the 
basic  documentation  for  the  VECTOB  II  model  [fief.  5].  Basi¬ 


cally,  Bonder  says  that  the  attrition  rate  coefficient  ^ is 
represented  in  the  parallel  acquisition  node  by  A  ^  *  Qij  af 
where  Q^is  the  probability  that  a  given  group-i  weapon  is 


firing  at  a  group- j  target. 


can  be  expressed  in  terms 


of  the  probability  that  a  group-j  target  is  visible  and  has 


been  detected  by  a  grcup-i  firer.  This  probabiltiy  is 
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Qij.  -  s./ri  (1  -  s.,)  for  j  >  2  . 

L3  !3  j«i  i3 

Assuming  homogeneous  forces,  this  equation  reduces  to 

Q.  .  =  S.  .  ,  and  the  i-1  subscripts  can  be  eliminated 
13  13 

since  there  is  no  distinction  among  types  of  force  elements. 
His  resulting  expression  for  S  is  identical  to  that  of  equa¬ 
tion  4  ,  where  s,.-and  P_  are  defined  the  same. 


13 


Bonder* s  expression  for  the  conditional  casualty 
rate,  a  •  is  written  as  the  inverse  of  the  expected  time  to 
kill  an  acquired  target.  This  expression,  when  further  bro¬ 
ken  down  contains  elements  dependent  on  the  weapon  system 
characteristics ,  target  characteristics,  and  a  number  of 
other  subsystem  performance  parameters.  These  parameters 
were  determined  empirically  outside  the  operational  environ¬ 
ment  of  the  battlefield  and  molded  together  in  this  pre¬ 
determined  structure  to  produce  the  value  of  the  attrition 
coefficient.  This  derivation  of  the  attrition  coefficient 
is  the  point  where  Bonder's  method  and  Clark's  method 
differ. 


Clark's  methodology  is  based  on  the  assumption  that 
the  process  of  attrition  in  combat  is  a  nonhomogeneous,  two 
independent  type  Poisson  Process  where  the  time  between 
casualties  has  an  exponential  distribution  with  parameter  \  , 
called  the  total  conditional  casualty  rate.  Likewise,  the 
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time  between  BED  casualties  and  the  tine  between  BLOE 
casualties  each  has  an  exponential  distribution  with  parame- 
ters  \r  and  respectively ,  where  X  =Xr*Xb.  From  the 
theory  of  Poisson  Processes,  the  probability  of  a  BED  casu¬ 
alty  occurring  is*r/x  ,  and  the  probability  of  a  BLOE  casu¬ 
alty  occurring  is^b/^  •  With  this  basis,  Clark  is  able  to 
develop  a  likelihood  functon  in  terns  of  the  conditional 
casualty  rates  a  and  0  ,  of  which  the  casualty  rates  are 
functions,  ie.,  Ab=  f  (a)  and  Ar=  f  (£) .  [Bef.  6]  It  should 
be  noted  that,  in  order  to  aaintain  consistency  of  terminol¬ 
ogy,  the  equations  attributable  to  Clark's  work  have  been 
adjusted  to  reflect  the  definitions  presented  in  this  the¬ 
sis. 

Clark  also  uses  the  concept  that  attrition  is  a 
function  of  target  acquisition,  which  depends  on  the  prob¬ 
ability  of  a  firer  acquiring  a  target.  This  concept  pro¬ 
duces  two  more  parameters  which  must  be  estimated:  p,  the 
probability  that  a  BLUE  target  is  not  acquired  by  a  BED 
firer;  q,  the  probability  that  a  BED  target  is  not  acquired 
by  a  BLUE  firer. 

For  the  functional  form  in  Clark's  COilAN  model,  the 
equations  can  be  written  as: 
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(8) 


-a(i-p"l  with  m(o)  =  m9 
fr*  */3[l-  qn  ]  with  n(o)=  n<> 
where  the  parameters  are  defined  by: 

the  conditional  attrition  rate  function  for  BLOE 

forces 

fP=  the  conditional  attrition  rate  function  for  BED  forces 
<*=  conditional  kill  rate  for  a  BED  weapon  against  acquired 
BLOB  targets 

0-  BLOE  weapon  conditional  kill  rate  corresponding  to  a 

p  =  probability  that  a  BLOE  target  is  unacquired  by  an 
individual  RED  firer 

q  *  probability  that  a  BED  target  is  unacquired  by  an 
individual  BLOE  firer 

a  *  nuaber  of  BLOE  firers  surviving 

n  *  nuaber  of  RED  firers  surviving 

These  functional  forms  depend  directly  on  the  abil¬ 
ity  of  firers  to  detect  targets.  To  see  this,  consider  just 
the  conditional  attrition  rate  function  for  the  BLOE  force, 

.  The  definition  of  p  implies  that  the  value  of  (1-pm)  is 
the  probability  that  a  BED  firer  has  at  least  one  acquired 
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BLUE  target.  Therefore,  the  attrition  rate  for  the  BLUE 


force  is  the  Itill  rate  for  acquired  targets  times  the  prob¬ 
ability  of  having  an  acquired  target  times  the  number  of 
firers. 

It  should  be  noted  that  Lanchester's  classical 
attrition  rate  functional  forms  are  special  cases  of  Clarkes 
attrition  rate  functional  forms,  and,  under  appropriate  con¬ 
ditions,  Clark's  equations  reduce  to  Lanchester's  equations. 
When  p  and  q  are  zero,  or  very  nearly  zero, 

fb  *  -  a  (1-0m)  n  =  -  a  n  (9) 


and. 


fr  *  -  0  ( 1-0n)  m  =  -  pm  (10) 

which  are  Lanchester's  classical  aimed  fire  equations.  When 
individual  targets  become  increasingly  difficult  to  acquire, 
p  and  q  assume  values  close  to  one  and  it  can  be  proven  that 


in  the  limit: 


at  ( 1-p)  an 


(11) 


fr  a  -  3  (1-q)  an 


If  the  attrition  rate  coefficient  a  is  set  equal  to  a(1  - 
p)  and  the  coefficient  b  is  set  equal  to  3(1  -  q),  then  the 
two  equations  above  are  identical  to  Lanchester's  equations 
for  area  fire. 

In  his  doctoral  dissertation,  Clark  developed  the 
likelihood  function  fox  a  heterogeneous  aodel,  however  since 
the  O&MSTkC  aodel  assumes  homogeneous  forces,  the  necessary 
equations  have  been  reduced  to  a  simpler  form  than  that 
derived  by  Clark.  By  differentiating  the  logarithm  of  the 
likelihood  function  with  respect  to  p,  the  following  equa¬ 
tion  is  derived: 


**  &  <:  -  L) 


1  V.  >  *  1 

•V.*» 


where 


Cj.  *  0  if  casualty  k  is  8ED 


5h 


1  if  casualty  k  is  BLOB 


i  *  number  of  BLOB  survivors  prior  to  kth  casualty 
k-1 

^  *  nusber  of  HBD  survivors  prior  to  kth  casualty 


t,  -  t.  ,  »  tiae  between  casualty  k  5  casualty  k-1 

K  »  1C—  l 


Likewise,  by  differentiating  the  logarithm  of  the  likelihood 
function  with  respect  to  a  ,  and  by  solving  for  a  ,  a  second 
equation  in  teras  of  £  and  &  is  produced: 


Differentiating  with  respect  to  q  and  0  produces  equivalent 

A  A 

type  equations  which  will  allow  the  estimates  for  q  and  &  to 
be  deterained  in  the  same  manner  as  follows. 

Generally,  the  two  equations  (13)  and  (iu)  can  now 
be  solved  simultaneously  to  deteraine  the  values  for  p 
and  a  .  Unfortunately,  this  is  easier  said  than  done.  The 
complexity  of  these  equations  (aore  so  in  the  heterogeneous 
case)  does  not  allow  explicit  equations  for  p  and  a  .  The 
remainder  of  this  chapter  will  discuss  a  method  for  deter- 
aining  the  estimates  of  p  and  a  from  these  equations  and  a 
way  to  extract  the  necessary  data  from  a  high  resolution 
Sonte  Carlo  siaulation  run  to  compute  these  estimates. 
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E.  ESTIMATING  PARAMETERS  FROM  A  HIGH  RESOLUTION  MODEL 
As  mentioned  before,  equations  (13)  and  (14)  must  be 
solved  simultaneously  to  produce  values  for  a  and  p.  Gener¬ 
ally,  these  equations  are  too  complicated  to  produce  such 
values  easily  so  other  techniques  aust  be  used.  The  aethod 
suggested  by  Clark  is  to  enter  a  value  for  p  in  equation 
(14)  and  solve  for  a.  This  value  and  the  selected  value 
for  p  are  then  put  in  equation  (13)  and  the  equation  is 
evaluated.  Notice  that  both  values  are  estimates  since  we 
are  coaputing  the  maximum  likelihood  estimates  of  p  and  a  . 
The  values  of  p  and  a  which  give  a  value  for  equation  (13) 
closest  to  zero  are  the  best  maxiaua  likelihood  estimates 
for  p  and  a  .  If  equation  (13)  equalled  zero  exactly,  the 
true  values  of  the  maximum  likelihood  estimates  would  be 
known. 

In  using  these  equations,  the  input  received  from  the 
high  resolution,  Monte  Carlo  simulation  must  include: 
the  time  between  casualties; 

the  number  of  RED  and  BLUE  survivors  prior  to  the  kth 
casualty  ; 

the  total  number  of  casualties  for  RED  and  BLUE  for  the 
battle ; 

the  type  of  casualty  (RED  or  BLUE)  of  the  kth  casualty. 


In  Clark's  theory  for  a  homogeneous  force  battle,  the 
conditional  attrition  rates  and  the  probabilities  that  a 
target  is  not  acquired  are  assuaed  to  be  constant  over  the 
entire  battle.  This  aeans  that  equations  (13)  and  (14)  are 
summed  over  all  casualties  in  the  battle  and  the  estiaates 
for  p  and  a  ,  which  are  determined  from  thea,  are  also 
assuaed  to  be  constant  over  the  entire  battle. 

In  applying  this  theory  to  STAB  output,  the  authors 
obtained  questionable  results  for  the  values  of  the  parame¬ 
ter  estimates.  The  estimates,  p  and  q,  were  computed  to 
have  values  of  0.99  and  0.93  respectively.  These  values 
imply  that  both  BLOB  and  RED  forces  experienced  difficulty 
in  acquiring  targets.  Analysis  of  the  Fulda  Gap  terrain  and 
consideration  of  the  respective  locations  of  the  opposing 
forces  support  the  estimated  values  for  p  and  q.  However, 
the  stated  value  for  p,  0.99,  is  a  truncation  of  the  actual 
value  computed.  If  proper  rounding  procedures  had  been 
used,  $  would  have  been  set  equal  to  1.0,  which  aeans  that 
the  RED  forces  would  not  detect  any  BLUE  targets.  Nondetec¬ 
tion  implies  that  no  BLUE  attrition  will  occur.  But  this 
contradicts  the  STAR  battle  which  did,  in  fact,  produce  BLUE 
casualties  amounting  to  about  two-thirds  of  the  total  BLUE 
force. 


57 


The  estimates.  Sand  0,  were  competed  to  be  2.4  and  0.03, 
respectively.  The  relatively  large  value  of  a  implies  that 
when  the  RED  forces  are  able  to  detect  BLUE  targets,  the 
BLUE  force  will  suffer  heavy  casualties.  There  appears  to 
be  an  unusually  large  difference  between  the  values  for  a 
and  $  which  are  not  reflected  by  the  values  for  p  and  g. 

This  disparity  is  most  likely  attributable  to  the  restric¬ 
tion  that  all  parameter  values  be  held  constant  over  the 
entire  battle.  Experimentation  with  different  parameter 
values  in  DANSTAC  indicates  that  a  value  of  about  0.40  more 
closely  approximates  the  attrition  observed  in  the  STAB  bat¬ 
tle.  Experience  with  these  computations  reveals  that  this 
is  not  necessarily  the  case  and  that  a  more  likely  assump¬ 
tion  is  that  the  parameters  are  constant  over  different  time 
intervals  of  the  battle.  In  general,  a  simple  battle  can  be 
divided  into  phases  where  the  intensity  of  the  battle  is 
constant  within  each  phase  but  different  among  the  phases. 
Basically,  the  phases  are: 

movement  to  contact  (low  intensity) 
initial  contact  (mid  intensity) 
main  battle  (high  intensity) 

withdrawal  from  contact  or  annihilation  (mid/low  inten¬ 
sity) 
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Moving  from  one  phase  of  the  battle  to  the  next  should  cause 
a  change  in  conditional  casualty  rates  (  aand/3  )  as  well  as 
the  probability  of  not  acquiring  a  target  (p  and  g) . 

The  question  now  becoaes  how  to  deteraine  the  breakdown 
of  the  battle  by  phases  using  the  output  froa  a  high  resolu¬ 
tion,  Monte  Carlo  siaulation.  That  is,  where  do  the  changes 
in  phases  occur  in  the  siaulation?  One  possible  aethod  is 
to  look  at  the  nuaber  of  casualties  as  a  function  of  tine. 
Another  is  to  look  at  the  nuaber  of  shots  fired  as  a  func¬ 
tion  of  time  (see  Pigure  11) .  In  both  cases,  time  is  seg¬ 
mented  into  time  intervals  (eg.  every  100  seconds)  and  a 
histograa  of  the  casualties  or  shots  is  plotted  for  each  of 
these  intervals.  Having  done  this  for  several  battles,  the 
results  identify  fairly  clear  choices  for  the  end  of  one 
phase  of  the  battle  and  the  beginning  of  the  next  phase. 
However,  more  analysis  of  these  two  methods  reveals  that 
they  might  not  truly  identify  the  phases  of  the  battle  and 
that  an  alternate  method  should  be  considered,  primarily  one 
based  on  range.  This  concept  is  discussed  in  detail  in 
Chapter  VI. 

Using  the  nuaber  of  casualties  per  time  interval,  a  dis¬ 
torted  picture  of  the  intensity  of  the  battle  might  be 
obtained.  The  interaction  of  the  probability  of  acquiring  a 
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?igura  11:  Histograa  of  Shots  par  Tine  Interval 

target  with  the  attrition  rata  of  a  firer  light  offset  each 
other  and  produce  results  which  conflict  with  the  true  bat¬ 
tle  intensity.  In  ether  words,  if  the  acquisition  probabil¬ 
ity  is  high  but  the  attrition  rate  is  low,  the  nuaber  of 
casualties  in  a  tiae  interval  light  indicate  a  low  intensity 
phase  of  the  battle  when,  in  fact,  it  is  a  high  intensity 
phase.  Therefore,  this  lethod  should  be  used  with  caution, 
if  at  all.  The  aethod  of  using  the  nuaber  of  shots  per  tine 
interval  has  siailar  faults.  Twenty  shots  by  one  firer  in  a 
tiae  interval  is  auch  different  in  teras  of  intensity  than 
one  shot  by  twenty  firers.  Again,  this  aethod  should  be 
used  with  caution. 
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The  question  is  then:  what  method  should  be  used?  One 
method,  which  is  used  by  the  Combined-  Arms  Combat  Develop¬ 
ment  Activity  (CACDA/CASAA)  ,  is  measuring  the  number  of 
shots  per  firer  per  time  interval.1  Although  it  also  has 
some  faults,  it  does  provide  a  better  indication  of  the 
intensity  of  the  battle.  It  is,  for  the  most  part,  indepen¬ 
dent  of  the  conditional  casualty  rates,  although  the  number 
of  shots  would  indicate  a  dependency  on  target  acquisition. 
This  dependency  is  reduced  by  averaging  the  number  of  shots 
over  the  number  of  firers.  In  effect,  this  measure  can  be 
considered  to  be  unbiased  by  the  changes  in  conditional 
casualty  rates  and  probabilities  of  not  acquiring  targets. 
Thus  one  is  able  to  get  a  feel  for  the  intensity  of  the  bat¬ 
tle  from  the  histogram  of  the  number  of  shots  per  firer  per 
tike  interval.  This  is  not  necessarily  scientific  reason¬ 
ing,  but  if  seems  to  work. 

Once  the  phases  of  the  battle  are  delineated  for  a  high 
resolution,  Monte  Carlo  simulation  run,  the  next  step  is  to 
determine  how  to  treat  the  'dead*  time  from  the  beginning  of 
a  phase  to  the  first  casualty  and  from  the  time  of  the  last 
casualty  to  the  end  of  the  phase.  For  phases  in  the  middle 


^Conversation  with  Mr.  Ernest  Boehner  from  CASAA  on  5 
December  1981 
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of  the  battle,  each  of  these  two  tines  is  a  portion  of  the 
tine  between  two  casualties  and  should  thus  have  a  bearing 
on  our  maximum  likelihood  estimates  from  equations  (1 3)  and 
(1«)  • 

To  resolve  this  situation,  the  time  of  the  last  casualty 
within  a  phase  is  used  as  the  end  of  the  phase.  For  exam¬ 
ple,  if  a  phase  is  supposed  to  end  at  1200  seconds  and  casu¬ 
alties  occur  at  1140  and  1225  seconds,  the  phase  will  be 
terminated  at  1140  and  the  new  phase  begun  at  1141.  This 
also  allows  this  time  between  casualties  to  be  accounted  for 
in  the  next  phase.  This  method  produces  some  variation  in 
phase  lengths  when  running  several  replications  of  the  bat¬ 
tle,  but  this  variation  is  minimal,  especially  since  the 
beginning  and  end  of  phases  are  being  judged  from  something 
less  than  a  scientific  method. 

Several  replications,  as  many  as  affordable,  should  be 
used  to  compute  the  maximum  likelihood  estimates  for  the 
four  parameters.  By  using  a  large  sample  population  (the 
casualties),  a  more  accurate  estimate  of  the  parameter  can 
be  obtained.  Likewise,  when  the  heterogeneous  model  is 
used,  a  larger  sample  population  is  desirable,  since  the 
number  of  engagements  between  two  particular  weapons  systems 
will  be  much  less  than  the  number  of  engagements  in  the 
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homogeneous  case.  This  could  produce  values  that  are 
statistically  insignificant. 

k  postprocessor  to  the  high  resolution  model  vas  devel¬ 
oped  to  make  the  numerous  calculations  required  to  produce 
these  estimators.  The  postprocessor  uses  an  optimization 
routine  to  determine  the  maximum  likelihood  estimates  from 
the  given  equations.  Detailed  information  on  the  postpro¬ 
cessor  is  provided  in  Appendix  E. 


17.  ALTERNATE  ED NCTIOHAL  FQBMS 


Different  »ilitary  situations  have  been  hypothesized  to 
yield  different  functional  for as  for  Lanchester-type  aqua¬ 
tions.  The  term  "Lanchester-type  equations"  refers  to  any 
differential-equation  model  of  coabat  and  the  term  "func¬ 
tional  fora"  relates  to  the  structure  of  the  attrition  rate, 
i.e.,  how  the  attrition  rate  equation  is  written  and  which 
parameters  are  considered  explicitly  in  the  equation. 
Henceforth,  a  standard  shorthand  notation  will  be  adopted  to 
refer  to  the  functional  forms  used  in  homogeneous-force  Lan¬ 
chester-type  coabat  models.  For  example,  Lanchester's  equa¬ 
tions  of  modern  warfare  are  expressed  as  P|F  functional 
forms  (or,  simply  F|F  attrition)  because  they  are  propor¬ 
tional  to  only  the  number  of  enemy  firers.  Taylor  [Ref.  7] 
provides  a  good  explanation  of  this  convenient  shorthand 
notation. 

Generally,  the  classical  Lanchester  functional  forms  for 
expressing  the  attrition  rate  for  forces  using  aimed  fire  or 
area  fire  are  too  simple  to  be  used  "as  is"  in  a  detailed 
analytical  combat  model.  The  assumptions  required  for  these 
basic  equations  to  hold  are  too  limiting  and  do  not  allow 
for  all  of  the  processes  that  are  known  to  occur  in  coabat 
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to  be  considered  in  the  force  attrition  model.  Some  of  the 
■ore  coaaon  functional  foras  for  attrition  rates  that  have 
been  considered  and  used  in  Lanchester-coabat  models  are 
shove  in  Table  2. 

◦f  special  interest  is  the  T|T  functional  fora.  This 
functional  fora  is  particularly  useful  in  describing  attri¬ 
tion  in  target  rich  environments.  Peterson  [Ref.  8]  hypoth¬ 
esized  that  the  equations, 

dx/dt  =  -ax  and  dy/dt  =  -by, 

characterize  the  early  stages  of  a  small  unit  engagement  in 
which  the  vulnerability  of  a  force  dominates  its  ability  to 
acquire  enemy  targets.  It  is  recommended  that  Peterson's 
T|T  functional  form  be  considered  in  the  heterogeneous  model 
described  in  Chapter  71  to  model  attrition  in  the  initial 
phase  of  the  battle.  The  remaining  sections  of  this  chapter 
vill  discuss  several  not-so-comaon  functional  forms  that 
have  been  formulated  and  used  in  "recently"  developed  combat 
models. 

A.  HELHBOLD-TYPE  EQUATIONS 

An  alternative  functional  fora  for  homogeneous-force 
attrition  rates  was  proposed  by  R.  Helabold.  Helmbold 
believed  that  when  opposing  force  sizes  were  grossly  une¬ 
qual,  the  larger  force  experienced  "inefficiencies  of  scale" 


Table  2:  FUNCTIONAL  FORM  NOTATION 


FUNCTIONAL  FORM  DIFFERENTIAL  EQUATION 


F  |  F 


dx/dt  *  -ay 
dy/dt  =  -bx 


FTj  FT 


dx/dt  =  -axy 
dy/dt  *  -bxy 


F|  FT 


dx/dt  =  -ay 
dy/dt  =  -bxy 


T  |  T 


dx/dt  =  -ax 
dy/dt  =  -by 


(F+T)  |  (F+T) 


dx/dt  =  -ay  -  cx 
dy/dt  =  -bx  d y 


in  producing  casualties.  Limitations  on  space  available  for 
maneuver,  reduced  reaction  tine,  increased  command,  control, 
and  communication  problems,  etc.  result  in  the  reduction  of 
casualty  producing  effectiveness.  These  limitations  may 
well  prevent  the  larger  force  from  using  its  full  destruc¬ 
tive  capability  against  the  smaller  force.  There  is  no 
definition  for  "grossly”  unequal  force  sizes,  but  the  fol¬ 
lowing  hypothesis  is  proposed.  The  validity  of  using  Helm- 
bold’s  equations  in  a  combat  model  is  directly  related  not 
only  to  the  force  ratios  (i. e. ,  x/y  or  y/x)  but  to  the  rela¬ 
tive  size  of  the  twc  opposing  forces  with  respect  to  their 
surrounding  terrain.  For  example,  a  5  to  1  force  ratio  may 
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well  be  a  situation  where  Helabold's  equations  apply  if  both 
of  the  opposing  force  sizes  are  relatively  large  with  res¬ 
pect  to  the  terrain  available  to  then  for  movement,  conceal¬ 
ment,  etc.  Conversely,  a  force  with  a  30  to  1  force  ratio 
aay  not  exhibit  inef f iciences  of  scale  if  the  opposing  force 
sizes  are  saall  in  relative  sizes  and  with  respect  to  their 
surrounding  terrain.  Consequently  the  decision  to  use  Helm- 
bold's  functional  fora  for  attrition  rates  rather  than  some 
other  fora  (like  Lanchester's  equation  for  modern  combat) 
really  depends  on  the  tactical  "situation". 

Helabold's  functional  form  introduces  a  modification 
that  accounts  for  the  reduced  fire  effectiveness  of  the  lar¬ 
ger  force.  In  mathematical  fora  his  equations  would  read: 
dx/dt  =  -a  g  (x/y)  y  with  x(Q)  =  x* 

dy/dt  =  -b  h  (y/x)  x  with  y(0)  *  y0 

where  a  and  b  represent  time-dependent  attrition-rate  coef¬ 
ficients,  and  g  (x/y)  and  h(y/x)  represent  functions  than 
modify  fire  effectiveness  of  individual  x  and  y  firers 
respectively.  According  to  Helabold's  hypothesis  the  two 
functions  g(x/y)  and  h  (y/x)  must  satisfy: 
g(1)  =  h(1)  *  1 
g(u)  *  h(u)  =  E(u) 

E(u)  is  a  strictly  increasing  function  of  u. 
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So,  the  generalized  Helabold-type  conbat  equations  can  be 
written  as: 

dx/dt  =  -a  E(x/y)  y 
dy/dt  =  -b  E(y/x)  x 

Por  the  case  in  which  E(u)  *  u  with  c  a  constant 
greater  than  or  equal  to  zero,  the  following  equations  for 
Helabold-type  coabat  are  obtained: 


dx/dt  =  -a 

|«w 

(x/y)'  y 

with 

x  (0) 

dy/dt  =  -b 

(y/x)l“w  x 

with 

y  (0) 

W  is  coanonly  called  the  "Weiss  parameter*'  with  w  =  1-c. 

The  above  equations  represent  the  special  exponential 
Case  of  Helobold  coabat.  If  a(t)  and  b(t>  represent  constant 
attrition  rate  coefficients  denoted  by  a  and  b  respectively, 
then,  when  W  =  1,  the  above  equations  reduce  to  Lanchester's 
Equations  for  Modern  Combat.  When  W  =  1/2  Lanchester's 
Area-fire  Equations  are  obtained. 

B.  HOW  TO  CHOOSE  THE  RIGHT  FUNCTIONAL  FORK 

A  very  logical  question  that  readers  of  this  study  may 
ask  at  this  time  is,  "Now  that  I've  seen  some  of  the  Lan- 
chester-type  attrition  rate  functional  forms,  which  one 
should  I  choose  for  ay  analytical  model?"  The  answer  to 
this  question  is  that  there  is  no  answer  to  this  question. 
The  choice  of  functional  form  is  purely  arbitrary  and  is 
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usually  left  to  the  discretion  of  the  aodel  builder.  The 
following  steps,  however,  aay  help  in  this  choice. 

1.  Carefully  consider  the  scenario  that  the  coebat  aodel 
is  intended  to  portray.  The  scenario  aay  well  indicate 
that  soae  particular  subset  of  the  functional  foras  is 
preferable. 

2.  List  the  assuaptions  that  aust  be  Bade  for  each  of  the 
functional  foras  in  the  subset  to  hold  and  see  which 
assuaptions  are  acceptable  to  you  and  your  aodel  and  which 
are  not. 

3.  observe  the  paraaeters  involved  in  each  functional 
fora  and  eliminate  those  foras  for  which  the  parameter 
values  or  estiaates  are  unobtainable  or  for  which  you  are 
not  willing  to  plug  in  arbitrary  values. 

Az  this  poin*,  the  final  selection  of  functional  fors?  •-? 
purely  one  of  subjective  judgement. 

There  is  no  requirement,  however,  that  a  model  builder 
must  accept  one  of  the  already  developed  functional  forms. 
There  is  no  doubt  that  other  foras  can  be  formulated,  but 
the  formulation  will  require  extensive  research,  analysis 
and  work.  A  good  check  to  see  if  an  attrition  rate  equation 
is  reasonable  is  to  insure  that  the  Lanchester  classical 
equations  are  special  cases  of  the  new  functional  form. 
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Once  the  final  functional  form  has  been  selected  or 
formulated,  it  may  be  enhanced  as  necessary  to  reflect  more 
of  the  complexity  and  realism  of  coabat.  Some  additional 
operational  factors  that  can  be  considered  are: 

(Jnit  breakpoints  and  battle  termination  criteria; 

Suppress ion ; 

C3; 

Reinforcements  or  r eplaceaents. 

The  main  assumption  that  allows  these  operational  factors  to 
be  "added"  to  the  attrition  rate  equation  is  that  synergis¬ 
tic  effects  are  not  considered  in  Lanchester-type  force-on- 
force  coabat  models.  No  synergistic  effects  imply  that  the 
Principle  of  superposition  holds,  therefore,  the  effects  of 
the  operational  factors  are  independent  and  additive. 


?.  FDTUBE  gMflAEC&flgflXS 

There  are  auaerous  enhanceaents  to  the  analytical  aodel 
that  can  be  aade.  The  enhanceaents  that  are  recoaaended, 
for  the  aost  part,  are  not  ained  at  iaproving  the  inner 
workings  of  the  current  aodel,  but  at  building  successive 
aodels  that  are  sore  detailed,  require  fewer  liaiting 
assuaptions,  and  are  generally  aore  coaplicated. 

For  follow-on  aodels  that  retain  the  hoaogeneous  force 
concept,  it  is  recommended  that  the  aodel  have  the  ability 
to  allow  the  user  to  choose  alternative  functional  foras  for 
the  force  attrition  rates.  This  should  be  done  to  better 
describe  the  situation  of  the  battle  at  a  particular  tine. 

A  likely  nethod  would  be  to  change  froa  one  functional  fora 
to  another  during  the  battle  in  order  to  more  realistically 
portray  the  attrition  occurring  then.  This  enhancement 
would  require  considerable  effort  in  developing  the  maximum 
likelihood  estimates  for  the  various  parameters  used  in  each 
alternative  functional  fora,  additional  work  with  the  STAfi 
postprocessor,  and  significant  change  to  the  current  DAMSTAC 
aodel.  Questions  which  arise  froa  this  suggestion  are: 

Which  forn(s)  should  be  chosen? 
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How  do  you  know  when  to  change  functional  fora? 

These  questions  are  partially  answered  in  Section  I?-C. 
Further  research  into  this  area  is  required. 

The  movement  subroutine  of  the  aodel  needs  to  be  modi' 
fied  to  aore  realistically  portray  the  movement  of  units  on 
the  battlefield.  Sven  though  a  relatively  low-resolution 
time-step  operation  is  occurring,  the  resolution  of  the 
movement  of  units  can  be  refined  to  vary  the  rate  of  move¬ 
ment  froa  one  time  step  to  the  next,  if  necessary.  Cur¬ 
rently,  a  change  in  rate  of  movement  occurs  only  at  the 
initial  contact  between  forces  and  at  breakpoints.  The  rate 
of  advance  of  the  BED  forces  should  be  allowed  to  vary  based 
on  decision  logic  which  considers  attrition  rate,  terrain, 
etc.  BLOB  forces  should  be  given  the  capability  to  move  to 
alternate  or  fall  back  positions  when  the  battle  situation 
dictates.  Decision  logic  must  be  enhanced  to  allow  BLUE 
forces  to  move  to  alternate  or  fail-back  positions.  It 
should  be  noted  that  this  "pace  of  battle"  question  has 
never  been  fully  explained  or  understood,  however,  the  move¬ 
ment  routine  logic  in  DAHSTAC  should  parallel  that  of  STAR. 

The  concept  of  units  having  different  formations  is 
ignored  in  the  current  model.  Units  are  identified  by  a 


point  location  only.  Some  work  needs  to  be  done  to  allow 
the  units  to  take  on  different  foraations  for  road  aarches, 
attacks,  etc.  In  order  to  do  this,  additional  attributes 
such  as  frontage,  depth,  orientation  and  geoaetric  configu¬ 
ration  should  be  considered  to  describe  the  formation  of  the 
unit.  &  postprocessor  for  STAR  could  take  the  individual 
tank  locations  froa  the  STAR  data  base  and  deteraine  the 
above  mentioned  attributes  for  input  into  the  DABSTAC  aodel. 
Attrition  rates  on  both  sides  must  be  adjusted  for  the  vari¬ 
ous  formations  of  an  attacking  unit. 

The  inclusion  of  indirect  fire,  suppression,  smoke, 
intelligence,  and  close  air  support  are  all  important  to 
making  the  aodel  a  usable  analytical  tool.  Additions  of 
communications  and  logistics,  primarily  ammunition  and 
petroleum,  oil,  and  lubricants  (POL) ,  are  also  important 
enhancements  that  must  be  included  in  the  aodel. 

As  was  seen  in  section  III-E-3,  Pva,the  probability  that 
a  target  is  visit  a  and  acguired  in  the  steady  state  is 

*  (  x  *  „  )  (j  «15* 

Both  of  the  factors  represent  a  particular  element  of  the 
attrition  process,  one  the  terrain  affects  on  line  of  sight. 
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and  the  other,  the  acquisition  process.  It  should  be  appa¬ 
rent  that  these  two  factors  can  be  separated  froa  each  other 
with  the  only  coaeon  variable,  y  r  the  rate  of  transfer  froa 
the  visible  state  to  the  invisible  state.  Kith  this  factor- 
ization,  if  the  line  of  sight  parameters  for  a  particular 
terrain  type  are  known  and  the  attrition  rate  coefficients 
have  been  determined,  the  target  acquisition  capability  of  a 
weapon  system  against  another  weapon  system  can  be  deter¬ 
mined.  This  would  allow  the  target  acquisition  parameter  to 
be  used  for  other  terrain  types  to  compute  attrition  rate 
coefficients  without  having  to  estimate  the  coefficients 
froa  a  high  resolution  simulation. 

The  concept  works  in  this  manner.  Buns  are  made  on  the 
high  resolution  model  (STAB)  using  a  particular  terrain. 

The  attrition  coefficients  are  determined  using  the  maximum 
likelihood  estimates  described  in  Section  III-D.  With  pre¬ 
determined  data  for  the  terrain  coefficients,  the  line  of 
sight  parameter  can  be  factored  out  of  the  attrition  rate 
coefficient,  A.  This  would  leave  the  value  of  the  acquisi¬ 
tion  factor  which  should  theoretically  be  constant  for  the 
firer-target  type  combination  used  in  the  scenario.  With 
this  value  and  the  terrain  coefficient  for  other  types  of 


terrain,  the  attrition  rate  coefficient  can  be  determined 


for  a  new  scenario  with  the  sane  firer-target  combination, 
without  running  the  high  resolution  aodel  again.  This 
method  assuaes  that  the  terrain  affects  and  the  target 
acquisition  process  are  independent  of  each  other.  The 
iapact  of  this  concept  is  enoraous  in  terns  of  tine  and  cost 
saved  froa  not  having  to  make  multiple  runs  of  the  high 
resolution  model. 

One  concept  which  is  not  strictly  an  enhancement  but  is 
an  important  consideration  is  the  idea  of  using  simulated 
coabat  data  froa  a  source  other  than  STAB,  specifically, 
data  generated  at  the  National  Training  Center  from  actual 
field  exercises  could  conceivably  be  used  to  compute  attri¬ 
tion  coefficients  in  auch  the  same  manner  as  STAB.  What 
this  provides  is  not  only  an  alternate  source  of  data,  but 
also  a  different  aodel,  since  it  would  include  many  human 
factors  that  would  not  likely  be  accounted  for  in  STAR. 
Depending  on  the  intended  use  of  the  model,  one  source  or 
the  other  may  be  desirable. 

It  is  necessary  that  a  progression  froa  homogeneous 
forces  to  heterogeneous  forces  be  made  if  DAMSTAC  is  ever  to 
be  really  useful.  This  enhancement  entails  a  significant 
amount  of  research,  aodel  building,  and  statistical  analysis 
in  order  to  iapleaent  the  concept  of  heterogeneous  forces 
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into  the  analytical  aodel.  Target  priorities,  fire 
allocation  algorithms,  and  a  library  of  attrition  rate  coef¬ 
ficients  sust  be  developed  for  various  scenarios  and  force 
sixes.  Changing  the  aodel  to  a  heterogeneous  type  requires 
soae  changes  in  the  equations  (13)  and  (14),  and  thus  soae 
significant  modifications  to  the  postprocessor  for  the 
aodel. 
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71.  HBTEROGEHEOUS  J1QDE£ 

DAHSTAC,  in  its  current  configuration,  aodels  combat 
between  two  opposing  homogeneous  forces,  but  actual  coabat 
consists  of  many  different  wea pon- system  types  operating 
together  as  coabined  aras  teams.  For  ezaaple,  there  may  be 
infantry,  tanks,  antitank  weapons,  artillery,  etc.  on  each 
side.  A  necessary  extension  or  Modification,  therefore,  to 
DAHSTAC  is  the  consideration  of  coabat  between  heterogeneous 
forces.  This  chapter  briefly  indicates  how  the  basic  ideas 
and  methodologies  used  previously  to  model  coabat  attrition 
for  hoaogeneous  forces  can  be  extended  and  adapted  to  the 
heterogeneous  force  structure. 

A  key  assuaption  that  aust  be  made  in  order  to  model 
coabat  between  heterogeneous  forces,  where  tanchester-type 
differential  equations  are  used,  is  that  the  attrition 
affects  of  various  different  enemy  weapon  systems  types 
against  a  particular  friendly  target  type  are  additive.  In 
other  words,  synergistic  effects  between  weapon  system  types 
are  not  considered  (Ref.  9].  Also  inherent  in  the  modelling 
of  heterogeneous  force  coabat  is  the  need  to  introduce  an 
allocation  factor.  This  allocation  factor  may  be  defined  as 
the  fraction  of  the  firers  of  a  particular  type  that  engage 
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a  particular  target  type.  The  allocation  factor  concept 
will  be  discussed  in  greater  detail  later  in  this  chapter. 

One  theoretical  concept  for  aodelling  the  attrition  pro* 
cess  between  two  opposing  forces  is  to  consider  attrition  as 
a  nonstationary  Markov  process  where  the  states  of  the  pro¬ 
cess  are  defined  as  the  nuabers  and  types  of  surviving  com¬ 
batants  and  the  ranges  between  opposing  coabatants.  The 
transition  probabilities  between  the  Markov  states,  under 
this  concept,  depend  on  the  ratio  at  which  surviving  coabat¬ 
ants  or  weapon  systeas  inflict  casualties  on  a  particular 
type  of  opposing  combatant.  Thus  the  nonstationary  Markov 
process  model  provides  a  means  of  relating  the  sequence  of 
states  and  the  times  of  the  state  transitions  observed  dur¬ 
ing  a  battle  to  attrition  rate  coefficients.  In  turn,  this 
allows  for  the  estimation  of  attrition  rate  coefficients  on 
the  basis  of  the  observed  sequence  of  states  and  transition 
tines.  Indrighetti  [Ref.  10]  developed  such  a  nonstationary 
Markov  model  that  related  weapon  system  effectiveness  to  the 
tine  seguence  of  casualties  observed  in  a  two-sided,  hetero¬ 
geneous  force  land  combat  simulation.  There  is  a  flaw,  how¬ 
ever,  in  kndrighetti' s  work  that  readers  should  note.  He 
implies  that  his  methodology  is  easily  extended  to  time-de¬ 
pendent  attrition  rate  coefficients.  This  is  not  true. 
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Andrighetti 's  thesis  does  contain  a  good  discussion  on 
the  theoretical  background  which  allows  for  the  maximum 
likelihood  estimation  of  paraaeters  used  in  attrition  rate 
equations.  An  iaportant  observation  made  by  Andrighetti  in 
his  thesis  is  that  the  difference-differential  equations 
obtained  froa  the  Chapman- Kolmogorov  equations  of  the 
nonstationary  Harkov  model  are  susceptible  to  a  general 
solution  only  in  special  cases  as  Gordon  Clark  showed  in  his 
doctoral  dissertation.  Clark  used  the  Chapaan-Kolnogorov 
equations  to  develop  expressions  for  time  state  probabili¬ 
ties  and  the  expected  survivors  in  a  heterogeneous  force 
battle. 

The  extended  version  of  the  DAHSTAC  model  which  shall  be 
called  DAHSTAC  II  should  be  a  heterogeneous  force  analytical 
model  that  uses  STAH  output  to  compute  parameter  values. 
These  parameter  values  should  be  allowed  to  change  as  the 
distance  or  range  between  opposing  forces  vary  during  the 
battle.  In  DAHSTAC,  an  attempt  was  made  to  allow  the  param¬ 
eter  values  to  change  with  respect  to  time.  However,  it  was 
discovered  that  the  maximum  liklihood  estimators  developed 
for  the  critical  parameters  in  the  attrition  rate  equations 
could  not  be  used  for  time  dependent  variables.  Two  feasi¬ 
ble  alternatives  that  allowed  for  variable  parameter  values 
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were  considered.  Both  alternatives  are  recommended  for 
future  evaluation.  The  first  alternative  looks  at  dividing 
the  battle  into  phases  based  on  battle  intensity  as  measured 
by  the  number  of  shots  divided  by  the  nuaber  of  live  firers. 
(Refer  to  Section  III-E)  The  second  alternative  for  divid¬ 
ing  the  battle  into  phases  is  based  on  the  range  between 
coabatants.  Here  also  the  paraaeters  can  change  froa  phase 
to  phase.  This  technique  will  work  assuaing  that  the  com- 
batants  on  each  side  nove  as  a  single  force  and  that  units 
are  located  by  their  centers  of  aass.  For  exaaple,  a  battle 
can  be  divided  into  three  distinct  tine  intervals  or  phases 
of  conbat:  an  artillery  prep  phase  during  which  all  conbat 
vehicles  and  weapon  systems,  other  than  artillery,  are  not 
participating;  a  moderate  to  long  range  battle  phase;  and,  a 
terminal  shorter  range  battle  phase.  Parameter  values  can 
then  be  computed,  using  procedures  to  be  discussed  later, 
for  each  battle  phase.  This  second  alternative  is  used  in 
the  COh.’-.NE¥  model  and  is  the  one  recommended  for  DAMSTAC  II, 

A.  THE  ATTRITION-RATE  FUNCTIONAL  FORM 

The  attrition  rate  functional  form  that  will  ultimately 
be  chosen  for  the  heterogeneous  force  model  is,  of  course, 
the  prerogative  of  the  modeller.  The  authors  of  this  the¬ 
sis,  however,  would  like  to  recommend  two  related  attrition 
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rate  functional  foras  for  consideration.  The  first,  which 
shall  be  labeled  as  Functional  Fora.  1  (FF1) ,  is  based  on  the 
assumption  that  during  the  course  of  a  battle  fire  is 
directed  at  the  highest  priority  target  surviving  and 
acquired.  Also  it  is  assuaed  that  the  relative  priority  of 
target  weapon  types,  as  considered  by  each  firer,  remains 
constant  throughout  the  engagement  and  that  this  relative 
priority  is  identical  for  each  firing  weapon  type.  Func¬ 
tional  Fora  1,  developed  by  Gordon  Clark  for  use  in  the 
COHAN  model,  can  be  written  as  follows: 


f“  (inn) 


E 

U*i 


<*ijk  n.  P5m*(l  -  P*"1) 


mj 


(16) 


where 

(m,n)  =  conditional  casualty  rate  of  BLOE  weapons  of 
type  j  during  phase  i  of  the  battle  given  strengths  of  m 
and  n. 

til 

f^,  («#n)  =  analog  of  f^  for  RED  weapons  of  type  k. 
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nw  =  the  number  of  surviving  BED  weapons  of  type  It 

»  the  number  of  surviving  BLUE  weapons  of  type  j 

J 

Jcill  rate  of  a  BL0E  firer  of  type  j  for  an  acquired 
BED  target  of  type  k  in  the  i'th  tine  interval. 

I 

kill  rate  for  a  red-f irer/blue-target  combination 
corresponding  to 


q;  =  the  probability  that  a  specific  BED  weapon  is  unac¬ 
quired  by  an  individual  BLUE  firer  in  the  i'th  phase  of 
the  battle. 


p;  a  the  probability  for  a  specified  bl ue-target/red-firer 
combination  corresponding  to  g; . 


mj  -  the  number  of  surviving  BLOE  weapons  of  higher  prior¬ 
ity  than  a  type  j  weapon. 
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n'k  =  the  number  of  surviving  BED  weapons  of  higher  prior¬ 
ity  than  a  type  k  weapon. 

The  reader  should  note  that  since  the  assumption  is  Bade 
that  a  firer  will  engage  the  highest  priority  target  he  can 
acquire,  the  probability  of  a  BLUE  firer  engaging  a  RED 

f 

weapon  of  type  k  is  q^ll -qjS  .  Similarly,  the  probability  of 
a  RED  firer  engaging  a  BLUE  weapon  of  type  j  is  ffl(1-f\"$  . 
These  two  probabilities  may  be  regarded  as  the  allocation 
factors  ft  r  Functional  Fora  1. 

All  of  the  critical  parameters  necessary  to  compute  the 
conditional  casualty  rates,  f^  (m,n)  and  f^  (m,n) ,  can  be 
obtained  either  directly  or  indirectly  (using  maximum  like¬ 
lihood  procedures)  from  the  output  of  a  high  resolution 
simulation  such  as  STAB. 

If  one  is  not  willing  to  assume  that  target  priorities 
remain  constant  throughout  the  battle  then  FF1  will  not  be 
sufficient.  The  terms  in  Clark's  functional  fora  which 
result  from  the  assumption  of  constant  target  priorities  are 
qj'*  and  p-^i  where  q|^  is  the  probability  that  all  surviving  RED 
weapons  of  higher  priority  than  type  k  have  not  been 
acquired  and  p^"*1  is  the  corresponding  probability  for  type  j 
blue  weapons. 
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One  way  to  extend  ClarJc's  formulation  in  order  that  the 
fixed  priority  may  be  dynamically  overridden  is  by  replacing 
the  terms  qN  with©jk  where  9}k  is  defined  as  the  probability 
that  BIOS  weapon  type  j  engages  BED  weapon  type  k,  given 
that  BLUE  weapon  type  j  acquires  at  least  one  BED  weapon 
type  k.  Note  that 

^k(1“<3,l,‘)  =  P(j  engages  k|at  least  one  k  acquisition) 

P  (acquire  at  least  one  k)  =  P(j  engages  k) 


Functional  Form  2  (FF2)  may  be  written  as; 


f“  (m,n)  =  £  nk  (1-p*4)  (17) 

f“k  (m,n)  =  ®j  9;k  (1-qJS  (18) 

where  5k^  is  the  RED  firer  BLUE  target  probability  corre¬ 
sponding  to  9jk  and  all  other  parameters  are  as  defined  for 
FF1.  Just  as  p,q,a,£  are  assumed  constant  over  the  entire 
phase  of  the  battle,  <5  and  0  ,  for  each  j  and  k,  are  also 
assumed  constant  over  the  entire  interval. 
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Serious  consideration  should  be  given  to  incorporating 
PP2  or  soae  similar  functional  form  into  DAHSTAC  ZZ  because 
STAB  does  allow  target  priorities  to  change  during  a  battle. 
Target  priorities  in  STAB  are  based  on  the  types  of  targets 
on  the  firer's  detected  list  and  the  type  of  ammunition  the 
firer  has  available  for  use. 

Zt  may  also  be  useful  to  use  a  different  functional  form 
during  the  initial  phase  of  combat  where  attrition  rates  are 
more  a  function  of  the  exposure  or  vulnerability  of  targets 
than  they  are  a  function  of  the  firer's  capability  to 
acquire  targets.  Clark  recommended  the  use  of  Peterson's 
Logarithmic  Law  to  model  attrition  during  the  initial  phase 
of  the  battle.  [Ref.  11] 

B.  MAXZHOH  LIKELIHOOD  ESTZHATOBS 

The  ability  of  the  DAHSTAC  IZ  model  to  provide  insight 
into  the  interactions  represented  in  the  STAB  simulation  is 
based  on  the  estimation  of  conditional  casualty  rates  (  a 
and  0  )  and  other  parameters  from  STAR  data.  These  estimated 
parameters  should  be  constant  within  a  phase  of  the  battle, 
but  independent  and  unrelated  to  results  in  neighboring  time 
intervals.  The  independence  of  parameter  values  between 
time  intervals  implies  that  the  estimators  for  a  specific 
time  interval,  or  phase  of  the  battle,  can  be  defined  by 
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analyzing  a  sample  of  data  froa  the  high  resolution  simula¬ 
tion  observed  during  that  single  -time  interval. 

The  aaxiaua  likelihood  estiaators  (ALE'S)  for  the 
unknown  parameters  a ,  0,  p  and  g  used  in  FP1  were  developed 
by  Gordon  Clark.  It  aay  be  of  interest  to  the  reader  to 
note  that  simplified  versions  of  these  ALE's  were  used  in 
the  DAMSTAC  modal  for  homogeneous  forces,  however,  since  we 
are  now  concerned  with  heterogeneous  forces,  the  "una¬ 
bridged"  versions  are  presented.  The  aaziaua  likelihood 

estimate  for  a  as  a  function  of  p  is: 

* 

V0) 


3j4 


(1  -  pm,i  )(ty  -  tY.») 


for  >  0 


(19) 


^  0  otherwise 

where  the  estimator  for  the  parameter  a  for  a  given  p. 

The  ALE  for  p  is  found  by  solving: 

3ln  L(o)  Y  (1  -  CyXnu^  -(m^  + 

a?  “  0(1-?^) 

(2°) 

-Y^T. (  at- ( o )n^  o^ti  (tY  -  tY-4  )  . 

Y««  j.«  t*i 

(m^|  -(mYj  )r5m,,i  ))  =  0 

The  equations  for  (g)  and  q  are  similar. 

The  estiaators  obtained  froa  the  above  equations  result 
froa  the  data  froa  a  single  replication  of  the  coabat  simu¬ 
lation.  lore  precise  estimates  are  obtainable  by  using  a 
larger  3aaple  consisting  of  a  number  of  observations  during 
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the  interval  in  question  from  a  number  of  independent 
replications  of  the  STAB  siaulation.  This  aodification 


entails  an  additional  suaaation  (£  *  where  x  =  the  total 

x 

number  of  replications)  in  each  of  the  estimating  equations. 
The  data  needed  for  the  computations  of  the  HLE's  include; 
y  *  total  number  of  casualties 

ty  »  elapsed  time  since  the  beginning  of  the  interval 
until  the  y'th  casualty 

tyt,=  the  end  of  the  battle  phase  or  the  elapsed  time  since 
the  beginning  of  the  phase  until  the  battle  ended  if  it  is 
less. 

m^i  =  the  number  of  BLUE  survivors  of  type  j  just  prior  to 
the  y *  th  casualty 

n^  *  the  number  of  BED  survivors  of  type  k  in  the  battle 
just  prior  to  the  y'th  casualty 

=  1  if  the  y'th  casualty  is  BLUE  ;  0  otherwise 

f Y  *  firer  weapon  type  for  the  y'th  casualty 
3  target  weapon  type  for  the  y'th  casualty 
my  ®  the  number  of  BLUE  survivors  of  higher  priority  than 
type  j  just  prior  to  the  y'th  casualty 

■ii  *  E 

<**• 

n^  *  BED  analog  of  ntyi 
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Bki  *  the  number  of  BIOS  weapon-type  j  casualties  due  to 
RED  weapons  of  type  k 

One  technique  of  coaputinq  the  HLE's  for  a  and  p  is  to 
nueerically  deteraine  p  using  equation  (20)  and  then  substi¬ 
tute  this  value  into  equation  (19)  for  j  =  1,2,...,b  and  for 
k  =  1,2,...,r  (b  =  the  total  nuaber  of  BLUE  weapon  types  and 
r  =  the  total  nuaber  of  RED  weapon  types).  A  recoaaended 
technique  for  numerically  solving  for  p,  the  estiaate  for  p, 
can  be  found  in  Clark's  work  on  the  Dyntacs  combat  model 
[Ref.  12].  This  technique  is  similar  to  the  one  developed 
for  the  homogeneous  case  discussed  in  Chapter  III.  Basi¬ 
cally  the  procedure  is  as  follows.  Since  p  represents  a 
probability,  it  must  satisfy  the  inequality  G<p<1.  There¬ 
fore,  the  function  Ln  l(p)  can  be  evaluated  at  points  in 
the  interval  0Sp<1  until  two  points  ,  p,  and  p*  ,  are  found 
such  that  the  product  of  ^  Ln  L(pt)  and  Ln  L(pz)  is  less 
than  or  equal  to  zero.  Since  Ln  L(p)  is  a  continuous 
function.  ^  l.  MB  =0  for  so.,  p  satisfying  p,  s  P  S  P,  . 
This  interval  is  reduced  by  evaluating  ■ip”  Ln  L (p)  at  the 
midpoint  of  the  interval  (P|»P*)*  Repeating  this  process 
converges  to  a  solution  for  Ln  L(p)  and  this  solution  is 
a  local  maximum.  This  technique  is  similar  to  the  bisection 
method  for  solving  nonlinear  problems. 
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lith  respect  to  Functional  Fora  2 ,  the  introduction  of 
the  allocation  teras  adds  considerable  coaplexity  to  the 
maximization  of  the  likelihood  function.  One  way  to  circus- 
went  this  probles  is  to  coapute  a  "best  guess"  for  the  9^  ' s 
froa  STAB  saaples  (for  the  given  interval) .  This  best  guess 
can  thenbe  substituted  into  the  ©jk  's  prior  to  aaxiaization 
of  the  likelihood  function,  thereby  reducing  thee^'s  to 

constants.  The  aaxiaization  can  then  be  perforaed,  as  in 

» 

Clark’s  original  eguations,  with  respect  to  p,  q,a  ,  0  . 

C.  THE  MODEL 

The  coaputational  procedures  that  will  constitute 
DANSTAC  II  aust  be  developed  by  the  aodeller  who  chooses  to 
carry  on  the  work  initiated  in  this  thesis.  As  a  starting 
point,  however,  the  techniques  used  by  Clark  in  his  COMAN 
aodel  would  be  excellent  choices.  The  adequacy  of  Clark's 
procedures  in  DAHSTAC  II  given  a  particular  scenario  and 
force  six,  to  predict  attrition  in  a  STAB  siaulation  with 
the  saae  scenario  and  force  mix  is  questionable  and  would 
have  to  be  tested.  Detailed  conparisons  and  analysis  of  the 
output  froa  both  aodels  will  be  required  to  answer  the  ques¬ 
tion  of  aodel  coaparability  and  modifications  to  the  com¬ 
puter  code  and/or  coaputational  procedures  in  DAHSTAC  II 
aust  be  accoaplished  as  necessary.  The  coaputational  cycle 
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employed  in  COBAN  is  based  on  the  assumption  that  the  dis¬ 
tribution  tines  between  casualties  in  i  a  high  resolution 
simulation  is  exponentially  distributed..  It  uses  Konte 
Carlo  procedures  to  generate  the  time  until  the  next  casu¬ 
alty  and  also  to  find  the  casualty  weapon  type. 

The  following  steps  constitute  the  computational  cycle 
recommended  for  DABSTAC  II. 

Step  1.  The  first  step  in  the  cycle  sets  the  time  equal 
to  zero  and  sets  the  battle  phase,  designated  by  the  letter 
i,  equal  to  1. 

Step  2.  Next  the  total  conditional  casualty  rate, 
is  computed  from  the  equation 

for  1=1,2,—,  I  (21) 

where  fj^  (m,n)  and  (m,n)  are  the  conditional  casualty 
rates  in  phase  i  of  the  battle,  given  force  strength  i  and 
n,  for  BLOB  weapon  type  j  and  BED  weapon  type  k,  respec¬ 
tively.  A*  (m,n)  is  the  parameter  of  the  exponential  distri¬ 
bution  describing  the  tines  to  the  next  casualty. 

Step  3.  Now,  using  the  exponential  distribution  with 
parameter  A'(m,n),  one  Monte  Carlo's  for  the  time  interval 
to  the  next  casualty.  Designate  this  tine  interval  as  t'. 

If  the  next  casualty  occurs  in  another  tine  interval,  ie.  if 
t  ♦  t*  >  t;  where  t;  represents  the  end  of  the  i*th  battle 
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phase,  then  i  is  increased  by  1,  the  tine  is  set  equal  to  t 

{ 

>  and  the  program  is  directed  back  to  step  2  where  V*‘(m,n)  is 

selected  as  the  distribution  parameter.  If,  however,  t  ♦  t* 
<  t;  then  the  next  step  is  to  Monte  Carlo  for  the  casualty 
weapon  type. 

Step  4.  The  weapon  type  for  the  next  casualty  is  repre¬ 
sented  by  the  random  variables  and  Cr  which  are  defined 
as 

*  j  if  the  casualty  weapon  is  BLOB  weapon  type  j 
0  if  the  casualty  is  BED 

Cr  *  k  if  the  casualty  weapon  is  BED  weapon  type  k 
0  if  the  casualty  is  BLOE 

The  Monte  Carlo  procedure  uses  the  conditional  casualty 
weapon  probability  distributions  for  Cb  and  Ct  .  These 
probabilities  are  computed  by 

p*  <Cb  *  ;J(C  *  1,a,n)  *  fj!  (a,o)  /\(myn) 

P'  (Cr  *  k |  C  =  1  ,m,n)  ■  {n,n)/  ^m(n) 

for  i=  1,2,...,  I  and  when  variable  C  is  set  equal  to  1  if 
a  casualty  is  determined  to  occur  in  interval  i  from  the 
previous  steps.  The  probability  that  the  casualty  is  of  a 
particular  type  is  computed  for  each  BED  and  BLUE  type 
weapon.  The  sum  of  these  probabilities  should  equal  1. 
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Step  5.  The  next  step  is  to  Bonte  Carlo  using  the  Oni- 
f ora  (0,1)  distribution  to  determine  the  casualty  type.  An 
exaaple  aay  help  in  understanding  the  above  procedure.  Sup¬ 
pose  there  are  two  types  of  BLUE  weapons  and  two  types  of 
SED  weapons.  Assuae  the  conditional  casualty  probabilities 
for  each  weapon  are  computed  using  the  equations  shown  pre¬ 
viously  and  they  are  0.25  for  each  weapon  type,  bl,  b2,  rl, 
r2.  If  one  orders  the  weapon  types  in  soae  arbitrary  fash¬ 
ion  such  as  bl,  rl,  b2,  r2,  and  then  accuaulates  the  prob¬ 
abilities  on  the  interval  (0,1),  one  can  assign  an  interval 
to  each  weapon  type.  The  exaaple  would  yield  the  following 
int  arvals: 

bl  -  (0,0.25)  ; 
rl  -  (0.26,0.50)  ; 
b2  -  (0.51,0.75)  ; 
r2  -  (0.76,1.00)  . 

The  uniform  random  number  can  then  specify  the  casualty 
weapon  type.  If  the  random  number  is  0.60  then  the  casualty 
weapon  type  in  our  example  is  b2. 

Step  6.  After  the  casualty  weapon  type  has  been  found 
the  next  step  is  to  decrement  the  number  of  that  type  weapon 
by  1.  Now  increase  the  value  of  t  by  t* ,  the  time  to  the 
next  casualty.  If  t  is  greater  than  or  equal  to  which 
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represents  the  ending  tiae  of  the  last  phase  cf  the  battle, 
or  if  either  side  is  annihilated  ,  the  battle  has  ended  and 
the  program  terminates.  If  t  is  less  than  tx  then  go  to 
step  2. 

Observe  that  the  above  computational  procedures  are  sto¬ 
chastic  in  nature  and  different  replications  of  an  engage¬ 
ment  will  likely  yield  similar  but  different  results.  If 
thi3  method  is  chosen  for  DAHSTAC  II,  variance  reducing 
replications  for  each  engagement  should  be  performed. 

D.  THE  STAR  POSTPROCESSOR 

It  is  very  important  that  anyone  who  attempts  to  develop 
the  DAHSTAC  II  model  understand  and  be  able  to  manipulate 
the  STAR  simulation  model.  In  the  initial  attempts  to 
extract  data  from  STAR,  certain  subroutines  must  be  turned 
off.  These  routines  include  the  Smoke,  Engineer,  close-air 
support,  etc.  whose  effects  on  parameter  estimates  used  in 
DAHSTAC  II  have  not  been  thoroughly  analyzed.  The  number  of 
type?  of  combatants  played  in  the  simulation  will  have  to  be 
regulated  so  that  they  do  not  overpower  the  capabilities  of 
the  analytical  model.  It  is  recommended  that  no  more  than 
three  different  types  of  weapon  systems  for  each  opposing 
force  be  attempted  when  initially  formulating  the  heteroge¬ 
neous  force  model. 
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It  is  iaportact  that  the  finished  version  of  the 
analytical  aodel  be  as  sinple  and  transparent  as  possible 
and  not  encumbered  with  parameter  estimating  algorithms  and 
computations.  This  consideration  plus  the  amount  of  data 
that  must  be  extracted  from  the  STAB  output  require  the 
development  of  a  STAB  postprocessor  program.  As  an  illus¬ 
tration,  if  one  considers  a  simple  battle  divided  into  three 
battle  phases  with  only  three  types  of  weapon  systems  on 
each  side,  the  DAHSTAC  II  model,  as  it  has  been  described, 
requires  the  estimation  of  114  parameter  values  and  the  com¬ 
putation  of  six  attrition  rates. 

The  postprocessor  at  Appendix  E,  developed  for  the  homo¬ 
geneous  model,  should  provide  a  good  foundation  for  the  more 
complicated  postprocessor  needed  for  DAHSTAC  II. 
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?II .  CO NC LODI  MG  BEHABKS 

The  concept  of  using  differential  equations  to  represent 
coabat  began  with  Lanchester's  classical  equations  and  has 
developed  into  a  very  dynaaic  field  of  study.  The  theoreti¬ 
cal  background  of  the  fitted  parameter  method  has  been  well 

docuaented  by  Clark  and  Taylor.  However#  the  practical 
application  of  the  theory  lacks  this  same  degree  of  discus¬ 
sion  and  documentation.  This  paper  has  focussed  on  some  of 
the  techniques  of  applying  Clark's  theory  to  an  actual  high 
resolution,  Honte  Carlo  simulation.  The  problems  are  real 
and,  as  has  been  seen#  the  solutions  presented  are  not 
doctrine.  It  is  hoped  that  further  investigation#  elabora¬ 
tion  and  documentation  of  this  modelling  technique  will 
occur  zn  the  near  future. 

The  choice  of  which  functional  form  to  use  in  a  particu¬ 
lar  analytical  coabat  model  is  entirely  up  to  the  discretion 

of  the  aodel  builder.  In  choosing  an  attrition  rate  func¬ 
tional  form  the  model  builder  should  consider  the  scenario 
portrayed  in  the  model  and  the  assumptions  that  are  compati¬ 
ble  with  it. 

Finally,  it  is  recommended  that  the  entire  program  for 
the  aodel  be  converted  to  the  0C5D  Pascal  language  on  the 
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Apple  II.  The  basic  structure  of  Pascal  is  more  conducive 
to  coabat  modelling  siaulation  than  the  structure  of  the 
BASIC  language.  The  concept  of  entities  and  attributes 
found  in  SISSCHIPT  II. 5  can  be  closely  emulated  in  Pascal 
while  this  is  not  as  easily  accomplished  or  clearly  defined 
in  BASIC. 

The  aodular  structure  of  procedures  and  functions  in 
Pascal  is  also  an  advantage  in  the  enhanceaent  process  of 
aodel  building.  Modellers  can  use  one  procedure  to  repre¬ 
sent  a  particular  action  or  phenomenon  occurring  on  the  bat¬ 
tlefield.  Any  action  which  is  unable  to  be  modelled  can  be 
left  in  the  program  as  a  dummy  procedure  for  later  work. 

Additionally,  the  Pascal  language  allows  more  descrip¬ 
tive  variables  to  be  used  than  the  two-character  variables 
used  in  Applesoft.  This  would  make  the  program  itself  more 
transparent  and  easier  to  modify  and  update. 

The  DAMSTAC  model  is  clearly  in  its  infancy  stage. 
Through  continued  enhancement,  this  model  can  contribute  to 
the  ongoing  work  in  the  field  of  combat  analysis. 
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JcISX  0£  VARIABLES 

A  =  casualty  rate  for  RED  attritting  acquired  BLOB  targets 
AM  =  angle  in  radians  measured  counter-clockwise  from  due 
east  sade  by  the  line  between  the  current  unit 
location  and  the  next  coordinating  point. 

B  *  casualty  rate  for  BLUE  attritting  acquired  RED  targets 
BE  =  the  breakpoint  for  the  BLUE  force 
BH  =  movement  criteria  for  BLUE  units  as  a  fraction  of 
unit  strength 

BF (I)  *  number  of  RED  targets  being  fired  at  by  BLUE  unit  I 
BP (I)  =  aost  recently  passed  coordinating  point  for  BLUE 
unit  I 

BR(I)  =  rate  of  movement  for  BLUE  unit  I  in  fl/SEC 
BT(I,J)  »  attrition  (casualties)  caused  by  RED  unit  J  on 
BLUE  unit  I  in  one  on  one  battle 
BX (I, J)  =  x-coord  of  coordinating  point  J  for  BLUE  unit  I 
BY(I,J)  *  y-coord  of  coordinating  point  J  for  BLUE  unit  I 
BO (I)  =  current  total  BLUE  force  strength  for  time  period  I 
BO (0)  =  initial  BLUE  force  strength 

B1  =  number  of  coordinating  points  for  BLUE  (including 
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initial  location) 


B1 (I)  *  initial  strength  of  BLOB  unit  I 

DB (I)  *  nueber  of  casualties  or  aeount  of  attrition  for 
BLOB  unit  I 

DH  ~  distance  travelled  by  a  unit  in  time  increment  DT,  at 
speed  BT 

DI  =  distance  betveen  current  unit  location  and  nezt  turn 
point 

DR ( J)  =  nueber  of  casualties  or  aeount  of  attrition  for  RED 
unit  J 

DT  -  tine  increeent 

DT  3  x  component  of  distance  between  current  unit  location 
and  next  coordinating  point 

DT  *  y  component  of  distance  betveen  current  unit  location 
and  next  coordinating  point 

DI  *  x  component  of  distance  moved  in  time  increment  DT,  at 
speed  RR  (I)  or  BR(I) 

D2  »  y  component  of  distance  moved  in  time  increment  DT,  at 
speed  RB(I)  or  BR(I) 

H  -  number  of  BLOB  units 

NT  3  max  time  of  simulation 

N  =  number  of  RED  units 

N$  =  file  name  (used  recursively) 
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P  =  probability  that  a  BLOB  target  is  not  acquired  by 
a  BED  unit 
PI  =  3. 141592654 

PB  =  status  of  printer;  0  =  printer  is  off 

1  =  printer  is  on 

Q  =  probability  that  a  BED  target  is  not  acquired  by 
a  BLUE  unit 

RE  »  the  breakpoint  for  the  BED  force 

BF (J)  s  nuaber  of  BLUE  targets  being  fired  at  by  BED  unit  J 
RG  *  range  between  a  BLUE  unit  and  a  BED  unit 
RH  =  movement  criteria  for  RED  units  as  a  fraction  of 
unit  strength 

RP (I)  *  most  recently  passed  coordinating  point  for  RED 
unit  I 

BR(I)  *  rate  of  movement  for  BED  unit  I  in  8/SEC 
RT(J,I)  =  attrition  (casualties)  caused  by  BLUE  unit  I  on 
RED  unit  J  in  one  on  one  battle 
BX  (I ,  J)  =  x-coord  of  coordinating  point  J  for  BED  unit  I 

BY(I,J)  =  y-coord  of  coordinating  point  J  for  RED  unit  I 

RO (I)  =  current  total  BED  force  strength  for  time  period  I 

RO  (0)  a  initial  total  BED  force  strength 
R1  *  nuaber  of  coordinating  points  for  RED  (including 
initial  location) 
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HI  (I)  s  initial  strength  of  BED  unit  I 
SB  (I)  =  size  of  BLOB  unit  I 


'i 


SR  (I)  3  size  of  RED  unit  I 
T  *  simulation  tiae 
XB(I)  *  current  x-coord  for 
XR  (I)  »  current  x~coord  for 
IB  (I)  =  current  y-coord  for 
TR (l)  =  current  y-coord  for 


BLOE  unit  1 
RED  unit  1 
BLOE  unit  I 
RED  unit  I 
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1.  Data  which  have  been  hardwired  into  the  eodel  are: 

A  =  2.4 

B  *  0.03 
BE  *  0.3 
BT  *  0 
DT  *  30 
BT  =  1500 
P  =  0.99 

PT  ~  3. 141592654 
Q  »  0.93 
HE  *  0.3 

2.  Once  contact  is  aada  (range  <  2700)  ,  the  rate  of  advance 
decreases  to  2  meters/second. 

3.  The  option  exists  within  the  prograa  to  allow  the  aove- 
■ent  of  units  based  on  a  criteria  of  fractional  force 
strength.  Bhen  unit  strength  goes  to  50%  (eoveaent  crite¬ 
ria)  ,  the  unit  will  aove  to  an  alternate  location  (for  BLUE 
units)  or  retreat  along  its  path  of  advance  (for  BED  units) . 
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This  appendix  provides  a  brief  description  of  the  util¬ 
ity  prograas  associated  vith  the  Aggregated  Model  of  STAB  on 
the  Apple  II  Computer.  These  prograas  provide  the  capabil¬ 
ity  to  the  user  to  transfer  the  necessary  data  from  the  STAB 
model  to  this  model  with  the  minimum  effort. 

The  utility  programs  available  are: 

Force  Maker 
Results  Beader 

Both  prograas  are  interactive  with  the  user,  who  answers 
questions  provided  by  the  computer.  These  prograas  use  the 
commands  for  the  Apple  II  Disk  Operating  System  (DOS)  which 
allow  them  to  read/write  Sequential  Text  Files  and  Random 
Access  Files  onto  the  disk.  Familiarity  with  these  types  of 
files  and  the  DOS  would  improve  the  understanding  of  these 
utility  programs.  [Ref.  13] 

In  the  following  program  descriptions,  the  variable  name 
used  in  the  program  is  listed  in  parentheses  after  its 
descript  ion. 
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A.  FORCE  BAKER  PROGRAM 

This  prograa  will  write  data  files  containing  informa¬ 
tion  about  the  two  opposing  forces. 

Force  Baker  produces  a  Sequential  Text  File  containing: 
Force  size  (number  of  units)  (B) 

The  strength  of  each  unit  (SB  (I) ) 

The  number  of  coordinating  points  for  movement  of  each 
unit  (the  first  coordinating  point  is  always  the  initial 
location  of  the  unit)  (B1(I)) 

The  grid  coordinates  of  each  coordinating  point  for  each 
unit  (BZ  (I*  J)  *  B7  (I ,  J) ) 

The  user  has  the  option  to  make  a  file  for  the  BLUE 
force  or  the  RED  force.  After  all  the  data  has  been 
entered,  a  data  check  routine  will  print  the  information  for 
all  units,  so  that  a  check  for  errors  can  be  made.  After 
the  user  has  checked  the  data  and  identified  the  name  for 
the  file,  the  prograa  will  record  the  file  on  the  diskette. 

A  copy  of  the  program  is  included  as  Listing  B  in  Appen¬ 
dix  F  and  a  sample  session  using  this  prograa  is  attached  as 
Figure  12. 
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1RUN  FORCE  MAKER 
MENU: 

1  MAKE  FILE  FOR  SLUE  FCRCES 
:  MAKE  FILE  FOR  -EE  PQRCSS 
0  -2LIT 

WHICH?  1 

*HAT  IS  THE  FORCE  SIZE  OF  SLUE  3 
FOR  UNIT  1 

HOW  MANY  COORDINATING  POINTS  (INCLUDING  START  LOCATION/'’ 

ENTER  GRID  CCCRDS  FOR  COORDINATING  POINT  #  1 
(-COORD  IS  SAlaC 
r -COORD  IS  100*33 

cnc  UNIT  0 

MOW  .MANY  COORDINATING  POINTS  (INCLUDING  START  LOCATION/'1 

INTER  GRID  OOCRDS  FOR  COORDINATING  POINT  4  1 
< -COORD  IS  34758 
r -COORD  IS  PS 100 


“OR  UNIT  3 

MOW  MANY  COORDINATING  POINTS  i INCLUDING  START  LOCATION/" 


INTER  GRID 
\ -COORD  IS 
v -COORD  IS 


CCCRDS  FOR  COORDINATING  POINT  4  1 

37*48 

•a -’503 


***  DATA  CHECK 


UNIT  FT  .(-COORD  ? -COORD 
:  1  341*3  100*33 

I  1  3*758  95100 

3  37*48  °“S97 

AN  CHANGES"’  (  ( N  >  N 

NAME  CF  FILE"  THESIS 

<*c r\»« '  * 

-  MAKS  FOR  SLUE  FORCES 

Z  MAKE  FILE  FOR  RED  FORCES 


WHICH" 


Figure  12: 


Saapla  Session  of 


Force  Saker  Prograa 
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B.  RESULTS  READER  PBOGRAH 

This  program  reads  the  data  file  of  tiae  intervals  and 
corresponding  BLUE  and  RED  force  levels  produced  by  the 
aodel.  Results  Reader  produces  two  sequential  text  files: 
the  first  file  with  suffix  ".COEFS"  contains: 

Conditional  Casualty  rate,  a  (A) 

Conditional  Casualty  rate,  0  (B) 

Probability  that  a  BLUE  target  is  not  acquired  by  a  BED 
firer  (P) 

Probability  that  a  RED  target  is  not  acquired  by  a  BLUE 
firer  (Q) 

The  second  file  with  suffix  ".RESULTS"  contains: 

The  number  of  data  points  for  BLUE  force  level  versus  time 
(H) 

The  number  of  data  points  for  RED  force  level  versus  tiae 
(N) 

Tima  of  the  battle  at  tiae  interval  I  (T) 

BLUE  force  level  corresponding  to  tiae  interval  I  (B0(l)) 
Tiae  of  the  battle  at  tiae  interval  I  (T) 

RED  force  level  corresponding  to  time  interval  I  (R0(I)) 
The  ".RESULTS"  file  contains  soae  duplication  of  data  for 
the  following  reasons.  With  the  ".RESULTS"  file  written  in 
this  format,  and  free  of  additional  data  such  as  that  in  the 
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".COPPS"  file,  it  caa  be  used  directly  by  the  Apple  Plot 
Prograa  to  produce  a  plot  of  the  force  levels  versus  tiae. 

A  listing  of  the  prograa  is  included  as  Listing  C  in  Appen- 
diz  P.  A  saaple  plot  of  the  DASSTAC  output  data  produced  by 
the  Apple  Plot  Prograa  is  included  as  Figure  13.  A  saaple 
of  the  output  of  D ASS I AC  froa  Besults  Header  is  included  as 
Pigura  14. 
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Pigure  13:  Saaple  Plot  of  Besults  Header  Prograa  output. 
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In  addition  to  the  basic  computer  aodel,  work  vas  done 
to  input  the  STAB  terrain  data  base  into  the  Apple  II  compu¬ 
ter.  This  effort  vas  done  during  the  initial  phase  of  the 
thesis  work  based  on  an  erroneous  assumption  that  the 
authors  made  concerning  the  functional  fora  of  the  attrition 
rates.  This  assumption  was  that  the  terrain  of  the  battle¬ 
field  and  line  of  sight  coaputations  must  be  explicitly 
played  in  the  model.  In  the  given  functional  forms  for  the 
attrition  rates,  line  of  sight  and  terrain  are  accounted  for 
by  the  probability  values  of  p  and  g. 

Before  this  assuaption  was  found  to  be  false,  four  util¬ 
ity  programs,  and  a  model  that  computed  the  elevations  of 
the  units  were  developed  and  are  included  in  this  document 
for  possible  future  use.  If  further  work  is  done  with  a 
aodel  on  the  Apple  II  computer,  the  ability  to  input  the 
STAR  terrain  into  the  Apple  II  may  become  an  important  part 
of  the  aodel.  Examples  of  models  that  might  use  the  STAB 
terrain  in  the  Apple  II  are  snail  aggregated  models  with 
functional  forms  that  explicitly  play  line  of  sight,  or 
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snail,  high-resolution  aodels  that  are  suitable  on  the  Apple 
II  coaputer. 

A.  THE  TB BRAID  UTILITY  PROGRAMS 

The  utility  prograas  that  allow  the  STAR  terrain  to  be 
input  in  the  Apple  II  coaputer  are  listed  below. 

flap  Maker 
Hill  Maker 
Hoods  Maker 
Forest  Maker 

These  four  programs  have  two  portions  each:  a  file  maker 
and  a  file  reader.  In  general,  the  file  maker  portion 
allows  the  user  to  input  a  new  file  or  update/change  an  old 
file.  The  file  reader  allows  the  user  to  print  out  on  the 
screen  or  to  a  printer  the  selected  file. 

For  the  Apple  II  combat  model,  a  STAB  terrain  file  (eg. 
FOLTERR  file)  is  divided  into  its  four  main  parts: 

the  list  of  hill  reference  numbers  in  each  grid  square; 
the  hill  parameters; 

the  list  of  forest  reference  numbers  in  each  grid  square; 
the  forest  ellipse  fitting  parameters. 

Each  part  is  created,  respectively,  by  one  of  the  four  util¬ 
ity  prograas:  Map  Maker,  Hill  Maker,  Hoods  Maker,  and  Forest 
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Halter.  When  these  parts  are  being  created  on  the  Apple 
Disk,  the  saae  general  file  naae  aust  be  used  to  identify 
the  data's  terrain  (eg.  FULDA)  .  Each  of  the  utility  pro- 
grans  will  add  its  own  suffix  to  the  general  file  naae.  For 
exaaple,  for  the  Fulda  terrain,  the  four  utility  prograas 
will  produce  files  naaed: 

FULDA. HAP 
FULDA. HILLS 
FULDA. HOODS 
FULDA. FO BESTS 

This  suffix  identifies  the  file  as  a  particular  portion  of 
the  terrain  data,  and  will  allow  the  Apple  coaputer  to 
recognize  it  from  the  other  files  with  the  same  terrain  naae 
when  reading  the  file  for  the  coabat  nodel. 

ME  »ake£  ££ogran 

This  prograa  will  read  and  write  data  files  contain¬ 
ing  a  catalog  by  grid  square  of  the  hills  which  influence 
the  terrain  in  that  grid  square.  This  prograa  produces  an 
equivalent  file  to  that  in  STAB  called  LIST. H.  Hap  Haker 
produces  a  Sequential  Text  File  containing: 

The  size  of  the  battlefield  (noraally  10  by  10)  (L,H) 

The  lower  left  grid  square  in  1000  ueter  units  (LX, LI) 

The  base  elevation  of  the  battlefield  (HL (0,0,0)) 
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For  each  grid  square  I,J: 

(1)  The  nuaber  of  hills  influencing  the  terrain 

(2)  The  reference  nuaber  of  these  hills  (HL(I,J,K)) 

The  user  has  the  option  to  aake  a  new  file  or  access 
an  old  file  for  updating  or  changes.  For  new  files,  the 
program  will  start  with  the  lower  left  grid  square  and  pro¬ 
ceed  froa  west  to  east  then  south  to  north.  The  data  can  be 
sawed  after  the  information  for  a  grid  square  has  been 
entered.  For  old  files,  the  program  will  access  it,  print 
out  the  file  on  the  screen  and  ask  for  the  starting  grid 
square.  The  grid  square  index  is  used  to  identify  in  which 
grid  square  to  start.  The  grid  square  index  is  an  integer 
froa  1  to  L  in  the  east-west  direction,  and  1  to  a  in  the 
south-north  direction  of  an  L  by  H  sized  battlefield.  The 
lower  left  grid  square  is  1,1;  the  lower  right  is  1, 1 ;  the 
upper  left  is  1,W;  and  the  upper  right  is  L,a. 

&  copy  of  the  program  is  attached  as  listing  0  in 
Appendix  F  and  a  sample  file  is  attached  as  Figure  15. 


2.  sill  aaias.  zsaasia 

This  program  reads/writes  data  files  containing  the 
fitting  parameters  for  each  hill  on  the  battlefield.  Hill 
Maker  produces  a  Bandos  Access  Text  Pile  containing: 

The  nuaber  of  hills  on  the  battlefield  (N) 

Por  each  hill  I: 

(1)  X  coord  of  the  center  location  of  the  hill  (H(l,1)) 

(2)  T  coord  of  the  center  location  of  the  hill  (H(I,1)) 

(3)  The  elevation  of  the  hill  top  aeasured  froa 
xero-sea  level  (H(I,3)) 

(4)  The  orientation  angle  of  the  ellipse  aeasured  in 
degrees  count er-clockwise  froa  east  to  the  major  axis. 

(H  (I  ,  4 ) ) 

(5)  The  eccentricity,  defined  as  the  ratio  of  major 
axis  length  to  ainor  axis  length  (H(I,5>) 

(6)  The  spread,  defined  as  the  distance  in  meters  aeas¬ 
ured  along  the  major  axis  from  hill  center  to  the  contour 
line  which  is  50  aeters  down  froa  the  peak  (H(I,6)) 

(7)  The  maximum  height  of  the  normal  curve  describing 
this  hill  mass  <H(I,7)) 

(8)  The  vertical  distance  aeasured  froa  the  peak  beyond 
which  this  hill  is  not  considered  in  computations  (H(I,8)) 
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A  copy  of  the  program  is  included  as  Listing  E  in 
Appendix  ?,  and  a  saeple  file  is  attached  as  Pigure  16. 

3.  Hoods  $ake£  ££0££4£ 

This  program  Mill  read  and  write  data  files  contain¬ 
ing  a  catalog  by  grid  square  of  the  woods  which  influence 
the  terrain  in  that  grid  square.  This  program  produces  an 
equivalent  file  to  that  of  dap  Maker.  Hood s  Maker  produces 
a  Sequential  Text  File  containing: 

The  size  of  the  battlefield  (normally  10  by  10)  (L,H) 

The  lower  left  grid  square  in  1000  meter  units  (LX,LY) 

For  each  grid  square  I,J: 

(1)  The  number  of  forests  influencing  the  terrain 
<FL(I,J,0)) 

(2)  The  reference  number  of  these  forests  (FL(I,J,K)) 

Hoods  Maker  contains  the  same  options  as  those  found 

in  Map  Maker.  A  copy  of  the  program  is  attached  as  Listing 
F  in  Appendix  F ,  and  a  sample  file  is  attached  as  Figure  17. 

4.  Forest  Maker  Program 

This  program  reads/writes  data  files  containing  the 
fitting  parameters  for  each  forest  on  the  battlefield.  For¬ 
est  Maker  produces  a  Random  Access  Text  File  containing: 

The  number  of  forests  on  the  battlefield  (N) 

For  each  forest  I: 
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(1)  l  coord  of  the  center  location  of  the  forest 


<F(I,1>> 

(2)  T  coord  of  the  center  location  of  the  forest 
(Pdri)) 

(3)  The  height  of  the  trees  in  the  forest  (F(I,3)) 

(4)  The  orientation  angle  of  the  ellipse  describing  the 
forest  aeasared  in  degrees  counter-clockwise  froa  east  to 
the  aajor  axis  (F(I,4)) 

(5)  The  length  of  the  seai-aajor  axis  of  the  ellipse 
<F(I,5)) 

(6)  The  length  of  the  seai-ainor  axis  of  the  ellipse 
(F  (1,6)) 

Forest  Baker  contains  the  saae  options  as  those 
found  in  Hill  Baker.  A  copy  of  the  prograa  is  included  as 
Listing  G  in  Appendix  Ff  and  a  saaple  file  is  attached  as 
Figure  18. 

2.  THE  TERRAIN  BODEL  PROGSAB 

A  copy  of  the  aodel  which  explicitly  plays  terrain  is 
included  in  Appendix  F  as  Listing  H.  A  line  drawn  along  the 
edge  of  the  prograa  delineates  the  terrain  aodel  portion 
which  was  added  to  the  basic  aodel.  The  additional  varia¬ 
bles  required  to  run  the  terrain  aodel  are  also  included  at 
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the  end  of  this  appendix.  Professor  Jaaes  K.  Hartman  [Bef. 
14]  provides  the  details  of  the  Methodology  for  line  of 
sight  and  elevation  using  the  STAB  terrain. 

equivalent  STAB  variables  are  in  (  ) 

A 1  =  computing  parameter  for  a  hill  (A) 

B1  =  computing  parameter  for  a  hill  (B) 

CR  =  critical  value  parameter  for  a  hill  (CHIT.  H) 

Cl  -  conversion  of  orientation  angle  to  radians 
C2  =  (spread  of  a  hill)  squared 

PI  =  value  of  the  elevation  of  a  point  due  to  a  hill  (FI) 

H  (I# J)  -  the  J-th  parameter  for  hill  I;  the  parameters  are 

1  =  x-coord  of  the  center  of  the  hill  (XC. H) 

2  *  y-coord  of  the  center  of  the  hill  (TC.H) 

3  *  peak  elevation  of  the  hill  (PEAK. H) 

4  =  orientation  angle  of  the  hill  (ANG.  H) 

5  =  eccentricity  of  the  hill  (ECC.H) 

6  =  spread  of  the  hill  (5PBD.H) 

7  =  max  height  of  the  hill  (HT.H) 

8  =  cut  height  of  the  hill  (COT.H) 

H  (I« 9)  *  status  of  hill  I; 

0  =  hill  I  has  not  been  checked  for  LOS 
1  =  hill  I  has  been  checked  for  LOS 
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HL(1,J,K)  =  the  number  of  the  K-th  hill  in  grid  square  I , J 
which  influences  the  terrain  in  that  grid 
square  I  =  1  to  L;  J  *  1  to  W;  K  =  1 
to  HL  (I, J,0) 

HL(I,J,0)  «  the  total  number  of  hills  in  grid  square  I,J 

which  influence  the  terrain  in  that  grid  square 
HL  (0,0,0)  a  the  base  elevation  of  the  terrain  map 
H2  =  temporary  variable  identifying  the  hill  being  checked 
for  line  of  sight 

L  a  number  of  grid  squares  on  the  terrain  map  in  the 
x  direction  (NGRIDX) 

LX  a  x-coord  of  lower  left  grid  square  of  terrain  map  in 
1000  meter  units  (X.LO.BDHI) 

LI  =  y-coord  of  lower  left  grid  square  of  terrain  map  in 
1000  meter  units  (X.LO.BDHI) 

M$  =  name  of  files  containing  terrain  data 

NH  =  total  nuaber  of  hills  on  the  terrain  map  (NHILLS) 

PI  *  computing  parameter  for  a  hill  (PXX.H) 

P2  *  computing  parameter  for  a  hill  (PXX.H) 

P3  *  computing  parameter  for  a  hill  (PYY.H) 

QI  *  quadratic  function  of  the  ellipse  represnting  the 

hill  (QI) 
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S  =  number  of  grid  squares  on  the  terrain  sap  in  the  Y 


direction  (NG8IDY) 

X  =  current  x-coord  of  unit  whose  elevation  is  being 

coaputed  (in  1 00.  meter  units)  (X) 

XB (I)  *  current  x-coord  for  BLUE  unit  I 
XB(I)  =  current  x-coord  for  BED  unit  I 
XS  =  distance  in  x  direction  between  unit  location  and 
center  of  hill  mass  (XS) 

XI  -  relative  number  of  grid  squares  in  x  direction  of  a 

BLUE  unit  location  with  respect  to  the  lower  left 
grid  square 

X2  *  relative  number  of  grid  squares  in  x  direction  of  a 
BED  unit  location  with  respect  to  the  lower  left 
grid  square 

13  =  number  of  grid  squares  in  x  direction  between  a  BLUE 
unit  and  a  BED  unit 

Xh,Y4  -  current  grid  square  being  checked  for  line  of 
sight  in  LOS  subroutine 

Y  *  current  y  coord  of  unit  whose  elevation  is  being 
computed  (in  100  meter  units)  (I) 

YB (I)  »  current  y-coord  for  BLUE  unit  I 
IB  (I)  =  current  y-coord  for  RED  unit  I 
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IS  =  distance  in  y  direction  between  unit  location  and 
center  of  hill  aass  (IS) 

11  =  relative  number  of  grid  squares  in  y  direction  of  a 

BLOB  unit  location  with  respect  to  the  lower  left 
grid  square 

12  =  relative  nuaber  of  grid  squares  in  y  direction  of  a 

RED  unit  location  with  respect  to  the  lower  left 
grid  square 

73  =  nuaber  of  grid  squares  in  y  direction  between  a  BLUB 
unit  and  a  RED  unit 

Z  =  temporary  variable  for  the  elevation  at  a  point  on  a 
hill  (Z) 

ZB  (I)  =  elevation  of  BLUE  unit  I 
ZR  (J)  =  elevation  of  RED  unit  J 
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AR2S2Q1Z  S 


SYNOPSIS  OF  TIE  POSTPROCESSOR 
The  following  appendix  is  a  summary  of  the  work  done  by 
other  students  as  a  part  of  a  joint  project  for  OA  4655. 

The  authors  gratefully  acknowledge  the  work  done  by  Captains 
Aabrose  R.  Hock  and  Steven  L.  Maddox  in  developing  this 
postprocessor  to  reduce  the  aaount  of  work  required  in 
extracting  the  necessary  data  froa  STAR.  Their  work  is 
included  here  in  this  appendix  since  it  is  an  integral  part 
of  the  par  me  ter  estiaation  process. 

The  postprocessor  developed  for  this  thesis  is  a  Sim- 
script  language  program  that  provides  to  the  user  a  summary 
of  the  data  generated  from  the  Simulation  of  Tactical  Alter¬ 
native  Responses  (STAR)  Model.  Additionally,  the  postproces¬ 
sor  performs  some  routines  that  assist  the  user  in  the  data 
analysis  of  the  model.  The  STAR  Model  is  a  Simscript  lan¬ 
guage  program  that  simulates  combat  between  two  combined 
arms  teams  in  a  combat  environment  that  includes  field 
artillery  fire  and  electronic  warfare. 

Since  the  model  assumed  in  the  project  only  required 
data  generated  by  a  tank  vs  tank  battle,  it  was  necessary  to 
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change  certain  aspects  of  the  STAS  model.  Thus  it  was 
necessary  to  "turn  off”  the  field  artillery  and  Electronic 
Warfare  modules  within  the  STAB  aodel.  Also  the  non-tank 
units  needed  to  be  deleted.  Referring  to  figure  19,  the 
Electonic  Warfare  code(EWCODE)  was  deleted  through  the  use 
of  a  coaaent  card,  and  a  Electonic  Warfare  off  (EWOFF) 
subroutine  was  inserted. 


//STARTHES  JOB  (3102, 0234) ,* COMB  ARMS  BATTLE ', CLASSIC 
//SIM. SYSPRINT  DD  DUMMY 


//SIM. SYSIN 

DD 

// 

DD 

//* 

DD 

// 

DD 

// 

DD 

// 

DD 

// 

DD 

// 

DD 

// 

DD 

// 

DD 

DISP=SHR , 
DISP=SHH, 
DISP-SHR, 
DISP=SHR, 
DISP=SHR, 
DISP=SHR, 
DISP=SHR, 
DISP— SHR, 
DISP=SHR, 
DISP=SHH, 


DSN* 

DSN= 

DSH- 

DSH- 

DSH= 

DSN= 

DSN* 

DSN* 

DSN* 

DSN* 


MSS, 

MSS 

MSS, 

MSS 

MSS, 

MSS. 

MSS 

MSS, 

MSS 

MSS 


S3102. 
,  S3 102 . 
S3  102. 
S3102. 
S3102 . 

,  S3 102. 
,  S3102 . 
S3 102 . 
S31Q2. 
,  S3 1 02  . 


THESIS 

THESIS 

THESIS 

THESIS 

THESIS 

THESIS 

THESIS 

THESIS 

THESIS 

THESIS 


(WRKPREAM) 

(SNAPDOTR) 

(EWCODE) 

(MOVCOORD) 

(COHMWORK) 

(TEMPCOMM) 

(TEMPGRND) 

(EWOFF) 

(FAEXPER) 

(ARTYMESG) 


Figure  19:  Electronic  Warfare  Deletion 


Thus,  whenever  the  main  program  called  for  the  usage  of  EW, 
the  EWOFF  subroutine  automatically  returned  the  program  back 
to  the  point  where  the  EW  routine  was  called  and  the  simula¬ 
tion  continued  with  no  Electronic  Warfare  generated.  The 
field  artillery  within  the  STAR  model  was  cancelled  by 
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having  the  siaulation  create  FA  fire  at  tiae  equal  to  5000 
tiae  units (see  Figure  20) . 

Schedule  a  FABEGIN  in  5000.0001  tiae  units 

Figure  20:  Field  Artillery  Deletion 

Since  the  siaulation  was  scheduled  to  run  only  2500  tiae 
units,  no  field  artillery  was  generated.  Finally,  the 
internal  logic  of  the  STAB  aodel  was  used  to  '•destroy”  the 
non-tank  units  within  the  siaulation  by  changing  their  basic 
load  to  zero (see  Figure  21).  This  action  created  situa¬ 
tions,  where  the  tank  eleaents  did  not  consider  the  non-tank 
eleaents  as  lethal  or  dangerous  targets  and  thus  were 
ignored. 
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Figure  21:  Alteration  of  Basic  Load 


Output  froa  the  STAB  model  can  be  divided  into  three 
categories: 
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Update  list  of  attacker/defender  status. 

Shot  data 

Significant  event  listing  (i.e.  movement,  EH  and  FA 
events) 

Only  the  shot  data  from  the  STAfi  model  is  needed  to  provide 
the  necessay  data  for  computations.  Figure  22  is  a  sample 
of  this  output.  In  order  to  utilize  this  data,  the  internal 
output  control  of  the  model  was  used  to  place  the  shot  data 
in  a  mass  storage  location  within  the  computer  to  be  later 
manipulated  by  the  postprocessor. 

The  flow  chart  for  the  postprocessor  is  contained  in 
Figure  23;  a  copy  of  the  program  itself  can  be  found  in 
Appendix  F  as  Listing  I.  The  program  flow  is  very  basic  and 
needs  little  explanation.  From  the  mass  storage  location  of 
the  computer,  the  postprocessor  reads  pertinent  information 
into  the  SHOTLIST.  Once  all  of  the  data  has  been  read,  the 
program  processes  the  data  into  a  casualty  list (CASLIST) 
keying  from  a  status  of  "dead"  from  the  Shotlist.  Figure  24 
is  a  sample  of  the  casualty  list  output.  Once  CASLIST  is 
created  the  postprocessor  performs  histogram  computations 
and  graphs.  The  postprocessor  can  print  two  types  of 
histograms.  Figure  25  contains  the  Simcsript  Library  pro¬ 
gram  for  the  histogram.  Additionally,  the  HISIG  routine  in 
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Figure  22:  STAS  Output  Shotlist 
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FORTRAN  can  be  accessed  to  produce  a  histograa  of  the  data. 
After  the  histograa  data  is  generated,  the  postprocessor 
calculates  the  attrition  rate  coefficients  using  the  formu¬ 
las  discussed  in  Section  III-E  of  this  thesis. 

Stief.  description  o£  the  postprocessor  (see  this  appen¬ 
dix  for  variable  description): 


Preamble 

• 

lines 

1  -  9 

definition  of  permanent  entities. 

lines 

10  -  39 

definition  of  temporary  entities. 

lines 

40  -  54 

definition  of  data  to  be  collected  for 

histograms. 

Main : 

line 

2 

variable  definition. 

line 

3-4 

system  parameters  read. 

lines 

13  -  31 

creation  of  shotlist,  reading  data  from 

file. 

lines 

43  -  89 

creation  of  casualty  list,  file  "dead". 

lines 

90  -  125 

processing  of  histogram  data 

lines 

125  -  186 

processing  of  attrition  rate  coeffi- 

cients 

• 
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SHOT  LIST  DATA  FOB  CASUALTIES  FOB  THIS  STAB  SIMULATION 


MOM 

TIME 

TARGET  SIS 

-HPN 

TYPE  TIME  BETWEEN  CASUALTIES 

1 

322 

59 

1-7 

0 

2 

334 

46 

1-7 

12 

3 

46  8 

58 

1-7 

134 

4 

574 

28 

1-1 

106 

5 

581 

42 

1-7 

7 

6 

591 

27 

1-1 

10 

107 

1210 

72 

1-7 

6 

108 
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HOBBES  OF  BLOB  KILLED  *  28 

BL80  *  47,25000  BLSIG  *  299.259 

NOMBER  OF  BED  KILLED  =  84 

BDHO  =  13.24096  RDSIG  =  219.106 


HISTOGBAM  OF  BLOE  1HD  BED  CASULTIES 
(DATA  BASED  ON  TINE  BETWEEN  BLUE  CASOLTIES 
AND  TIDE  BETWEEN  BED  CASOLTIES) 
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Figure  25;  Siascript  Histograa 


liSl  2l  laiiabj.es  F2i  £ostE£Qcg§sor 
AHAT  =  estimate  of  the  attrition  coefficient  for  BED 
attriting  BLOE 

BHAT  -  estimate  of  the  attrition  coefficient  for  BLOE 
attriting  RED 

BLHST  -  Siascript  routine  for  the  histogram  of  BLOE  time 
between  casualties 

BLHO  =  Siascript  routine  for  calculating  the  mean  time 
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between  casualties 


BLSHTHST  =  Siascript  routine  for  the  histogram  of  BLUB 
shot  times 

BLSHTIG  =  siascript  routine  for  calculation  of  the 
standard  deviation  of  BLUE  shot  times 
BLSHTMU  *  Siascript  routine  for  calculation  of  the  mean 
of  BLUE  shot  times 

BLSIG  =  Siascript  routine  for  calculation  of  the 

standard  deviation  of  BLUE  tiae  between  casualties 
BLUE. SHOTS  =  time  of  the  previous  BLUE  shot 
BLUE.TIHE  =  time  of  the  previous  BLUE  casualty 
CAS  *  casualty 

CASHST  *  Siascript  routine  for  the  histogram  of 
tonal  time  between  casualties 
CASI3  =  Siascript  routine  for  calculation  of  the 

standard  deviation  of  the  total  tiae  between 
casualties 

CASLIST  =  array  containing  list  of  casualties 
CAS MU  *  Siascript  routine  for  calcualtion  of  the  mean 
for  the  total  time  between  casualties 
CBLUE  =  number  of  surviving  BLUE  force  elements 
CK  *  integer  used  in  the  calculation  of  the  attrition 
coefficient 
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1  =  BLUE  casualty 
0  =  BED  casualty 

CBED  =  number  of  surviving  BED  force  elements 

COBTIME  =  time  of  previous  casualty 

FHAME  =  firer's  naae 

FSWTYPE  *  firer’s  system-weapon  type 

J.BL  =  number  of  BLUE  casualties  (used  in  BLHST  routine) 
J.BD  =  nuaber  of  BED  casualties  (used  in  BDHST  routine) 

J BLHST  =  nuaber  of  BLUE  shots(used  in  BLSHTHST  routine) 
JCAS  a  nuaber  of  casualties  (used  in  CASHST  routine) 
JBDSHT  *  number  of  BED  shots (used  in  BD5HTHST  routine) 
JSHOT  =  number  of  shots (  used  in  SHOTHST  routine) 

JTOT  -  nuaber  of  casualties  (used  in  TOTHST  routine) 

LIM  *  array  of  time  periods  that  represent  phases  of 
the  battle 

LIMA  »  lower  time  limit  of  phase 
LIMB  =  upper  limit  of  phase 
LP  =  partial  derivative  of  p 
LQ  *  partial  derivative  of  g 
MAXBLUE  =  initial  BLUE  force  size 
MAXHED  3  initial  BED  force  size 
MIHP  =  ainiaum  of  p 
MIHQ  a  minimum  of  g 
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BK1  =  number  of  RED  elements  in  time  interval 


NK1  =  number  of  BLUE  elements  in  the  time  interval 
NLIB  »  number  of  time  intervals  to  be  evaluated 
HSHOT  *  number  of  shots  in  the  total  battle 
STABLE  *  number  of  tables  (battles)  to  be  evaluated 
SOB  =  counter  fo  number  of  casualties 
NUB3L.CAS  =  number  of  BLUE  casualties 
NOBRD.CAS  *  number  of  RED  casualties 
MOBCAS  =  casualty  entity  number 
P  *  probability  of  not  acqiring  a  target 
PHAT  =  estimated  p 
PROVAL  =  process  variable 
0  =  stop 
1  =  process 
QHAT  *  estimate  of  g 
RANGE  =  range  of  firer  to  target 

RDHST  *  Simscript  routine  for  the  histogram  of  RED 
casualt ies 

RDBO  =  Simscript  routine  for  calulation  of  the  average 
time  between  RED  casualties 

RDSHTHST  *  Simscript  routine  for  the  histgram  of  total 
red  shots 

RDSHTMU  =  Simscript  routine  for  calculation  of  the  average 
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red  shot  tine 


RDSHT5IG  *  Siascript  routine  for  the  calculation  of  the 
standard  deviation  of  the  RED  shot  tine 
RDSIG  *  Siascript  routine  for  calculation  of  the 

standard  deviation  for  HED  tiae  between  casualties 
BSD. SHOTS  «  tiae  of  previous  RED  shot 
RED. TIME  =  tiae  of  previous  RED  casualty 
ROHNO  =  run  number 
SHOT  =  shot 

SHOTHST  =  Siascript  routine  for  total  shots  fired 
SHOTLIST  =  array  of  shot  times 

SHOTSfJ  =  Siascript  routine  for  calculation  of  the  average 
tiae  of  shots 

SHOTSIG  =  Simscript  routine  for  calculation  of  the 
standard  deviation  of  total  shots 
STATUS  -  status  of  shot 
STATUS. CAS  =  status  of  casualty 
T.BTHH.CAS  =  time  between  casualty 
T.CAS  =  casualty  number 

TBLUE  =  number  of  BLUE  casualties  within  the  interval 
TGT.CAS  =  target  casualty 
TIME  =  tiae  of  shot 
TIHECAS  *  tiae  of  casualty 
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TIHELIM  =  naxtine  of  battle  when  a  casualty  occurred 

TMCAS  =  tine  of  previous  casualty 

THRDCAS  =  tine  of  previous  RED  casualty 

TOTAL. CASTIHE  *  total  casualty  tiae 

TOTAL. SHOT. TIES  =  total  shot  tine 

TOTAL. TIHE  *  total  casualty  tiae  for  the  creation  of 
T.BTWH.CAS 

TOTHST  =  Sinscript  routine  for  the  histrograa  of 
total  casualty  tiae 

TOTHO  =  Sinscript  routine  for  calculation  of  the 
average  casualty  tine 

TOTSIG  =  Sinscript  routine  for  calculation  of  the 

standard  deviation  of  the  total  casualty  tiae 
TRED  =  total  RED  casualties  within  the  tiae  interval 
TSOB  *  tiae  counter 
TSWTYPE  =  target  system-weapon  type 
TTBLOE  =  total  BLUE  casualties  for  XX  battles 
TTRED  =  total  RED  casualties  for  XX  battles 
TTCAS  =  total  casualties  for  XX  battles 
TTT  =  nuaber  of  shot 

TYPE. CAS  *  systea  weapon  type  of  the  casualty 
YAR  =  attrition  coefficient  variables 
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LISTING  Q£  COapaTBS  gHOGEAMS 


A.  DABSTAC  HQOEL  PHOGEA3 

10  REM  UPDATED  24  FEB  32 

20  REM  DETERMINISTIC  .Q66RE6ATED  MODEL  OF  STOP  ON  THE 

21  REM  APPLE  COMPUTER  <  OAHSTQC > 

30  DIM  X8<  5  :> ,  V6<  5  > ,  XRr  5  > ,  VW  5  > ,  8  V'<  5 . 1 8  > ,  8W  5 , 1 8  > ,  RX<  5  ,  1 8  > 
48  DIM  8P<  5  >  *PP<  5  > »S8< 5  > *SR<  5  > »RV<  5  •  1 0  ' 

50  OIM  8F<  5  },RF<  3  >,8T<  5 .5 >,RT<  5 >3 >,OB<  5  ".,OB<  5  > 

60  OIM  8R<  5  >,RR<  5  > 

70  DIM  81<3>*RU3> 

30  T  *  0;OT  a  30:  MT  a  1500 
30  D*  *  CHRS  <4> 

100  P  *  8.89:0  a  0.33 
U0  RI  *  3. 14 15-92654: PP  *  0 

129  A  *  2.4:8  *  1.2 

130  BE  »  0.3: RE  a  0.3 
140  BM  a  .3: PM  a  .5 
150  KK  *  0 

160  NK  *  MT  /  0 T 

16 1  REM  ZZaFLPG  FOR  TYPE  OUTPUT 

162  REM  ZZal  3  IVES  UNIT 

163  REM  LOCATIONS  fc  SIZES 

164  REM  Z2=0  SIDES  TIME  t,  TOTAL 

165  REM  RED  %  BLUE  FORCE  LEVELS 

166  44  3  0 

170  OIM  00<  NK  > »R0<  MK > 

160  HOME 

131  REM  R£ AO  DATA  FROM  FILE  FOR 

132  REM  BLUE  FORCE. 

183  REM  DATA  INCLUDES: 

134  REM  MaNUMSEP  OF  UNITS 

185  REM  QlaNUMBER  OF  COORD  RTF 

186  REM  6X»8Va»?RI0  OF  COORD  PT 
137  REM 

130  INPUT  "NAME  OF  FILE  ffOR  BLUE  HATO?  11  ;ns 
200  PRINT  OS i "OPEN  " ; NS; “.BLUE" 

210  PRINT  OSi“REAO  “  ;N*;*\BUJE" 

220  INPUT  M 

280  FOR  {  a  1  TO  M 

240  INPUT  Bi 

250  FOR  J  *  i  TO  St 

260  INPUT  B<«  l  .J  •:  INPUT  BV<  I \"> 

270  NEXT  :  NEXT 

280  PRINT  OS  {"CLOSE  NS,-" .  BLUE" 
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28 1  REM  SET  INITIAL  UN T T 

282  REM  LOCATION  EQUAL  TO 
233  REM  1ST  COORD  PT. 

284  REH 

298  FOR  I  =  1  TO  H:X8<  I  >  =  8XC  !  >1  >:V8<  I  =  =tV<  T  #1  ):RR<  I  >  a  l 
338  PRINT  81  s  FOR  JJ  =  1  TO  Bis  PRINT  BXC  I,jj  >,BVC  NEXT  JJ 

310  PRINT  “SIZE  OF  BLUE  UNIT  TNPMT  SBC  I  ' 

328  BUI)  »  SBC  I  :• 

321  REH 

322  REM  SUM  INITIAL  FORCE  SIZE 

323  REM 

338  88<  0  >  *  B0<  8  >  +•  SBC  I  > 

331  REM 

332  REM  RATE  OF  BLUE  MOUEMENT 

333  REM  EQUALS  0 

334  REM 

340  8R<  I  >  »  8 

339  NEXT 

360  PRINT 

361  REM 

362  REM  READ  DATA  FROM  FILE  FOR 

363  REM  RED  FORCE. 

364  REM  SAME  FORMAT  AS  BLUE 

365  REM 

370  INPUT  “NAME  OF  FILE  FOR  RED  DATA'?  “  :NS 
388  PRINT  OS ; “OPEN  " ;N*;u.REO“ 

398  PRINT  OS; “READ  ";NSs“.REO“ 

400  input  N 

+10  FOR  I  3  i  to  N 

428  INPUT  R1 

+30  FOR  J  3  i  TO  R1 

440  INPUT  RXC  I  .J):  INPUT  RVC  I  ,.J  ) 

+50  NEXT  :  NEXT 

+60  PRINT  OS; “CLOSE  ";N*;".RED" 

+70  FOR  I  *  i  TO  NsXRC  T  >  »  PXC  1,1  UVP<  I  ■  3  pv<  1 . 1  :•?  RPC  I  ^  a  \ 

+80  PRINT  "SIZE  OF  RED  UNIT  ";I INPUT  SRC  p. 

+98  Rlc  I  )  3  SRC  i  i 

300  R0C  0  >  a  ROC  0  >  +  SRC  r  :■ 

310  PRC  I )  3  5 
328  NECT  I 

530  PRINT  :  INPUT  "DO  VOU  WANT  THE  PRINTER  ON?  ■  V.'N  >  “.-AS 

535  IF  AS  3  -V"  THEN  PR  =  1 

540  HOME  :  UTAB  3 

5-50  IF  PR  s  t  THEN  PPM  1 

568  PRINT 

564  REM 

565  REM  PRINT  INITIAL  DATA 

566  REM 

578  PRINT  "TIME";;  POKE  36.7:  PR  TNT  "*»».UF";:  POKE  36.28:  PRINT  *PEO" 
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590  pr#  a 

-184  PEM 

585  REM  WHICH  OUTPUT  TVPE? 

588  REM 

598  IF  22  *  1  THEM  60SUB  1788?  60TQ  818 

888  60SUB  1498 

818  REM  MfilN  PR06RRH 

828  T  *  T  +  OT 

824  REM 

825  REM  CHECK  IF  Mflx  TIME 
828  REM  IS  EXCEEDED 

827  PEM 

838  IF  T  <  =  MT  THEN  848 

835  PRINT  "END  OF  8PTTLE  DUE  TO  TIME1* ;  60SUB  1540s  FHO 
648  60SUB  798:  REM  MOVEMENT  SUBROUTINE 
858  I30SU8  1170s  REM  ATTRITION  SUBROUTINE 
654  REM 

855  PEM  INCREMENT  THE  T THE 

658  REM  PERIOD  INDEX 

657  PEM 

888  KK  »  Kk  +  t 

678  88<  KK  >  a  0:R0<KK  >  a  0 

674  REM 

675  REM  SUM  SLUE  FORCE  UEUEL 

678  REM  FOR  THIS  TIME  PERIOO 

677  PEN 

698  FOR  I  a  i  TO  M 

698  88<  KK  )  a  80<  KK  )  SB<  I 

708  NEXT  I 

704  REM 

705  REM  SUM  RED  FORCE  LE"EL 

706  REM  FOR  THIS  TIME  PERIOO 

7 97  REM 

710  FOR  J  a  i  TO  N 

720  R0<  KK  >  a  R0f  kk  .<  *  SR*  j 

730  NEKT  J 

734  REM 

735  REM  WHICH  TVPE  OF  OUTPUT*? 

< i8  PEM 

740  IF  22  a  I  THEN  S09UB  1768:  on tci  -pp 
750  60SIJB  1490 
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752  REM 

753  REM  CHECK  IF  ENO  OP  PATTI  P 

754  REM  -  CRITERIA  POP  FORCE 

755  REM  LEUEL3  73  EXCEEDED 

756  REM  IF  CRITERIA  IS  EXCEEDED 

757  REM  SOTO  SAUE  RESULTS 
756  REM 

786  IF  80<KK:>  >  8E  4  B6<0>  THEN  770 

785  PRINT  :  PRINT  - ENO  OF  BATTLE  DUE  TO  ATTRITION  OF  BLUE" 

786  60SU8  1540?  ENO 

770  IF  R0<  KK  >  >  RE  4  R0<0>  THEN  730 

775  PRINT  :  PRINT  "END  OF  BATTLE  DUE  TO  ATTRITION  OF  RED- 

776  GGSUB  1540:  END 
730  SOTO  826 

736  REM  MOUEMENT  SUBROUTINE 
794  REM 

735  REM  FOR  RED  UNITS 
7:36  REM 

300  FOR  I  »  1  TO  N 
304  REM 

385  REM  SET  LAST  COORD  PT 

306  REM 

310  J  *  SRI  > 

313  REM 

314  REM  COMPUTE  COMPONENT?  OF 

315  REM  DISTANCE  BETWEEN  NEXT 

316  REM  COORD  PT  AND  CURRENT 

317  REM  LOCATION 

313  REM 

320  DV  »  RV<  I,J  M)-  VR<  I  •> 

330  OX  a  RX< I.J  +  l >  -  XRC I > 

333  REM 

334  REM  00  SEOMETRV  TO  SET 

335  REM  CORRECT  AN6LE  <  AN  > 

336  REM 

340  IF  OX  a  0  ANO  OV  >  0  THEN  ON  =  PT  -  'SOTO  330 

350  IF  OX  »  0  ANO  OV  <  0  THEN  AN  *  7  -  PT  7:  SOTO  330 

360  AN  a  HTN  ( OV  OX  > 

370  IF  OX  <  0  THEN  AN  =  PI  4  AN 

373  REM 

374  REM  COMPUTE  OH,  DISTANCE 

375  REM  TO  8E  TRAVELLED  IN 

376  REM  THIS  TIME  PERIOD 

377  REM 

330  0 1  a  OT  4  RR< I >  *  COS  <  AN  • 

330  02  a  OT  4  RR<  I  >  *  S IN  -  AN  ) 

300  OH  a  ?QR  <01  *  2  +•  02  *  2"< 
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904 

REM 

905 

REM 

COMPUTE  01.  DISTANCE 

Jw 

REM 

TO  NEXT  COORD  PT 

•907 

REM 

FROM  CURRENT  LOCATION 

908 

REM 

910 

01  » 

SQR  (OX  *  2  *  OV  -  2> 

911 

REM 

912 

REM 

IF  01  <  OH,  HOLE  UNIT 

913 

REM 

ONLV  THE  AMOUNT  TO 

914 

REM 

NEXT  COCRO  PT 

915 

REM 

IF  01  IS  USED.  UPOATE 

916 

REM 

last  cooro  PT 

917 

REM 

920 

IF  01  >  3  OH  THEN  950 

930  01  a  OX: 02  *  OV 

9*0  RR  I  >  *  j  +  i 

958  VR<  I  >  a  TNT  (.  VR<  T  >  *  02  :• 

900  XR<  1)3  I  NT  <  XR<  I  :•  +■  Ot 
970  NEXT 

974  REM 

975  REM  00  THE  SOME  FQR  BLUE 

976  REM 

980  FOR  I  *  1  TO  M 

■985  IF  BR(  1)30  THEN  1150 

990  J  »  0P<  I  ) 

1000  OV  »  8V<I,J  +  l)  -  V8<  I  ) 

1810  OX  *  bx^  I  ,j  +  i  -  :<B<  I  > 

1020  IF  OX  »  0  UNO  OV  >  0  THEN  ON  -  PT  /  ?•  ROTO  1000 

1830  IF  OX  *  0  ilNO  OV  <  0  THEN  ON  »  3  *  PT  ✓  2:  60TO  1060 

1040  AN  a  HTN  <0V  "  OX  ) 

1050*  IF  OX  <  0  THEN  UN  *  PI  +  W 
i860  01  *  OT  4  9R<  I ')  ■»  COS  -:0N  > 

1070  02  *  OT  *  8P(  I  4  8IN  <:  AN> 

1089  OH  a  sqr  (  01  '  2  +•  02  A  2  > 

1090  01  3  3QR  OX  ■'  2  ♦  OV  -  2> 

1100  IF  01  >  3  nn  THEN  1130 

1110  01  3  OX: 02  3  OV 

1120  BP'  I  :•  3  .j  m  i 

1130  V6<  I  >  a  [NT  VB<  I  >  4  02 

u«*0  xe<  i  >  *  i nt  <  :<p<  n  *  oi> 

1150  NEXT 
■160  RETURN 
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1170  REM  ATTRITION  SUBROUTINE 

1172  REH 

1173  REH  ZERO  OUT  ATTRITION 
U74  REH  UGLUES  FOR  wLL  UNITS 
1175  REH 

1 18id  FOR  I  *  1  TO  M:8F<  I  >  »  8sOB<I>  *  0*  NEXT 

1190  FOR  J  *  1  TO  N*.RF<J>  a  9fQRCJ>  *  0!  NEXT 

1192  REM 

1193  REM  CHECK  R0N6E  8ETHEEN 

1194  REM  OPPONENTS.  IF  RfiNGE 

1195  REM  <»  27O0  ATTRITION  HILL 

1196  REM  OCCUR.  REDUCE  ROTE  OF 

1197  REM  RED  MOUEMENT  TO  2  M.'S 

1198  REM 

1208  FOR  I  *  1  TO  M 

1210  FOR  J  *  1  TO  N 

1220  8T<  I  #  J  >  »  0sRT<J.I>  *  O 

1230  R6  *  SQR  <  <  X8<  I  >  -  XRC.JV)  -  2  *  'V«*I>  ~  VR<  J  V>  - 
1248  IF  R6  >  2708  THEN  1310 

1250  RR<  J  >  »  2 

1251  REM 

1252  REH  INCREMENT  THE  NUMBER 

1253  REM  OF  TOTS  FIRED  AT  0V 

1254  REM  EACH  UNIT.  COMPUTE 

1255  REM  ATTRITION  CAUSED  8V 

1256  REM  RED  ON  BLUE  fc  MICE 

1257  REM  UERSA 

1258  REM 

1279  QF<  I  >  *  8F<  I  >  +  1 

1280  RF<  J )  a  RF«:  j >  +  t 

1285  IF  RRtJ  )  a  -9  THEN  8TC  I  ,  J  >  =  Of  ROTO  1295 

1290  BTCl.J)  a  a  /  N  *  Cl  -  <P>  $B<  I>>  *  SR<  J  1 

1295  IF  9R<  I  >  a  3  THEN  RT<J,I  >  a  0.  ROTO  1310 

1300  RT(  J.I  >  *  8  /  M  4  <1  -  <Q>  *  SR<  J  V>  *  S8<'  I 

1318  NEXT  :  NEXT 

1311  REM 

1312  REM  SUM  TQTRL  ATTRITION 

1313  REM  FOR  EACH  UNIT 

1314  REM 

1328  FOR  I  *  1  TO  M 
1338  FOR  J  a  i  TO  N 

1340  IF  RF<  J )  a  0  OR  0F<O  *  0  TH*N  iT70 
1350  Oea>  *  06<  I  >  +  8T-:  I  ,J  >  ✓  RF-r.J'- 
1360  OR<  J )  a  0R<  J >  RT<J*I  )  '  9F<  t  > 

1370  NEXT  :  NEXT 
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1371  REM 

1372  REH  COMPUTE  NEW  UNIT 

1373  REM  STRENGTH.  SUBTRACT 

1374  wEH  UNIT  ATTRITION  eronM 

1375  REH  CURRENT  UNIT  STRENGTH 

1376  REH  IF  UNIT  STRENGTH  ROEG 

1377  REM  NEGATIUE,  MAKE  IT 

1378  .  REM  EQUAL  8. 

1378  REM 

1380  FOR  I  «  1  TO  rt 

1390  SBC  I  >  »  SB<  I  #  -  OBC  I 

1400  IF  3B<  I  >  <  8  THEN  SR'  T  =  0 

1410  IF  38C  I  >  <  *  BM  *  8 1C  I  >  THEN  8RC  I  >  =  S 

1429  NEXT 

1430  FOR  J  »  1  TO  M 

1440  SRC  J  >  *  SRC  j  >  -  0R<  JT 

1450  IF  SRC  J  'i  <  0  THEN  SRCJ  '»  =  8 

1460  IF  SRC  J-)  <  a  RM  *  Pt'  n  TWFN  oo<:  .*  >  «  -  * 

1470  NEXT 
1400  RETURN 

1490  REM  PRINTOllTxSUMMORV 
1500  IF  PR  a  1  THEN  PR#  1 
1510  PRINT  Tj:  POKE  36.7;  PRINT  (tftfvto; 

1515  POkE  36.20s  PRINT  R0CKK} 

1320  PR#  0 
1530  RETURN 

1540  REM  SPUE  RESULTS  IN  A  FTl.c 
1550  PRINT 

1555  INPUT  “00  vou  WANT  0  F tlf  nr  tuts  pn^'c  nqrp?  <v.^n  -  ,l  :us 
1560  IF  A#  a  "N*  THEN  1750 
1570  0#  a  CHRS  <4> 

1580  INPUT  “NAME  OF  FILE  FOR  FORCE  '  'EL  DATA'?  “;F* 

1590  PRINT  OS .“OPEN  “ jF*;“.C0EFS“ 

1600  PRINT  0S;“HRITE  “  ;FS;**.C0EFS* 

1610  PRINT  As  PRINT  8s  PRINT  P:  PRINT  0 
1620  PRINT  OSi “CLOSE  “ ;FS:“.COEFS“ 

1630  PRINT  OSi “OPEN  “ ;FS RESULTS “ 

1640  PRINT  OS;“HRITE  “  ;FS,-“.  RESULTS" 

1650  PRINT  NK  *  l:  PRINT  NIC  +  1 

1660  FOR  I  a  0  TO  NK 

1670  T  a  I  4  .3 

1680  PRINT  T:  PRINT  SAC  I  > 

1690  NEXT  I 

1700  FOR  I  *  A  TO  NIC 

1710  T  a  I  *  3 

1720  PRINT  T:  PRINT  R0C I- 

1730  NEXT  I 

1740  PRINT  OSi “CLOSE  “ :F*;“. RESULTS" 

1750  RETURN 
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17S8  REM  PRINTOUT/SUMMRRV 

1779  HOME  :  UTRB  3 

1739  IF  PR  *  1  THEN  PR#  1 

1739  PRINT  “TINE  IS  ";T 

i see  PRINT  :  PRINT 

1810  PRINT  "8LUE  FORCES" 

1320  PRINT  “UNIT" ;  TfiS<  7  >i“XCn0R0" ;  TOR<  13  >suVCOnRO“ i  TQ8<  23>f"SI2E" 
1830  FOR  I  ■  1  TO  M 

1340  PRINT  “  “;I ;  TP8<  7>?XB<  I  >;  TQ8<  !5'jV8<  !  >:  TQ8<  24>;SB<I> 

1850  NEXT 

i860  PRINT  :  PRINT  “RED  FORCES" 

1370  PRINT  "UNIT";  TR8<  7>j"XC00R0“ ;  TQRC  i5 >;“VCOORn" :  tqb<  23>;"3I2E" 
1380  FOR  I  *  l  TO  N 

1390  PRINT  -  ";I;  T08<  ?>;XR<  I  >;  TfiB<  15>;VR<I>;  TQgK  2# ';SR<  I  > 

1900  NEXT 

1918  PRINT  :  PRINT 

1820  PR#  8 

1330  INPUT  “HIT  RETURN  TO  CONTINUE" ;Zt 
1940  RETURN 
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FOHCS  BAKER  PBOGRAB 


I a  REM  FORCE  OATA-FTLE  MAKER 
28  DIM  6X< 5 *  18> #RV<  5  » 1  «:■  ,81  <  18  ' 

38  SOTO  178 
HOH£ 

58  INPUT  “NAME  OF  FILS’  “  ;NS 
88  OS  *  CHRS  4  :> 

.79  PR  I  NT  OS;  "OPEN  "  ;  NS;  “ .  “  ;  MS 
88  PRINT  OS; "WRITE  “;NS;“. ";NS 
88  PRINT  M 
188  FOR  I  a  1  TO  M 
118  PRINT  8t< I > 

128  FOR  J  »  1  TO  81U  ) 

138  PRINT  8X<I,J>:  PRINT  8Y'T,J' 

148  NEXT  J 
158  NEXT  I 

168  PRINT  OS ; “CLOSE  “ ;NS ; “ . “ ;MS 

178  HOME 

188  UT08  5 

198  PRINT  “MENU; “ 

208  PRINT  “1  MAKE  FILE  FOR  BLUE  FOprES" 

218  PRINT  “2  MAKE  FILE  FOR  RED  FORCES- 

228  PRINT  “8  QUIT- 

238  PRINT  ;  INPUT  “WHICH’  " tZ 

248  IFZM  THEN  6NO 

•250  IF  Z  ■  l  THEN  MS  a  “BLUE" 

268  IF  Z  =  2  THEN  MS  »  “RED- 

270  HOME 

290  PRINT  “WHAT  IS  THE  FORCE  SIZE  OF  “;MS;:  INPUT  "  “;M 

298  FOR  I  =  1  TO  m 

380  PRINT  :  PRINT  “FOR  UNIT  ":I 

318  PRINT  “HOW  MANY  COOROINOtinO  POINTS  " ; 

315  INPUT  “<  INCLUOTNS  START  LOCATIONS  •’  "  ;PH  I  • 

320  FOR  J  »  1  TO  61'.  I  :• 

338  PRINT 

340  PRINT  “ENTER  OR 10  COORDS  FOR  COORD T NOT INS  00 TNT  *  “ 
350  INPUT  "X-OOORD  IS  ";SX<  I ,  J  ■ 

360  INPUT  “V-COORD  IS  “  I  ,J  > 

378  NEXT  j 
380  NEXT  I 
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•■330  HONE 

400  PRINT  j  PRINT  TB8<  t2 >;“#**  OBTO  CHECK  ***“ 

410  PRINT  :  PRINT  :  PRINT 

420  PRINT  “UNIT  PT";  TR6*  9 >;“X-COOPO“ ;  TO0<  t? >;“V-COORO" 

430  FOR  I  *  l  TO  M 
440  FOR  J  *  1  TO  91(1  > 

450  PRINT  TQ6<  2>ili  t«8<  TRB<  19  >sRX<  T ,.j  >;  TQB(  13>;8V<I,J> 

400  NEXT  J 
4?0  NEXT  I 

*80  INPUT  “RNV  CHRN6E3'?  (V-'N>  “  ;Z* 

490  IF  2*  *  “N“  THEN  40 
300  PRINT 

510  INPUT  “WHICH  UNIT'5’  “;I 
520  INPUT  “WHICH  POINT'?  “ i J 
530  INPUT  “X-COORO  IS  “;BX(t,J> 

540  INPUT  “V-COORO  IS  “;8V<I, J  ) 

550  SOTO  390 
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C.  RESULTS  READER  P HOGS AH 


10  HOME  :  UTwS  3 
20  0*  »  CHR*  <  4 

>0  INPUT  “NOME  OF  FILE  FOP  FORCE  LS'-'EL  ORTQ*  ”  ;F * 
40  PRINT  0* • "OPEN  “;F#;“.COEFS“ 

50  PRINT  0*»“REflO  “ ;Ff;“. OOEFS” 

60  INPUT  0:  INPUT  9:  INPUT  P:  INPUT  Q 
79  PRINT  0*i “CLOSE  " sFf ;“.COEFS“ 

90  PRINT  0*; “OPEN  “ sFs i“. RESULTS “ 

90  PRINT  0*;“REP0  " ;F*;“. RESULTS” 

100  INPUT  N:  INPUT  N 

110  NK  a  N  -  1 

120  OIH  80< NK > ,R0< NK  > 

1 30  FOR  I  a  0  TO  NK 
140  INPUT  T:  INPUT  80<  I  > 

150  NEXT  I 

160  FOR  I  a  0  TO  NK 

i?e  input  t*  input  R0<t> 

130  N0<T  I 

100  PRINT  0*; “CLOSE  ”;F*;“. RESULTS” 

200  INPUT  “PRINTER  ON?  ”;«S 
210  HOME  :  UTPR  3 
220  IF  a*  a  “V"  THEN  PR#  1 
230  PRINT  -«  a  «;* 

240  PRINT  “8  =  '  ;6 

250  PRINT  “P  a  “ ;p 

260  PRINT  “0  a  -jQ 

270  PRINT 

230  PRINT  “TIME” ::  POKE  36.7?  PRINT  •‘SLOE” 

235  POKE  36.20:  PRINT  "RED" 

230  FOR  KK  a  0  TO  NK 
300  T  *  J0  +  Kk 

>10  PRINT  T;:  POKE  36. 7:  PRINT  80000; 

315  POKE  36.20:  PRINT  R0<  KK  :■ 

320  NEXT  KK 
730  PR#  0 

3-*0  ENO 
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HIP  SAK28  PBOGSAM 


10  REM  MOP  CPTOLOS  MPKER/REPOER 
20  DIM  HL<10*10»12> 

30  HOME  :  UT06  3 
40  PRINT  “MENU: " 

50  PRINT  “i  MOP  MAKER- 
60  PRINT  -2  MOP  READER" 

79  PRINT  “0  QUIT” 

30  PRINT  :  INPUT  "WHICH'7'  "iX 
30  IF  X  *  0  THEN  ENO 
100  ON  X  SOTO  120*300 
110  REM  MOP  CATALOG  MAKER 
120  HOME  :  UTQB  2 

130  PRINT  T08<  11  >**•##*  MOP  MAKER  **+* 

140  PRINT  :  PRINT 

150  INPUT  “NEW  OR  OLD  MOP'7  <  N/O  "■  “  s2* 

160  IF  Z$  =  "N"  THEN  II  =  lsJJ  =  1:  SOTO  220 
170  INPUT  “NOME  OF  MOP'7  “;M* 

130  0 S  a  CHfii  ( 4  '! 

130  30SUB  620 

200  INPUT  “ENTER  STOP  TINS  SR  10  SWORE  'Z,!>  “:r!  *JJ 
210  SOTO  250 

220  INPUT  "SIZE  OF  MOP  •:  1000  M  SOUOPES  v?  “;L*W 

230  INPUT  "ENTER  LOWER  LEFT  6RI0  SQUARE  <EG.  50*33'  "*LX*LV 

240  INPUT  “BOSE  ELEVATION  OF  MOP^  “;HL< 0*0*0 > 

250  FOR  I  *  II  TO  l 
260  FOR  J  a  jj  TO  W 
270  PRINT  :  PRINT 

230  PRINT  “FOR  SR  10  SQUARE  “  ;  I  ;“.“*  J 

290  INPUT  “NUMBER  OF  HILL  MOSSES'7  ";HL<  I *J*0  >:KK  a  HL<  I,J*0> 
300  IF  KK  a  0  then  380 

310  PRINT  “ENTER  THE  “;KK;“  HILLS  POP  RRIO  -WQPF  ;J 

320  FOR  K  *  1  TO  KK 

330  PRINT  “HILL  ";K*“  IS  INPUT  ul<I*J.k> 

340  NEXT  .< 

350  JJ  =  1 

360  PRINT  :  INPUT  "SOUS'7  >  V/N 
370  IF  Hi  =  "V"  THEN  400 
330  NEXT  J 
390  NEXT  I 
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400  HOME  :  UTA8  4 

419  INPUT  "NAME  OF  MAP  FILS'7  "  ;MS 

420  os  a  CHRS  <:4> 

430  PRINT  OS i "OPEN  "*MS;“.MAP" 

440  PRINT  OS; “WRITE  “ ;MS*“ . MAP" 

450  PRINT  L:  PRINT  W 
460  PRINT  LX:  PRINT  LV 
470  PRINT  HU  0*0*0  > 

400  FOR  I  »  1  TO  L 
460  FOR  J  =  1  TO  W 
500  PRINT  HL<I#J.0> 

510  IF  HL<I*J*0>  »  0  THEN  550 
520  FOR  K  »  1  TO  HL<I*J.0> 

538  PRINT  HLU,J*K  > 

540  NEXT  K 
550  NEXT  j 
560  NEXT  I 

570  PRINT  OS; "CLOSE  "*MSs".MAP“ 

560  PRINT  :  PRINT  :  PRINT  :  PRINT  "...  DONE":  °R I NT  :  »R 

580  INPUT  "DO  VOU  WANT  TO  CONTINUE  INPUTTING  HATA^  <  V.-N  > 

600  IF  II  a  "V*  THEN  200 

610  SOTO  30 

620  OS  =  CHRS  <4> 

630  PRINT  OS; "OPEN  ";MS;".MAP" 

640  PRINT  OS; "READ  " iMS;".MAP“ 

650  INPUT  L;  INPUT  W 
660  INPUT  LX:  INPUT  LV 
670  INPUT  HL< 0*0*0 > 

630  FOR  I  *  i  TO  L 
630  FOR  J  a  1  TO  W 
700  INPUT  HL-.  I*J*0> 

710  IF  HU I*J.0>  a  8  THEN  750 
720  FOR  K  =  1  TO  HL<I.J*0> 

730  INPUT  HL<I*0*JO 
740  NEXT  K 

730  NEXT  J 

760  NEiT  I 

770  PRINT  OS; "CLOSE  ,';MS;".MAP- 
730  RETURN 

730  REM  MAP  CATALOG  FILE  READER 

300  HOME  :  UTAB  2 

610  PRINT  TA8<  11  ■>;"***  MAP  PFQOER 
620  PRINT  :  PRINT 
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330  INPUT  “NRME  OP  HfiP‘  FILE’  " ; HS 
*40  INPUT  "PRINTER  ON'?  <VxN>  ";flS 
350  0*  *  CHftf  <4> 

360  PRINT  0* i "OPEN  " iMSj“.MRP" 

370  PRINT  USi"REftO  “  iM*;*4 .  MRP- 
380  input  l:  input  h 
•390  INPUT  LX:  INPUT  LV 
300  INPUT  HL< 0.0.0 > 

310  FOR  I  »  1  TO  L 
920  FOR  J  »  i  TO  H 
930  INPUT  HU  I.J.0) 

■340  IF  HL<I.J.0>  =»  0  THEN  980 
350  FOR  K  »  1  TO  HLU.J.0) 

:360  INPUT  HUI#J.K) 

370  N&XT  K 

380  NEXT  J 

990  NEXT  I 

1800  PRINT  QS; "CLOSE  “;MS;“.MGP“ 

1010  HOME 

1028  IF  RS  *  «v-  THEN  PR#  1 

1030  PRINT  "NAME  OF  MRP;  "iMS 

1040  PRINT  "LGHER  LEFT  i3RI0  S0UPRE;  "  ;l.X  "  *LV 

1050  PRINT  "SIZE  OF  MRP;  ";L;“  8V  ";W 

1060  PRINT  "8«SE  ELEU0TION:  -;HL< 0.0.0 > 

1070  PRINT 
1080  FOR  I  *  1  TO  L 

1090  FOR  J  »  1  TO  W 

1100  PRINT  r,J.HL< I.J.0* 

1110  IF  WUC I . J .0  >  *  0  THEN  1160 
1120  FOR  K  »  l  TO  HL<I.J.0> 

1130  PRINT  HLv I.J.K 

1 140  NEXT  K 

1150  PRINT 

1160  NEXT  J 

1170  PRINT  j  NEXT  I 

1130  PR#  0 

1130  PRINT  ;  PRINT 

1200  INPUT  "HIT  RETURN  TO  CONTINUE" sXS 
1210  8OT0  30 
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HILL  HAKES  PROGBAH 


>j  «£m  hill  ONTO  >=ILE  MWKEP.'REfiPEP 
20  QIM  rK  200*3 > »K<  3 > 

30  K<  1  >  *  5s  K< 2 >  »  9iK<3>  3  l4;K<-*:<  »  ?3 
40  K<.5>  *  22iK<6>  *  24*K<?>  *  29sK<3>  =  33 
50  REM  HILL  OftTR  FILE  MRKER 
30  HOME  :  UTOe  2 
7®  PRINT  “MENU:** 

30  PRINT  "1  HILL  MPKER" 

30  PRINT  “2  HILL  RE&OER* 

100  PRINT  “0  (SUIT" 

110  PRINT  :  INPUT  “WHICH'7'  “;X 

120  I F  X  *  0  THEN  ENO 

130  ON  X  SOTO  150 >320 

140  REM  HILL  OfiTR  PILE  MRKER 

150  HOME  :  UTRB  2 

160  PRINT  TR8<  ll>;“444  HILL  MQKER 
170  PRINT  ;  PRINT 

130  INPUT  “NEW  OR  OLD  PILE'7  •'  N/O  '  M:0* 

130  IF  W  a  “N“  THEN  23® 

200  SOSU8  560 

210  INPUT  “WHICH  HILL  00  VOU  WONT  T0  STQRT  WTTW7  “;II 
220  SOTO  2b® 

230  HOME  :  UTwB  4 
240  II  *  1 

250  INPUT  “NUMBER  OF  HILLS'7  “  ;N 
260  for  I  a  II  TO  N 

27®  PRINT  "ENTER  THE  OflTO  POR  HILL  “;t 
280  INPUT  “XCOGRO  IS  “;H<.  I,i> 

230  INPUT  “VCOORO  IS  “ ;H< I ,2 > 

300  INPUT  "PERK  ELECTION  IS  “  fH(  1 .3  > 

31®  INPUT  “ORIENTRTION  0N6LE  IS  “ ;H*  T .4  • 

320  INPUT  “ECCENTRIC ITV  IS  “iW  1,5' 

7.30  INPUT  "3PREP0  IS  “  iH<  1 ,6  > 

340  INPUT  “MwX  HEISHT  IS  “  iH<  1,7  ".' 

7.50  INPUT  “CUT  H£  I  SWT  IS  “,*wri,3> 

36®  PRINT  :  PRINT  :  INPUT  '’SOUS'7  <  V.- N  >  “  ,ZS 

379  IF  2*  *  “V-  THEN  390 

380  NEXT  I 

730  HOME  :  UTQB  4 
400  0*  a  CHR*  '  4  > 

410  input  -nohe  of  fils’7  “;n* 

420  PRINT  0* ; "OPEN  “ *NA .HILLS,  L34“ 

430  PRINT  0*;u  WRITE  “ ;N*:“. HILLS,  R®“ 

44®  PRINT  N 

450  FOR  I  a  1  TO  N 

■US®  PRINT  0*; “WRITE  “  ;N*s“ .  HU  L3,  ;T 
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470  FOR  J  *  1  TO  3 
480  PRINT*I*J> 

498  NEXT  J 
500  NEXT  I 

510  PRINT  OS; "CLOSE  “ ;NS;“. HILLS" 

520  PRINT  :  PRINT  .-  PRINT  "...  OONE" 

525  PRINT  :  PRINT 

530  INPUT  "00  VOtJ  WONT  TO  CONTTNUE  INPUTTING  OPTO*?  <V/N>  ";PS 


540 

IF  PS  3  “V"  THEN  210 

550 

GOTO  69 

560 

HOME 

570  OS  *  CHRS  <  4  > 

590 

INPUT  "NPME  OF  FILE?  " ;NS 

590 

PRINT  OS; "OPEN  “ ;NS;". HILLS- 

L34" 

600 

PRINT  OS; "REPO  ";NS;". HILLS* 

R0" 

610 

INPUT  N 

620 

FOR  I  *  1  TO  N 

630 

PRINT  OS; "REPO  “;NS;“. HILLS* 

R"  ;  T 

640 

FOR  J  =  1  TO  8 

650 

INPUT  * I*J  > 

660 

NEXT  J 

670 

NEXT  I 

690 

PRINT  OS; "CLOSE  ";NS;“. HILLS" 

690 

HONE 

7O0 

PRINT  "NUMBER  OF  HILLS:  " ;N: 

PRINT 

710 

M  *  l 

720 

FOR  I  *  1  TO  N 

739 

IF  H  a  11  THEN  PRINT  :M  =  1 

740 

PRINT  I; 

750 

FOR  J  3  1  TO  3 

760 

PRINT  TP6<  K<  J  )  );H<  I*.J  <s 

779 

NEXT  J 

730 

PRINT 

790 

M  3  M  l 

300 

NEXT  I 

310 

RETURN 
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320  REM  HILL  QOTP  FILE  REPOER 

330  PR  »  0 

340  HOME  :  UTP8  2 

350  PRINT  TP8<  10);“*#*  HILL  REPOER  ***“ 

360  PRINT  :  PRINT 

379  0*  *  CHRS  '4> 

380  INPUT  “NPME  OF  MPP'?  “;NS 

338  INPUT  “PRINTER  ON?  -'V/'N')  “;2# 

380  IF  2*  »  “V"  THEN  PR  «  t 

310  PRINT  OS;“OPEN  “;NS;“. HILLS,  1.34“ 

320  PRINT  OS; “REPO  “;NS;“. HILLS,  R0“ 

330  INPUT  N 

340  FOR  I  »  1  TO  N 

■350  PRINT  OS; “REPO  “;NS;“. HILLS,  R“;I 

360  FOR  J  a  1  TO  3 

370  INPUT  H< I , J  > 

380  NEXT  J 
330  NEXT  I 

1900  PRINT  OS; “CLOSE  “ ;NS;“. HILLS- 
1810  HOME 

1020  IF  PR  a  1  THEN  PR#  1 

1030  PRINT  s  PRINT  NS;“  HILL  OPTR“:  PRINT 

1040  PRINT  “NUMBER  OF  HILLSi  N 

1050  PRINT 

i860  M  a  l 

1070  FOR  I  a  1  TO  N 

1080  IF  H  a  tt  THEN  PRINT  :M  a  1 

1080  PRINT  I; 

tl00  FOR  J  a  t  TO  0 

1118  KK  a  k<J>  ♦  3  *  J 

1120  IF  PR  »  t  THEN  POKE  TR.KKj  PRINT  GOTO  11*3 

1130  PRINT  TQBc  K<J)>;H<  I,.J); 

1 140  NEXT  j 
1150  PRINT 
1160  H  a  H  *  t 
1170  NEXT  I 
1130  PR#  0 

1 1'30  INPUT  “HIT  RETURN  TO  CONTINUE"  «XS 
1200  GOTO  60 
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?.  aOODS  MAKES  PB0GBA3 


L0 

REM  LOGOS  CATALOG  MAKER/REAOER 

28 

OIM  FU 

10.10.3) 

30 

HOME  s 

UTA6  3 

40 

PRINT  “ 

MENU:" 

58 

PRINT  “ 

'1  WOODS  MAKER" 

60 

PRINT  “2  WOODS  READER" 

70 

PRINT  " 

■9  QUIT" 

80 

PRINT  : 

INPUT  “WHICH'?  "iX 

80 

IF  X  » 

0  THEN  ENO 

100 

ON  X  GOTO  120.770 

110 

REN  WOODS  CATALOG  MAKER 

120 

HONE  ! 

!  UTA8  2 

130 

PRINT 

TA8<  11 >;“***  WOODS  MAKER  *** 

140 

PRINT 

:  PRINT 

150 

INPUT 

“NEW  OR  OLD  MAP-?  <  N.'0  )  “  .-?* 

168 

IF  ZS 

a  “N"  THEN  IT  a  j* jj  a  1:  GOTO 

170 

INPUT 

"NAME  OF  MAP?  ";MS 

180 

03  =  CHR*  <4) 

130  G03UB  600 

800  INPUT  "ENTER  STORTING  BRIO  SWOPE  <I,J>  I.JJ 
210  GOTO  240 

220  INPUT  "SIZE  OF  MOP  <1000  M  SQUARES  V?  “:L,W 

230  INPUT  “ENTER  LOWER  LEFT  GRID  SQUARE  <  EG.  50.33)  "iLX.LV 

240  FOR  I  *  II  TO  L 

250  FOR  J  »  JJ  TO  W 

280  PRINT  j  PRINT 

270  PRINT  “FOR  GRID  SQUARE  "  ;  I ;  ’ V ;  J 

280  INPUT  “ NUMBER  OF  FORESTS'?  " ;FL<  I.J.0  ):KK  »  FL'il.J.Q) 

280  IF  KK  »  0  THEN  ?70 

S00  PRINT  “ENTER  THE  “;KK;“  FORESTS  FOR  ARID  SQUARE  “;t;“.";J 
310  FOR  K  a  l  TO  KK 

320  PRINT  "FOREST  “;K;“  IS  INPUT  FLU.J.K* 

■330  NEXT  K 
340  JJ  a  1 

350  PRINT  ;  INPUT  "SAUE'?  <V'N'";A* 

360  IF  A*  a  "V"  THEN  390 
370  NE<T  j 
380  NEXT  I 
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396  HOME  :  UT08  4 

466  INPUT  “NONE  OF  MRP  FILE?  “;H* 

416  0*  a  CNR*  <  4  ) 

426  PRINT  OS; “OPEN  “ ;M*;" .WOODS" 

436  PRINT  0*J “WRITE  “  .  WOODS* 

440  PRINT  L:  PRINT  W 
456  PRINT  LX:  PRINT  LV 
466  FOR  I  »  1  TO  L 
470  FOR  J  *  1  TO  W 
466  PRINT  FL<I#Jj>0> 

496  IF  FL< I»J#0>  »  9  TWEN  530 
500  FOR  K  »  1  TO  FL< I »J»9  > 

516  PRINT  FL< I # J#K  > 

520  NEXT  « 

536  NEXT  j 

546  NEXT  t 

550  PRINT  0*i“ CLOSE  "  sW*;'* .WOODS* 

560  PRINT  :  PRINT  :  PRINT  :  PRINT  *'...  DONE":  PRINT  :  PRINT 

570  INPUT  "00  VOU  W«NT  TO  CONTINUE  INPUTTING  OPTO?  < V  -N )  “ tZS 

580  IF  2*  a  “V"  THEN  200 

590  GOTO  30 

606  OS  a  CHRS  <  4 ) 

616  PRINT  0*i“OPEN  *;H*j'\ WOODS- 

620  PRINT  0*j“R£flO  WOODS'* 

630  INPUT  L:  INPUT  W 

640  INPUT  LX:  INPUT  LV 
630  FOR  I  a  t  TO  L 

660  FOR  J  a  i  TO  W 

676  INPUT  FL<I#J#0> 

660  IF  FL(I#J»0>  »  0  THEN  720 

690  FOR  K  a  1  TO  FL<  I.J.G  - 
700  INPUT  FLv I »J#K  > 

’10  NEXT  K 

720  NEXT  J 

T30  NEXT  I 

T40  PRINT  0*s "CLOSE  “ irtSi * . WOODS" 

730  RETURN 
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788  REM  WOODS  CPTPL06  FILE  REPOER 
778  HOME  :  OTPS  2 

738  PRINT  TP8<  1 L  >;“*+*  WOODS  READER  W 

7S8  PRINT  :  PRINT 

388  INPUT  "NAME  OF  MPP  FILE?  “  ;MS 

318  INPUT  “PRINTER  ON?  <VN>  “;PS 

320  OS  »  CHRS  <4> 

•338  PRINT  OF;  “OPEN  “  ;MS;  WOODS" 

340  PRINT  OS » “REPO  “;MSs“. WOODS* 

358  INPUT  L:  INPUT  W 

368  INPUT  LX:  INPUT  LV 
370  FOR  I  a  1  TO  L 

360  FOR  J  *  1  TO  W 

390  INPUT  FL<I*J»0> 

300  IF  FL< I»J«8>  »  0  THEN  940 
310  FOR  K  a  l  TO  FL< I,J,8> 

320  INPUT  FL< I,J,K> 

330  NEXT  K 
340  NEXT  J 
350  NEXT  I 

360  PRINT  OSi “CLOSE  “ iMSi“. WOODS* 

378  HONE 

380  TF  PS  a  -v“  THEN  PR#  t 

398  PRINT  “NPHE  OF  MPP:  “;HS 

1888  PRINT  “LONER  LEFT  8RI0  SQUARE:  “;LX;“  " ;L^ 

1910  PRINT  “SIZE  OF  MPP:  “it i"  8V  “;W 

1820  PRINT 

1038  FQR  I  »  i  TO  L 

1040  FOR  J  a  i  TO  W 

1850  PRINT  I,J,FU  I,J,0‘:. 

i860  IF  FL<  I #J»8>  a  0  THEN  tU8 

1870  FOR  K  =  1  TO  FLv  I  #J#0  > 

I860  PRINT  FL<l,J,K>j“  “i 

1890  NEXT  K 

1180  PRINT 

1110  NEXT  J 

1120  PRINT  :  NEXT  I 

1 130  PR#  8 

1140  PRINT  :  PRINT 

1158  INPUT  “HIT  RETURN  TO  CONTINUE” :XT 
11 68  SOTO  38 


157 


G.  FOHEST  SAKEB  ? BOG RAH 


10  REM  FOREST  MAKER/ READER 
20  OIH  F<  120.6 ).K<S> 

•30  K(  1 .)  *  5:K<2>  *  12:K<3>  «  2RsK?  A’'  *  2R 
40  K<5  *  30s  K<6>  *  36 
50  HOME  :  UTA6  2 
60  PRINT  “MENU:  " 

79  PRINT  “1  FOREST  MAKER- 

30  PRINT  "2  FOREST  READER “ 

30  PRINT  -0  QUIT" 

100  PRINT  :  INPUT  "WHICH?  ";X 
110  IF  X  *  0  THEN  ENO 

129  ON  X  SOTO  140,690 

139  REM  FOREST  OATA  FILE  MAKER 

140  HOME  :  UTA8  2 

150  PRINT  fR&<  10>;"***  FOREST  MAKER 
16®  PRINT  s  PRINT 

179  INPUT  "NEW  OR  OLO  FILE7  «'  N/'O  >  ";AS 

130  IF  As  a  "N‘‘  THEN  220 
1 90  GOSUB  530 

200  INPUT  "WHICH  FOREST  DO  VOU  WANT  TO  FT APT  NTTH?  * ;  H 

210  GOTO  250 

220  HOME  !  UTA8  4 

230  II  »  t 

240  INPUT  “NUMBER  OF  FORESTS'?  “;N 
250  FOR  I  a  ii  TO  N 

260  PRIHT  "ENTER  THE  OATA  FOR  FOREST  “;T 
279  INPUT  "XCOGRO  IS  "  ;Fv  I ,  i:> 

259  INPUT  "VCOORD  IS  “;F<  1,2 '• 

290  INPUT  “HEIGHT  OF  TREES  IS  "  ;F*  1 ,3  > 

200  INPUT  "ORIENTATION  AN6LE  IS  ":F«'l.4> 

310  INPUT  "SEMI-MAJOR  axis  LENGTH  IS  " ;F-  1 .5  ■ 

•220  INPUT  "SEMI-MINOR  AXIS  LENGTH  IS  "  ;P'  r  ,p  > 

330  PRINT  s  PRINT  ;  INPUT  ,*SAUE'?  <V/N>  "  .Z* 

240  IF  2*  *  "V"  THEN  3S0 
350  NE>T  I 

260  HOME  :  UTA8  4 
370  OS  *  CHRS  <  4  :■ 


380  INPUT  “NOME  OF  FILE?  11  ;NS 

390  PRINT  OS; “OPEN  NS .FORESTS,  L3S“ 

490  PRINT  0*;H  WRITE  “; NS; “.FORESTS,  RR“ 

410  PRINT  N 

4^0  FOR  I  »  1  TO  N 

430  PRINT  0S;“WRTTE  “ ;NS;“. FORESTS,  R“sT 
440  FOR  J  »  t  TO  8 
458  PRINT  F< I ,  J  > 

460  NEXT  J 
470  NEXT  I 

*00  PRINT  0*; “CLOSE  NS; “.FORESTS* 

490  PRINT  i  PRINT  :  PRINT  DONE" 

495  PRINT  :  PRINT 

500  INPUT  “00  VOU  HfiNT  TO  CONTINUE  INPUTTING  ORTO*?  <V/N>  “  ;RS 
510  IF  OS  a  “V“  THEN  2130 
520  'SOTO  59 
530  HONE 

540  OS  a  CHRS  <4> 


550 

INPUT  “NOME  OF  FILET  “ ;NS 

560 

PRINT  OS ; “OPEN  NS; “.FORESTS, 

139“ 

570 

PRINT  OS;-R£AO  NS; “.FORESTS, 

R0“ 

580 

INPUT  N 

590 

FOR  I  a  i  TO  N 

600 

PRINT  QS;“REPG  “;NS;“. FORESTS, 

R“:I 

610 

FOR  J  a  i  TO  S 

620 

INPUT  F< I , j  > 

630 

NEXT  J 

640 

NEXT  I 

S50 

PRINT  OS; “CLOSE  “;NS;“. FORESTS 

W 

660 

HONE 

670 

PRINT  NS;*  FORESTS  LIST";  PRINT 

680 

PRINT  “NUHBER  OF  FORESTS:  “;N: 

PRINT 

690 

N  a  1 

700 

FOR  I  a  1  TO  N 

710 

IF  H  a  11  THEN  PRINT  :N  a  i 

720 

PRINT  I; 

730 

FOR  .j  a  i  TO  6 

740 

PRINT  TP8<  K<  J  >  );F<  I  ,.J  • ; 

750 

NEXT  J 

760 

PRINT 

779 

M  a  M  *  1 

790 

NEXT  I 

790 

RETURN 
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;■ m  REM  FOREST  OfiTP  FILE  PEfiOER 
310  HONE  :  UTflB  2 
320  PR  *  3 

330  PRINT  T08<  10 >,'“***  FOREST  REPOER  «#*“ 

340  PRINT  :  PRINT 
330  0*  *  CNR*  <4> 

360  INPUT  “NONE  OF  MPP^  “;N* 

370  INPUT  “PRINTER  ON^  <V/N>  “ iZt 

880  IF  Z*  *  “V“  THEN  PR  *  1 

398  PRINT  0*i“OPEN  “  ;N*j**.  FORESTS.  L7P" 

900  PRINT  0*i“REP0  “iN#;'*.  FORESTS.  P.0“ 

310  INPUT  N 

920  FOR  I  *  1  TO  N 

330  PRINT  Q*#“REPO  “iN*;“. FORESTS.  P“;I 
940  FOR  J  *  l  TO  6 

350  INPUT  F<  I.J> 

360  NEXT  J 

370  NEXT  I 

380  PRINT  0*; “CLOSE  “ ;N*j“. FORESTS" 

390  HOME 

i860  IF  PR  a  1  THEN  PR#  1 

1010  PRINT  s  PRINT  NS."  FOREST  OPTO**:  PRINT 

1020  PRINT  "NUMBER  OF  FORESTS  “;N 

1030  PRINT 

1040  Hal 

1850  FOR  I  a  1  TO  N 

1060  IF  H  a  11  then  PRINT  :M  a  1 

1078  PRINT  I; 

1080  FOR  J  »  1  TO  S 
1990  KK  »  K<J>  ♦  3  *  J 

1180  IF  PR  *  1  THEN  POKE  36.KK:  PRINT  F/I.J>»s  GOTO  tl20 
1110  PRINT  Tfifi<  K<  J>>*F<  I.J>* 

1 120  NEXT  J 
1130  PRINT 
1140  H  a  H  +■  1 
1150  NEXT  I 
l 160  PR#  0 

1170  INPUT  “HIT  RETURN  TO  CONTTNUF" 

1 130  GOTO  50 
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3.  TEBBAIH  SO  DEL  BBOGHAfl 


18  REM  UPDATED  82  ccB  32 
REH  OAMSTAC  TERRAIN  MODEL 

DIM  XB<  5  >,V8<  5  ),XR<  5  > .VR<  *5  >,RX<  9.  1 3  ■  .8V<  5,  t8  '  .PX<  5, 10  > 
DIM  2B<5>»ZR<  5  >.BL'-  18,18, 12 > 

0 IM  8P<  5  > »RP<  5  ••  ,S8<  5  VSR<  5  ) ,  RV'  9 . 1 8  > 

OIM  BF<  5>.RF<  5  >.8T<  5,5>,RT<  5 .5  >.0B(  5 >,DR<  5  • 

DIM  8R<  5  >.RR<  5> 

OIM  81'  5  >»R1<  5  ) 

T  »  8s OT  *  30s  MT  »  1500 
OS  *  CHRS  <4> 

•  P  »  8.39:Q  a  8.33 
110  PI  a  3. 141592654: PR  a  0 
120  A  a  2.4:8  a  0.03 
130  BE  »  8. 3s RE  a  0.3 
140  BM  *  ,5:RM  a  .5 
150  KK  a  0 
160  NK  a  rtT  /  OT 
165  22  a  1 

I?8  OIM  80*NK>,R0<N»O 

171  HOME 

172  INPUT  "NAME  OF  MOP  FILE’  H;MS 

173  PRINT  OS OPEN  “;HS;“.MOP* 

1 74  PRINT  OSi-REAO  “;MS;*.MOP* 

175  INPUT  L:  INPUT  A 

176  INPUT  LXs  INPUT  LV 

177  INPUT  ML <:  0.8,8  ) 

173  FOR  I  »  l  TO  L 
173  FOR  J  a  1  TO  W 

180  INPUT  HLf.  1,4,8 '> 

181  IF  HL< 1.4.8)  a  0  THEN  134 

182  FOR  K  »  1  TO  ML1'  I.J.0  ) 

183  INPUT  HU  1. 4, JO 
134  NEXT  :  NEXT  :  NEXT 

185  PRINT  OSi “CLOSE  “.'MS .-".MAP** 

|  136  LX  a  LX  -  IsLV  a  LV  -  t 
137  PRINT  OS i “OPEN  " ;MS;“ .HILLS,  L34" 

188  PRINT  OS s “READ  “ :MS»". HILLS*  PR" 

139  INPUT  NHs  OIM  hKNH,3> 

130  FOR  I  a  l  TO  NH 

191  PRINT  OS : “REAO  “  ;MS;‘\ HILLS,  R“ ;l 

132  FOR  J  *  l  TO  3 

193  INPUT  H<I,J> 

134  NEXT  j  NEXT 

1 35  PRINT  OS; “CLOSE  ";MSs».HTLL3" 
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— i 
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I 98  HOME 

;99  INPUT  “NAME  OF  FILE  FOR  BLUE  QATQ?  - ;N* 

200  PRINT  OS i “OPEN  “ ;NS;“. BLUE “ 

210  PRINT  DS;"READ  “ iNS;“ , SLUE" 

220  INPUT  M 

230  FOR  I  *  1  TO  H 

240  INPUT  81 

250  FOR  J  *  1  TO  81 

260  INPUT  8X< I #J  >:  INPUT  BVC I #  J  > 

270  NEXT  :  NEXT 

290  PRINT  OS, “CLOSE  “;NSi“.BLUE“ 

290  FOR  I  *  1  TO  Ms  X8<  I  >  *  BXC  I  #  1  >s  V8C  I  >  =  BV<  I  #  1  >s  BP‘  I ')  =  1 
300  PRINT  81:  FOR  JJ  a  I  TO  91:  PRINT  9X<  I  #.jJ  >  ,8V<  I  #  JJ  >:  NEXT  JJ 
310  PRINT  “SIZE  OF  SLUE  UNIT  ";l5t  INPUT  SWT'* 

320  81<  I  >  a  S8<  I  > 

330  80<0)  «  80<  0  >  +  SBC  I  > 

>40  8Rt  I  >  *  0 

341  X  a  x8<  I  ):V  a  V8<  I  > 

342  X4  a  1NT  <  X8<  I  >  X  1000*  -  LX: ¥4  a  TNT  «*  V8<  T  >  '  1000  *  -  l_V 
>43  S0SU8  i960 

>44  i8(I)  a  2 

>45  PRINT  “ELEVATION  OF  SLUE  “;I.*“  IS  “;Z 

>46  FOR  II  a  1  TO  HUX4,V4,0>:W2  a  NL<X4,V4,Il  >:H<W2,0>  a  0:  NEXT 
>30  NEXT 
360  PRINT 

>70  INPUT  "NAME  OF  FILE  FOR  RED  OQTCr?  «  ;NS 
380  PRINT  OS i “OPEN  “;NS;“.REO" 

390  PR*  JT  OSi“REAO  “;NS;“.RED“ 

400  INPUT  N 

*10  FOR  I  *  1  TO  N 

420  INPUT  R1 

430  FOR  J  a  i  TO  R1 

440  INPUT  RX<I»J>:  INPUT  RV<  I  ,J  J 

436  NEXT  :  NEXT 

460  PRINT  OSi “CLOSE  “ #NS;“.REO“ 

470  FOR  I  »  l  TO  N:,XR<  I  :<  S  RX(  T,1  -:VR<  T  >  a  RVC  1,1  ->:RP<  I  ■  =  1 
480  PRINT  “SIZE  OF  RED  UNIT  TNPt.iT  SRC  I  - 

490  Rlc  l  :>  *  3R<  I  > 

500  R0<  0  >  a  R0<  6  )  ♦  SRC  I  > 

510  RR<  I  )  a  5 

511  X  a  :<R<  I  >:V  a  VR<  I> 

512  X4  a  INT  iXRCI  *  x  1000')  -  LX:  ¥4  a  TNT  •  vpC  I  •  x  1000  )  -  LV 

5 1 3  60SUB  1  '360 

514  2R<  I  .*  a  j 

515  PRINT  “ELEVATION  OF  RED  "ih“  IS  “  ;.T 

516  FOR  II  a  1  TO  NLCX4,V4,0  >:H2  *  MLCX4.V4. 1 1  >:WC  H2.0  >  =  0:  NEXT 
529  NEXT  I 
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■530  PRINT  ;  INPUT  “00  V00  WANT  tup  opruTPO  ON~  'V.-M7  "  ;AS 

535  IF  A*  =  "V"  THEN  PR  =  1 

54#  HOME  :  UT06  3 

550  IF  PR  *  1  THEN  PR#  t 

560  PRINT 

570  PRINT  “TIME'S:  POKE  38.7?  PRINT  '‘SLOE" ::  ^OKE  38,20:  PRINT  "RED“ 
580  PR#  O 

580  IF  32  3  I  THEN  GOSUB  1760:  GOTO  610 

800  GOSUB  1490 

610  REM  MAIN  PROGRAM 

620  T  a  T  ♦  OT 

630  IF  T  <  a  MT  THEN  640 

635  PRINT  “ENO  OF  BATTLE  DUE  TO  TIMF"*  GOSUB  1540:  END 

640  GOSUB  790 

650  GOSUB  1170 

660  KK  *  KK  4  l 

670  80CKIO  a  0;R0<  KK ')  *  0 

680  FOR  I  *  i  TO  H 

680  80<  KK  :>  a  80<  KK  >  4  6B<  I  • 

700  WENT  I 

710  FOR  J  a  t  TO  N 

729  R0<  KK  >  a  R0<  KK  >  +  SR<  J  > 

730  NEXT  J 

740  IF  12  3  1  THEN  GOSUB  17R0:  'GOTO  "60 
750  GOSUB  1490 

760  IF  80<KK  >  >  BE  *  88*  8  >  THEN  77R 

765  PRINT  :  PRINT  “ENO  OF  BATTLE  OUE  TO  ATTRITION  OF  BLUE" 

766  6GSU8  1540:  ENO 

770  IF  R0<  KK  ;•  >  RE  4  R0<0  >  THEN  780 

775  PRINT  :  PRINT  "ENO  OF  RATTLE  OUE  to  ATTRITION  Op  RED" 

776  GOSUB  1540*  ENO 
730  GOTO  6 20 

7:90  REM  MOUEHENT  SUBROUTINE 
300  FOR  I  a  1  TO  N 
610  J  »  ftFKi  > 

620  OV  a  av<  I  ,-j  ¥  \  )  -  VR<  T  > 

830  ox  »  RXf  i ,  j  +  r>  -  :<p,<  i  > 

840  IF  OX  a  0  UNO  OV  >  0  THEN  ON  =  E»T  /  GOTO  PRO 

850  IF  OX  a  0  UNO  OV  <  0  THEN  AN  =  j  ■*  PT  •  2*  GOTO  380 

860  AN  a  ATN  <  OV  /  0X> 

370  IF  OX  <:  0  THEN  AN  a  PI  *  AN 
880  01  *  QT  4  RR< I >  4  COS  < AN  > 

390  02  a  OT  4  RRtP  4  SIN  'AN':' 

980  Oh  3  30R  <01  *  2  4  02  2  7 

910  01  =  SQR  <  OX  ^  2  4  OV  -  2  > 

920  IF  01  >  »  OH  THEN  350 

930  01  »  OX: 02  3  OV 
940  RP  1  >  a  J  f  1 
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550  VR<  I  >  a  I  NT  c:  VR(  I  >  +  02  > 

960  XR<  I  >  *  INT  '  XR<  I  >  +  01  > 

961  X  a  XR<  I  >:V  a  VR<  I  > 

962  X4  *  INT  <XFKI>  1000  )  -  LX:V4  a  INT  <  VP<  T  "•  1000  )  -  LV 

963  SOSUB  I960 

964  2R< I  )  a  2 

965  PRINT  "ELEUOTION  OF  RED  “;I;“  TS  -  ;2 

966  FOR  II  a  i  TO  HUX4.V4,0  ):H2  a  WLO<4,V4.Il  )sKKH2<6>  *  0:  NEXT 
970  NEXT 

966  FOR  I  »  t  TO  H 
990  J  a  ep<  I  > 

1000  OV  a  8V<  I  ,J  +  1  >  -  V8<  I  > 

1010  OX  a  ex<I,J  +  1 )  -  XB<  I  > 

1029  IF  OX  a  0  6N0  DV  >  0  THEN  ON  a  PT  /  2*  GOTO  I960 

1030  IF  OX  »  0  ONO  OV  <  0  THEN  ON  a  3  *  PT  /  2s  SOTO  I960 

t040  ON  a  OTN  (OV  /  OX) 

1050  IF  OX  <  0  THEN  ON  a  Pi  +  £N 
i960  01  *  OT  *  BR<  I  >  *  COS  (ON) 

1079  02  *  OT  *  0R<  I  >  *  SIN  (ON!) 

I960  OH  a  SQR  (01  *  2  +  02  *  2> 

1090  01  *  3QR  (OX  a  2  ♦  ov  a  2) 

1100  IF  01  >  a  OH  THEN  1130 
1110  Cl  =  OX: 02  =  OV 

1129  Bf*  I )  »  j  +  i 

1 130  VB<:  I )  a  INT  <  VB<  n  +  02  ■) 

1 140  X8< I )  a  INT  ( XB< r  )  +  01  ) 

1150  NEXT 

1160  RETURN 

1 179  REM  OTTRITION  SUBROUTINE 
1 186  FOR  I  a  i  TO  HjBF< I >  a  0: QB< I  >  »  p:  NEXT 

1190  FOR  J  a  i  TO  N:RF<rj  )  a  0:OR(  .J  '  »  6:  NEXT 

1290  FOR  I  =  l  TO  M 

1210  FOR  J  a  i  JO  N 

1229  BT(  I #.J  »  a  0:RT<  J,I  )  a  0 

1230  R6  a  3QR  ( <  X8<  I ')  -  XR<  J  >  >  *.  3  +.  .r  VB<  T  •  -  VR<jy>  2> 

1240  IF  RG  >  2706  THEN  1316 

12 50  RR(  I)  a  2 

1251  REM  CHECK  LINE  OF  SIGHT 

1252  XI  a  fNT  ( X8<  I  >  /  1800  )  -  LX:  VI  *  TNT  >'  VB<  I  >  /  1000  )  -  LV 

1253  X2  =  INT  ( XR<  J  >  /  1000  )  -  lX;V2  a  INT  r  VR(.j  )  x  1000  >  -  LV 

1254  X3  *  X2  -  X 1 : V3  a  V2  -  vt 

1255  FOP  II  a  13  TO  OBS  <X3> 

1256  X4  *  XI  +  II  *  SON  f X3  > 

1257  FOR  Jl  =  0  TO  OBS  (  V3 

1258  V4  a  vi  +  Jl  *  SON  V3 

1261  RRIHT  "CHECKINS  SRID  SQUPRE  " ;X4;“," ;V4S“  =0R  BLUE  “;ti“  RED 
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1262  80SUS  i960 

1263  NEXT  :  NEXT 

1264  FOR  I 1  *  1  TO  NH:H< 1 1 ,6  -  *  0:  NEXT 

1265  REH  LOS  IS  OK 

1276  6FC  D»8FU  ')  +  1 
1280  RFC  J  >  *  RFC  J )  ♦  1 

1290  8T(I*J>  «  Q  »  <  l  -  CP-'  *  SB<  I>>  4  SRC  i> 

1380  RTCJ.p  »  8  *  <t  -  <0>  SRC  J  >  >  *  S0CI> 

1319  NEXT  s  NEXT 

1320  FOR  I  *  1  TO  M 
1338  FOR  J  *  1  TO  N 

1340  IF  RFC J)  »  0  OR  6FC I  )  *  O  THEN  t?7« 

1350  08c  I )  a  08C  I  >  +  9TC  I  .J  >  ✓  RFC  J  > 

1360  ORCJ)  »  ORCJ)  *  RTCJ,I  )  '  BFC  I  ) 

1370  NEXT  s  NEXT 

1388  FOR  I  a  i  TO  H 

1380  SBC  I )  »  SBC  I  )  -  OBC  I  •) 

1400  IF  SBC  I  >  <  0  THEN  SRC  I  )  a  0 

1410  IF  38c I >  <  a  BN  4  Bt< I >  THEN  SRC  I >  =  3 

1420  NEXT 

1430  FOR  J  *  1  TO  N 
1440  SRC  J  >  a  3RCJ)  -  ORCJ) 

1450  IF  SRC  J  >  <  a  THEN  SRCJ)  =  8  * 

1400  IF  SRC  I >  <  *  RM  *  RtC I >  THEN  RRC  pa  -  3 

1470  NEXT 
1480  RETURN 

1490  REH  PRINTOUPSUMMGRV 
1500  IF  PR  a  i  THEN  PR#  1 

1310  PRINT  Ts:  POKE  36*7:  PRINT  80CKK>j;  POKE  36,20:  PRINT  R0OOO 
1520  PR#  0 
1530  RETURN 

1540  REH  SPUE  RESULTS  IN  R  FILE 

1550  PRINT  :  INPUT  "00  VOU  WONT  R  FTLE  OF  THIS  RUN'S  ORTR^  CV/N)  ";HS 
1360  IF  P#  a  “N"  THEN  1730 
1570  0#  a  CHR#  C  4) 

1580  INPUT  '‘NRHE  OF  FILE  FOR  FORCE  LEUp.  OPTQT  “  ;F* 

1390  PRINT  0#i“0PEN  “ ;F*j“.COEFS“ 

1600  PRINT  D#; “WRITE  “ ;F*;“. COEFS* 

1610  PRINT  R«  PRINT  8:  PRINT  Pt  PRINT  0 
1620  PRINT  OS i“ CLOSE  •  ;F*j“.COEFS“ 

1630  PRINT  OS*“OPEN  “ iF*t\ RESULTS" 

1640  PRINT  0*; “WRITE  “ jF#s“. RESULTS* 

1650  PRINT  NK  +  1:  PRINT  HK  +  1 

1660  FOR  I  a  i?  TO  NK 

1670  T  a  r  *  3 

1680  PRINT  T:  PRINT  90C I) 

1690  NEXT  t 
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1700  FOR  I  *  0  TO  NK 

1710  T  *  I  *  3 

1720  PRINT  T;  PRINT  R€K  I > 

1730  NEXT  I 

1748  PRINT  0*; "CLOSE  " ;F*;“. RESULTS" 

1750  RETURN 

1760  REH  PR I NTOUT/SUMMGRV 

1770  HONE  s  UTPB  3 

1780  IF  PR  ■  1  THEN  PR#  1 

1780  PRINT  "TINE  IS  “;T 

1380  PRINT  :  PRINT 

1310  PRINT  "BLUE  FORCES" 

1320  PRINT  "UNIT";  TQ8<  ?>;"XCOORD“;  TO0<  i.5>s"VC0nR0“ ;  Tfi8<  23). “SIZE" 
1830  FOR  I  s  1  TO  H 

1840  PRINT  "  “  ;I ;  T08<  7>;X8<I>;  TRB<  t5>jVMI>i  TPB<  24>;SB<I> 

1850  NEXT 

I860  PRINT  s  PRINT  "RED  FORCES" 

1370  PRINT  “UNIT";  T«B<  7>;"XC00R0";  TQB<  15 > ; "VCQORO" ;  TfiB<  23>;"SIZE" 
1880  FOR  I  «  1  TO  N 

1330  PRINT  **  *;I;  T08<  7'»*XR<  I>*  TQB<  15>;VR<  I  >;  TQB<  24>;SR<I  .» 

1300  NEXT 

1310  PRINT  :  PRINT 

1 920  PR#  m 

1330  INPUT  "HIT  RETURN  TO  CONTINUE" tZS 

1  SCTI  |Dki 

1350  REM  ELECTION  COMPUTATION 
1360  CR  »  -  39393999 
1370  2  »  HL<0#0#0> 

1380  FOR  HI  »  1  TO  HL<  X4*V4#0  "> 

1390  H2  a  HL< X4,V4,Ht  ) 

2000  IF  H<  H2#0 )  »  1  THEN  2150 

2010  XS  *  X  -  H<H2*1)  *  190; VS  a  v  -  H<M2.2>  *  100 

2029  PI  *  LOS  <  H<  H2#7')  /  <H<  H2,?>  -  50 )> 

2039  Bl  *  PI  4  <H<H2#5»  *  2 
2840  Cl  *  H<H2,4>  *  PI  ✓  180 
2050  C2  *  H<H2,8>  *  2 

2060  PI  a  -  <oi  *  <  COS  <C1  •  >  A  2  4  bi  #  /  SIN  <C1>'  2>  '  C2 

2070  P2  a  -  <01  *  <  3IN  <  Cl  > )  ~  2  +  81.  *  <  COS  <01  v>  2>  C.? 

2080  P3  a  <2  *  COS  <C1>  *  SIN  <Ct  >  *  <  01  -  Pt  ,V>  y  C2 

2090  IF  H<  H2 , 7 >  H< H2 »8 )  THEN  CR  a  LOS  <<I-KH2#7>  -  H<H2*8V>  -  H<H2#7 

2100  01  *  PI  4  XS  A  2  +•  P2  4  vs  2  4  P3  *  XS  *  vs 

2110  H<:  H2>0  *  1 

2120  IF  QI  <  CR  THEN  2150 

2130  FI  a  H<  H2#3>  4  HCH2,7>  +  <  EXP  •  01  >  -  1  > 

2140  IF  FI  >  z  then  Z  a  ft 
2150  NEXT 
2160  RETURN 
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I.  POSTPROCESSOR  PROGSifl 


— 

4 

UI 

— 

VI 

a. 

X 

UJ 

> 

VI  4 

-1 

— 

■» 

UI  vi> 

VO  CO 

X 

VI 

-1  UI 

UI4 

VI 

4 

VI 

BVI  _J_) 

-J— 

CU 

a 

UJ 

<UI  B< 

BX 

« 

-I  VI 

— -J  <UI 

44 

4 

z 

BUI 

Xco  —X 

</* 

— > 

2 

4  J 

44  X 

Ui 

X 

— 

— 

—O 

>—  <1/1 

•a 

<J 

UJ 

IS 

X4 

XCO  >4 

ca 

>4 

X 

—  • 

4— 

X<UI 

< 

UI 

4  - 

col— 

>XV) 

UI>_i_JO. 

ft* 

co 

XX 

Zv> 

4 

4UI 

0  B4 

a£ 

UI 

UI 

u.  a 

U< 

X>-J 

UIX4UI  • 

< 

-J 

UM 

f- 

• 

UI  B 

—UJ— XI— 

</>  > 

B 

UI4 

44f-—  - 

044 

cOZOX  4 

Ui 

4 

•ft 

m 

> 

— 

> 

> 

UI 

> 

UI 

UI 

aui 

UI 

a 

a 

•X 

a 

• 

• 

a— 

• 

a 

0 

— 

a 

i— 

1— 

VI VI 

1— 

VI 

CO 

4 

CO 

UJUJO 

HI  LUZX  • 
tuujx  r— ♦--) 


I  .  ■■■  X  H/)H 

i-co— vi_i4a 
co  •<(•« 

<ui  •  i—O 
30000-U. 


h-  ui  H/iui/n>  lu£m  <MUi«/)x  -j  xv>— »  BQ  Of—  *x 

Zui-j  -co— *— 4t/Hu»— a.x  u  i— >4B  co  uiuix  -  -  -XZO 

—Zcq  UJ  0'*X2J<  Ult/iz  <  0-14  zzzxz«o 

00—4  x<h^«  -  -t/)uj0>  44444xui*- 

co— —  z  a  couoi—  cum  XaO—  xuif-4  uiXuiXuioxvj 

—  4  «X  <  -X  -04  V)v)ZjMWia«  4-lZ— XUIXOXOXO  — 

PI>  J<  «lUV)MI-UX4<UJd  44J>Z  >03— XLLlX  O  O  l-UJX 

—  U4>  a.  44  U  04iU  4  4  4  <0<IW-UiKtUWI 

o  LUH*  4>4UK<JUJU0liJHQ  »V)  -  X— vixuZXviX  v»X— ui 

—  xq-j  t-  uu/i  X3*-ae-K34</i«n  -<s>Lua:<  «-—«-x  X 

o  Xf-4  -30X  0>4ft-Z4  XU4Xl-UlO<  XZ  X  X  V1H- 

x  4  UJ  UlVIf—  4f-BZ<— >>X  •  <4-llJi:»LUUJ— U.VO  VI  V0UI4 

00  X  -X  X»-  f-  ctf-  «OI-VI>XBft»  XO  0<1JJ<UJ<X  vi  UJ 

—  UI  —  VI  -V1Q  V1V1X  Z04  OT4ZX— ujirt  X  X  1-04  Ul-J 

x  v«-  «x<  t-<o<utouj  -  <uiO  t-o  »  —uj  1-1—4x01-01-2  x  -jb 

u  —vO  UI—  Z  4Z3CUI  •  oui  «  ■<  4  OXZ  UJX  X  XVW/I—  B4 

a  _j—  3— VI  4  -a  -CjOZXUJKOUIXvO  »</*XI—  41/IUIO— Irtl0-Jjfl0lrth*440  4— 


• 

4  «X 

HI 

-14 

— 

U4> 

0 

UI— 

— 

XQ-I 

0 

XI-4 

X 

4  ui 

VI 

X  -X 

>— 

UI 

X 

V1»- 

•t< 

CJ 

—VI 

UI— 

a 

-J— 

O— VO 

X 

— -1 

J  *4 

a  vi 

BO 

• 

X4 

XUJUJ 

i/lij 

4XX 

r-* 

X  — 

ao 

44 

•*» 

O 

•HJJV? 

~4 

cov) 

I-X4 

m 

ZZ 

z— O 

Z  «U3 


—  a£4 

-O  -  4> 


(M  .*  j*CO_>»—  «V|»— 4 

O  OOUJO  •-JUIVIX 

•HO*  —  OUI4— 

m.J  XZHI3M-- IH- 

— oa  uiui— ui-io—  o 
mo  — >«3i-— uix 

flOOJ  vivoz  Z>— VO 

QXZ  »UJUIUIUIU|4 

-j—  vivo  zzz  x> 
vo  vi  —— —  vcux 
>*ui  uiuj  zu.  u.  xx  Zuj 

UV0-JXXUJUJUJUJ4UJ> 
acu  .cuk-kZOOOXOUJ 

QXXX  4  O 


ZZ  zviui— <uix—  xc^x  —tcjuif-xm  as  b>b  q£lu 

C4  44  O  03  mot—  — 4Z  X  XX  hX>X>XOXO  UlO 

4  4Z  *>V)  Z  co>— VIU— t-t-  -  OBOCQO  OO  Out 

4UJI— 44  —X  X<  X  VOVIZ  Z  ZOZ<U  UJt— 

Vi  viacx  «v)i_j  »  »hc  —  •<<  o  o  o  —  — z 

4  4  >  4-VJV04  «<N  -C7X  OUJO  UJOUJOUHMUJ  ZM 

x  —ujo«uuua:aj  -ZO-m-XcmxcmX  XO  -* 

hHixz  uxa.<z  *  — <z-jcaat-  i—  *-o»— t-  z 

V)  VH-CO— VOI—  «0>>— O.XXOCQXX  act-  Z4 

4  — — ILI-Ztof-Zt*  XJUZOh — VIVM-VI—VO  VOO  4 

(J  _J  —  -  24  4  404—  VO 

—  UIUI  UIUI  UUJUJUI—UJ  UIUI  UIUI  UIUJ  O  O  — ’  V04 

>  azzzz^zzzazzzzzzzuj-i— a— —voi-  4 

X  x— — 4— — 4— — — — i— cot— xi— Cvo<co  — 

U<  vOu.U_LL.U.t_HLU.U.>U-U.LLCLCL.U.CL<  «V0  •1/IKTuX  — 

>  uiuj  >1  uiuiui  uiuiujujuiujuj-joxiX'^*— |co 

UI  UlQQQOmOQQUIQQQOaOQJ^U^Q^a^J 

—  —  —  0-1  — I  — J  — J  z*— 

LL  u.  U.  O-JO-lO-lO-jOO— U.U. 

UI  UJ  UJ  U<24Z<Z4ZZ<  UIUI 

q  a  a  <i— 4>— 4t-<t-4ujxoQ 


XUI— 4 

X 

X 

ft* 

»• 

X  cOUJ 

X 

X 

u. 

u. 

u. 

NNVflt 

UI 

UJ 

UJ 

Ui 

UJ 

NNsa 

X 

a 

a 

a 

167 


o  » 

UJ  * 
QC  * 
xto* 


* 

* 

* 

* 

u 

X 

h- 

o 

z 

UJ 

at 

b- 

LO 

O 

UJ 

ac 

♦ 

♦  -j 

♦  X 

♦  — 

♦  b- 

♦  — 

♦  z 

4  — 

♦ 

♦  to 

♦  O* 

X* 


£ 


a 

QCUJ 
a  to 
UJ 
ULUU  UJ 

at  to 


S  X 

«0  —4 

»  I 
<M  -r 
X  X 
to  to 
h  >*  n  > 
X  X 
-1-3 
V--JV--J 

axax 
x  x 

to—  to- 

-3-3 

UJUJIUIU 

O.ZO.Z 
>-  >• 


to 

> 
UJX 
to  3 


a 

_  x 

30  UJ  to 

tux  — 

ZUJ  %  UJ 

O'  <>J  X 

>-•  o  < 

at  X  Z  —  3030 

<a  tu  u»  *-  to^-tot- 

h- t-  at  O  u.  tu 

to  -  -  x  o  a 

3—  to  atoafO 

a  h»  -*  o  o 

to  x-  a  uj  <n  <m 

-ja 

uj  at  a 


a 

X 

to 

UJ 

X 


a 

< 

UJ 

at 


UJ 


to 

> 

x 


<  UJ< 

-  to  to 
a.  xx 
◦  uu 
o  •  • 

-J  O-J 

-  ac<a 

XX 

a  oo 


to 


h.  x 

to  3 

—  *to 
*x 


QO 


a 

o 


X  -4-4 

a 

<  QO 
tu  ao 
Q  XX 
X 


at  n 
x 

xx 
XV- 
-O  V- 

ujz  a 
WZX  OUJ  <x 
X33  -JQC  Z 

XaCZ  «*-  -4 

X  xto 
-XZ  <  at 
~  xuj  a 

3  U. 

03  toX-JtO 

a  uiv-<ouj<ti 

UUJ  z—  X-4 
I2V-JJ-  O 
UJ— *,3  X-JH- 
UJ—  —  t» 

•  -vozv-oz 

a^*3t-Mjo« 

Z  ♦»--JZ»- 

Z>— -4.  —  ai»— 

oz-va.-*  au<ouj 

»  —  wOj 


X 

a  3  <03  <o 

— 

b-  *  O 

o  * 

to 

>  4— 4nQ  — 

o  * 

h» 

— 

h-  1  «  1  « 

a 

X  UJ 

z 

tu 

3  (Mxn> 

— 

•o  11 

-J  X 

a 

to—  3  »X  41 

ac 

— 

X  X 

z 

H-V-  z  Z 

to 

— 

Z  V- 

O  v- 

< 

•ooxox 

h- 

—  to«o 

oar 

ac 

-Xv-OV-O 

— 

3 

— o 

OUJ  at 

1—0  -J  -J 

»— 

X 

h-to_jZ 

UJ 

4> 

a— jo-jo 

O 

to 

go  — 

-V-  — 

— 

V- 

533252 

X 

to 

— 

to 

xxuju. 

tov-X 

1-0  X 
toztoui 

a 

— *-o  ar 

— 

3 

h -JC 

—  xat 

X 

ojxujv-ujj— 

UJ 

V— 

XV-  x 

jtOUO 

to 

Xv—  UJ  UJ 

X 

XI- 

tjazat 

h—  • 

— 

Xto  to -J  to -J 

— 

t— to 

xx— o 

O  —to 

uu 

Z  —*  — 

1—0— 

tuto  a 

X—  — 

X 

h-Q  h*>  h— 

V 

-  -J 

Ztoct 

tOI-Q 

OUJOX 

•odtxa 

uj  too. 
U-4— X 
a  ujui 
at  a.a.3 
a.— >» 

-  Xt-to 
<03— J 
Q  tOX 
v-atu.3 
ui  ar 


CJUJX 

t— tu—  —  x— a— a  to 

-Zatv— u-v-aujj-  xv-  x*-«m 

Q—  a  h  o<atatogto_ifrt 
UJ  tOX->  -XUJ  X  X  H— 
—  qh<jio  .—to  at  to  <o  to  to  too  u 
o  zo-J—zv—  a— — x> 

X  VUZUJUJXOUJQXUJXUJXtO  • 

-o.oxaxuja.j-Q.j-  z 

<0«0tU>— J«0>UJ— >»  >  ZX 


to 


QOae  tu  — Q.  OV-— 

XXXOcu  0X0  ZX  H- 
uitotuXro  i— Xu.  —  «o-40 

— ct— otto—  to— too  —  X 
»  UJ  II  UJJ*  j— too  to 

>toa<oz>-  v-  u.z  o3Q- 

•a>a—  -3  3ujoo  xj— <tu 

Z-JJ--J  ZUJ  coz  —  tOX  x 
Xuj  3ujj->xz  — ccv  v-  x 


to 

X 

tj 


—33  O— to— ao-  atoutat  Xto  Z  X 

.fjuvju.  3030  -uu  -iu.i-u.x-j  t-  aaa  u  h- 

xuj  to  at  v*  to  tMtooctoacv-az  O  toOO  v— Xa  xu.  uj 
H.QQ.U.  KtLHHU.hU.hQU*  OfSI  lO  UH  X  O  UJ—  X 


ac-v  *v 

at 

m 

Q — 

oa  toz 

to 

3 

X 

x  — 

UJ 

UJ 

uj  at 

Z 

at 

at  a 

m 

— 

x  XO<—  v- 
uj  ax^u-tuQ 
at  to— ox 
O—  U.U.  uu 
X 


to  toxac 


mac  -  v-zvo 


CLO.  U.UJU.  Uito  a  a. 

a_j 

atov-z  — •  ac 

— — Q— O  i—  —  r— — 

NMh0V)U)MiUZV)0 

stst  LUUJ  INUJ3C 

^tUUJ 

uj— at  XV- xu. 

toio  JJQQhVI 

to  -4 

UJO.1-  u 

—  ox 

a  je 

tu  UJ 

as  uj 

at 

a-  -  -  -  z 

—  X 
tu 
at 
o 


168 


to 

> 

< 

3 

-J 

4 


-j 

> 


-U 


tOi 


i/l 


$ 

4— 


4_j 

4 


i/l 

2  _  _ 

y-  x  a 


i/> 

< 

o 


i/> 

< 

o 

-J 

CO 

X 


l/l 

3 


t/lUJ 


o»-» 


I  </> 
I  —  -I 
i/l  < 


44*4— 4t-«/» 
— 00*000  3 
*-ui  •*— 4— tux 
OXJ 

uj* 


—  ZLLK/1  UJ 

>-434—  X 
QXuomm 


l-i-o  3  </>U 
O— at 
XUJtU  M  I_|t03*  UJ 

too.)-  uj4.J*i/> 
UJ 


< 

— 

** 

3 

t/l 

»* 

-1 

4 

** 

4 

O 

•  • 

** 

i/l  i/l 

** 

UJ 

l/l 

UJ< 

** 

X 

> 

— o 

— 

4 

1— UJ 

II  II 

3 

* 

z 

JXi/l 

• 

— 1 

* 

a 

4—3 

OO 

UJ 

4 

* 

— 

3HI 

— 

3 

/I* 

>— 

i/>  t- 

iflin 

-1 

UJ 

3* 

< 

4— 

_ia 

co 

X 

X 

-j 

um-  « 

oa  at 

l-  II 

3 

4  l/l* 

u. 

H* 

X 

Z04* 

1-0*  * 

a 

m 

i/l  l/l 

— 

UJ— O* 

3** 

a 

<4 

t/l 

UJI/f— « 

X** 

t/i 

UJ 

OO 

3<</l* 

i-** 

-j 

at 

• 

at 

1—04* 

-  ** 

4 

aa 

4 

UJXO* 

oo  *  * 

t— 

LL 

at  uj 

>- 

CQ3  •* 

— ** 

a 

i/> 

o 

Xac 

t/l 

zz 

l/ll/l** 

t- 

< 

3  t- 

UJ  X 

-IQ** 

a 

l/l 

ZX(-o 

i/l 

XXH 

caac 

UJ 

a 

-J 

3— 

— 

— *-«£> 

X 

at 

lTI 

4 

-Z-i 

X 

>■■■'«  • 

•  -I  « 

t— 

X 

X 

O 

/l  i/i 

t— 

2H 

33 

>— 

X 

a 

t/14  4 

XX33 

't— 

I 

i— 

40*0 

at 

UJUJ  » 

JQZX 

Qt/l—OJ  UJ 

u<4toX  X 
at04— ♦►■* 

S  —(_}»-•* 
l-o—  •*>- 
H  4t/lQ*UI 
g<uj*M 
— UIOOt*UJ>- 
1/1XUJ  »attO 

a*o  —  — woz*<—  4—  xx*  — 

It  ►-♦UJ— */»-«/>»- 4—  U-J  —-J 

1/1 


—  a*  — •  *»-x*—  uj»-uj<omuj  *  a 

i/lUJ||i-*Z>-ZUJ— *tOX3at<l-XX  • 

<X  3  3VH<4<<MOU  -►-*-« 
umujXikihui  <«gz>> 

•>-X>-  to  mct-J*— *-  «/*/>  3*  » 

t-  — -HJJ  —  •  <#lfl  It  — 4UI  U  » 

«  H3XH>1HhI>  »4  It  UX  uj*  toi/i 
It  l/l  —  4 
1» 


o  •* 

•  _l*Z 
Qai*  — 
OCX* 
XX  ui 
22  il  4 

z  u 


0.2— 
>— t/l 


i/l 


toac-to 

coat 

—  «. 

<«** 

oo** 


i zo*  —a** 


Z' 

a. 


/> 


— «/»OZtOIJJ—  tXVKOS  *0-111  •«/» 

Ul<  — 4_J</lUJl/lX*  —4tjx>— UJ— 1  .J4  won  3*  It  < 

<ujjo  4X4>-*i-o— */>43<n— ujcoo  4  »  *  o 

0—4  —■*0—0—  *x  •  OUU  UJ*t/l 

— tO*--4-|  *3  tt  II  i/XUJt/iaJI-l—  3Z&  •X'-OZ  <2 

1/140  4  H  Z  •>—  40X  Zh*2«av»«4«  *  U— 

40—*—  X-J  lUUJtOO  •4——  •  U>UJ  • 

OUJ  ZOMK4U.ZX4  UUZUjUJUJaCUJIkbO-Jt-  *Q  UlOtO 

xxt— •*— 4®*"— o*-at'4JJa3M4«  <oat— xac< 


1  /1 4  X 

< 

to—  •* 

coat** 

400 

3 

>  at  uj 

XX 

</>U  < 

cn 

too.  a. 

33  » 

<  «UJIU 

< 

> 

ZZ  'I 

Ot— 3 

o 

•*— 

Q 

•  -JU.i/1 

H-l— 

X  UJQ 

I-X0JO3 

at 

UJt-0  » 

—  -IUJ 

1-  X 

o 

o— * 

—  -I_i 

II  —  X</»- 

u. 

at_jto* 

3  — -J 

X3UJ 

4</14* 

It— 

/i  Z.y~<-r> 

4 

1-40 

l/l  5£ 

3—tUOC 

ZH--IO. 


O  «Q  -Jt-XO> 

•►-M-tL 


4UJ 

oz 


SlU 

HZ 


LL.UJUJ  >-l— UJLUUJ 


0—1 


at 

a. 


O  l-UJX  -XU 

OtOCtl-l-l-O  JJI3J 
UJt—  Z  *ZUja.o;  ♦  i— 

-J  O  JUX- **-4 

— Oujujuj  at  uju.— 03  #4 

u-o-i_j.ja.at— tc-izi-  h- 
Z  IIH-Z 

H»— l/ltO— 

IU0C4— Of 

•  jCLUJOi 
•  a.  • 

o  a  i— 

uu  a 

at  — «  -j  -  —  - 


4 

a 


o 

x 

/i 


o— 

/I 

24  Z— 
UJ— o*  — o 

X  •*-! 

— t/M—  ♦ 

1-40* 

O*- 


i-o  UJ  UJ 

UJZ  30 
-IUJ 

»ac 

. J 

(M  U.U. 

OJ  CC 


XX 

30 

z< 


at 

a 

u. 


O.Z  at  at 

*  —  UJUJ 

*  stat  33  aa 

*  1/10. 


22 


169 


w 

• 

* 

* 

%n 

• 

* 

3 

• 

* 

X 

• 

• 

* 

* 

• 

* 

3 

UI 

* 

O 

3 

■ 

* 

X 

«« 

X 

o 

H* 

o 

II 

UI 

• 

a. 

<-» 

M 

(A 

UI 

L9 

N 

UI* 

x 

►- 

M 

UI 

►- 

<* 

o 

o 

UI 

UI 

X 

UI  < 

o* 

1— 

3 

< 

M 

3  £ 

•» 

X 

M 

u 

►- 

x  ar 

UI* 

H' 

lo 

o 

<* 

*mrn 

* 

a. 

i— 

o* 

# 

M 

a 

M.  LU 

•* 

• 

m 

LO0CZ 

• 

M  _J  OO 

Z  UI 

V9 

• 

<iuo 

f 

z 

PM 

Z  H  0>> 

in 

(A 

• 

H- 

m 

UI* 

O 

1-  H  II 

1—  tl  4* 

X 

• 

_i 

Ul 

uiXO 

3* 

M 

ui  O 

OILI  xz 

•J 

u>* 

< 

X 

LU30 

X  • 

H— 

X  •  r’I'O 

•  X  03 

_ i 

a) 

3* 

t— 

«Z^ 

►-* 

O 

ui  m  atoc 

H-m  < 

5 

(A  • 

X* 

o 

a 

• 

z 

<  M  OI« 

H*  b- 

(/» 

H- 

►- 

►— 

3UIP*> 

3 

o  i-  3>a* 

•  UI 

«c 

Uw» 

M 

«b 

«X<S> 

M* 

u. 

•  <  _j 

M  4p4 

o 

# 

o 

00 

3*-aC 

UI 

**  X  GDI  1  '■ 

UI  1/1  | 

(/I 

QZ* 

—3 

Ui 

</>  M 

<  H 

< 

a.  xujuj 

<  >  -* 

UJUJ  — 

IU»—# 

ulX 

'  X 

<  .LL 

O 

X  <— >— i 

o*<s 

at*-* 

>— ♦— 

►» 

3 

UUI 

3 

u» 

I-  H*  X 

-#*QZ 

^•JUJ 

ao 

.4 

■* 

•UI 

-X 

< 

UlM  LU 

aui*  _iui  h 

J(OM 

< 

ui3t 

ZmX 

3oo 

LUZ  -JO  II  000-***<toc-* 

3  K 

LU 

• 

3 

UI 

OH" 

X< 

UI 

M  Q 

o*  =  * 

^Z4 

(AZ 

3 

MUI* 

UIO 

UI 

1—  UI  M  H  || 

•  •*  LU  Z 

<UQ 

<m4 

X 

< 

l-Z* 

QUIZ 

< 

M* 

-jz*o  -* 

1— UI  _i  H 

OUMA 

1— 

O 

4m« 

uiZ** 

O 

<m«  >MZ(MUIUI3  II  O  1- 

-»<r 

4} 

a 

33* 

l/)M 

o 

3o*MO  cc.ce. 

*—l—  >— UJ 

_ 

QHJ 

UJ**# 

►— 

z 

LO 

<306 

■* 

*; 

UI  •  M  OLO_J 

UiUi 

• 

UI 

a>  3 

UI* 

3 

<t— H»  C2  UI»>UII-U>  < 

acoo 

-4>- 

ui* 

UI 

u 

LUO 

<• 

UI 

o  *Z>--i 

<UI<m40 

UU 

CQZ 

<* 

UIXO 

o* 

< 

«H>«MCO  Mt— O  U  — 

QUi« 

*fc»* 

u* 

h— 

Z 

mmQ 

• 

o 

ztz  OLUUJ 

i—  «  UJ 

zz 

cC 

UI 

v» 

p— 

m  OUImIOZI-UI 

</» 

<MZ(A 

a. 

P-* 

ai 

« 

<  H* 

II  ui 

uui 

uiac  ►-wr»QC 

3 

HU3 

ui 

a. 

MM 

X  3 

<3 

►—CL  Ui  «X 

jKUIA  • 

X 

LU  OlX 

(A  • 

XI 

301 

UI 

VKXX 

ouiX 

OX 

M  4  4*0000  «^ui 

H* 

352*- 

MM 

►— 

-KO 

ZX 

i-m 

h-  - 

mp4 

— LUH* 
ZM 

Oh* 

Z  •  OX  <UIUI03 
OP*  »-  il  ll  UZ3-I* 

</> 

fld  LUUI 

z 

ZUJ 

Z>N 

l-OZW 

1— UI 

UI 

H  o  II 

LU 

QQU1 

30 

ui 

3 

ZOUI 

LU3UJ 

UUUI 

o  *0  p*^-XJ<«I- 

Z 

LLU 1  z 

t— 

uih*  wituu. 

a 

UIM  M 

z<z 

az 

oo*^HM«au  i—o  ui 

ai/Nu« 

O 

z_iuu 

zo 

op- 

2MI- 

-O*- 

M 

-<►— *  II  Q:tf  <<<M(/IX 

<X-i 

HiM 

~«<z 

— « 

—030 

■JUI4 

XJ 

II  *  a.  .Z»Ui  3LL 

XCDm 

M  * 

-O— 

~x 

»  —1* 

JJU< 

< 

< 

<  ►— <M 

_J  It  * 

UI  U1_l 

1  < 

3  *  00  JU3W-<U<V 

0£fM 

«*  m  h  H-p-f-acz 

ac< 

< 

LU— *< 

mac 

<  H 

Lunaa^oi 

3 

>  H  X  LUUJUJOm  u. 

01-01- 

Z  U(N 

i  a 

> 

Z  LL  <ZH-3I- 

□H 

at  ujJli*  — 

z 

o<zz 

ac« 

M*— 

t-Q 

Ml-  >U 

mZ<Z 

H*Z 

HIM  1— 

Q. 

h-0<M 

LU 

_jZ_ia-Zi- 

LU 

OZlUOi 

JM|-M 

l/)M 

y~  Mice 

X 

u>*»  ac 

►-« 

MA 

>-« 

—xxu. 

acooe 

M« 

Zee  30 

Q. 

•*  a. 

ZO 

(Mae»~)Xae»« 

ZO 

(MOO— ^afMQ.H-0. 

xa. 

—a  ll 

X 

Mlb 

agM^x 

MU. 

a. 

LL 

a.  h» 

a. 

a. 

M. 

* 

xi 

* 

•  M 

IA 

UI 

oo 

•  m 

170 


171 


LIST  OF  REFERENCES 


1.  Hardison,  D.  C..and  others.  Special  Study  Group, 
Departaent  of  the  Aray,  Review  of  Aray  Analysis.  Voluae 
I  -  Main  Report,  section  3-5,  April,  79/9” 

2.  Lucas,  G. ,  "Generalized  Differential  Equations  of  Coa- 

bat",  A Jy 11a bag  af  fat  gcoaoajc,  £er sogneland 

Force  Effectiveness  sfufl^gs,  p. zETji HGeceaBer  T9 71 . 


Ibid.,  p.  263. 

G  - 

p.  14,  The 


4.  Clark,  G.  M. ,  T£e  qoabat  Analysis  Model.  Ph.D.  Thesis, 
Ohio  State  university,  Coluabus,  Ohio,  1969. 


5. 


6.  Clark,  G.  M.,  Tfee  Coabat  Analysis  Model.  Ph.D.  Thesis, 
jj.^149.  The  Ohio  Sfafe  university,  Coluabus,  Ohio, 


7.  Taylor,  J.  G. ,  Forge- on  Fojfie  Attyitiop  Modelling.  Mil¬ 
itary  Applicatigns  Sect ion  of  the- Operations  Research 
Society  of  Aaerica,  p.  29,  Arlington,  Virginia,  1980. 

8.  Ibid.,  p.  33. 

9.  Ibid.,  p.  77. 


172 


13. 

14. 


ipple  coapater  Inc. ,  Hi£  £25  flftaflalzBiSfe  SRSEaliaS  §!§' 
isi#  ^980. 


U  &nd  iia®  of  Light 
£1,  Navalrpostgraduate 
paper. 


173 


BIBLIOGRAPHY 


Apple  Computer  Inc. ,  £hg  Applesoft  fraflguage  Manual.  1978. 

Kiviat,  P.  J. t  Villanueva,  R.  and  Markowitz,  H.M.,  SIMSCBIPI 
II^5c|rograa»ing  Lan^u^ge.  2nd  Edition,  Consolidated'  Analy- 

Taylor,  J.  G.,  Lanche star-type  fifi^gls  o£  Warfare,  draft  edi¬ 
tion,  Voluae  2,  Haval  Postgraduate  School, Monterey,  CA, 
1980. 


INITIAL  DISTB IBOTION  LIST 


Defense  Technical  Inforaation  Center 
Caaeron  station 
Alexandria,  Virginia  22314 

Library,  Code  0142 
Naval  postgraduate  School 
Monterey,  California  93940 

Departaent  Chair aan.  Code  55 
Department  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  California  93940 

Chief 

TRADOC  Research  Eleaent  Monterey 
Naval  Postgraduate  School 
Monterey,  California  93940 


Associate  professor  Jaaes  K.  Hartaan 
Code  55Hh 

Departaent  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  California  93940 


Associate  Professor  S.  H.  Parry 
Code  55Py 

Departaent  of  Operations  fiesearch 
Naval  Postgraduate  School 
Monterey,  California  93940 

Headquarters 

(JS  Army  Training  and  Doctrine  Coaaand 
ATTN:  Director,  Studies  and  Analysis 

Directorate,  Mr.  S.  Goldberg 
Port  Monroe,  Virginia  23651 


Associate  Professor  A.  L. 
Code  53Zh 

Departaent  of  Matheaatics 
Naval  Postgraduate  School 
Monterey,  California  93940 


Schoenstadt 


No.  Copies 
2 


Professor  Jaaes  G.  Taylor 
Code  55Tw 

Departaent  of  Operations  Besearch 
Naval  Postgraduate  School 
Monterey,  California  93940 


175 


1 


10. 


Office  of  the  Coaaanding  General 
OS  Aray  Training  and  Doctrine  Coaaand 
ATTN:  General  Gl^ng  Otis 
Port  Honroe,  Virginia  23651 


11.  Headquarters  1 

OS  Aray  Training  and  Doctrine  coaaand 
ATTN:  ATCG-T  (BG  HorelllT 

Port  Honroe,  Virginia  23651 


12. 


Office  of  the  Coaaanding  General 
OS  Readiness  Coaaand 
ATTN;  General  Donn  A  Starry 
NacDill  APB,  Florida  33621 


1 


13.  Hr.  falter  Hollis  1 

Deputy  Onder  Secretary  of  the  Aray 
(Operations  Research) 

Departaent  of  the  Aray.  The  Pentagon 
Washington,  D.  C.  20310 


14.  LTG  Howard  Stone  1 

Coaaanding  General 
OS  Aray  Coabined  Aras  Center 
Fort  Leavenworth,  Kansas  66027 


15.  Director  1 

Coabined  Aras  Coabat  Development  Activity 
ATTN:  ATZL-CAC- A  (Hr.  Lee  Pleger) 

Port  Leavenworth,  Kansas  66027 


16. 


17. 


Director 

Coabat  Analysis  Office 
ATTN:  Hr.  Kent  Pickett 

OS  Aray  Coabined  Aras  Center 
Port  Leavenworth,  Kansas  66027 


Coaaand  and  General  Staff  College 
ATTN:  Education  Advisor 
Room  123,  Bell  Hall 
Port  Leavenworth,  Kansas  66027 


1 


1 


18.  Dr.  Wilbur  Payne,  Director  1 

OS  Aray  TRADOC  Systems  Analysis  Activity 
White  Sands  Hissile  Range,  New  Hexico  88002 


19. 


Headquarters,  Departaent  of  the  Army 
Office  of  the  Deputy  Chief  of  Staff  for 
Operations  and  Plans 
ATTN:  DAHO-2D 

Washington,  D.C.  20310 


1 


20.  Coaaander  1 

05  Aray  Concepts  and  Analysis  Agency 
ATTN:  HOCA-WG 

8120  Woodaont  Avenue 
Bethesda,  Haryland  20014 


176 


1 


21.  Director 

US  Aray  Material  Systems  Analysis  Activity 
ATT H:  DRXS7-<PM  JBr.  Bill  Nieaeyeri 

Aberdeen  Proving  Grounds,  Maryland  21005 


22.  Director  1 

US  Aray  Might  Vision  and  Electro-Optical 
Lab 

ATTN:  DEL-NV-VI  (Hr.  Prank  Shields) 

Port  Belvoir,  Virginia  22060 


23.  Director  1 

USATRASANA 

ATTN:  Mr.  Ray  Heath 

White  Sands  Missile  Range,  New  Mexico  88002 


24. 


Director 

Coabat  Developments,  Studies  Division 
ATTN:  MAJ  W.  Scott  Wallace 

OS  Aray  Armor  Agency 
Port  Knox,  Kentucky  40121 


1 


25.  Coamandant  1 

US  Army  Field  Artillery  School 
ATTN:  iTSA-CD-DSWS 
Port  Sill,  Oklahoma  73503 


26. 


Director 

Coabat  Developments 
US  Army  Aviation  Agency 
Fort  Rucker,  Alabama  36 


36  2 


1 


27.  Director  1 

Coabat  Developments 
US  Aray  Infantry  School 
Fort  Benning,  Georgia  31905 


28.  Director  1 

Combat  Developments 
ATTN:  ATZA-CDE  (CPT  James  Mudd) 

US  Aray  Engineer  School 
Port  Belvoir,  Virginia  22060 


29.  Director  1 

Coabat  Developments 
ATTN:  ATSA-CDF-S 

US  Army  Air  Defense  Agency 
Fort  Bliss,  Texas  79905 


30.  Commander  1 

US  Aray  Logistics  Center 
ATTN:  ATCL-OS  (Mr.  Caaaeron/CPT  McGrann) 

Port  Lee,  Virginia  23801 


177 


1 


31.  Director 

Combat  Developments 
ATTN:  LTC  Harnagel 

OS  Army  Signal  School 
Port  Gordon,  Georgia  30905 


32.  Deputy  Onder  Secretary  of  the  Army  1 

for  Operations  Research 
Room  2E261,  Pentagon 
Bashington,  DC  20310 


33.  MAJ  Jeffrey  L.  Ellis,  OSA  1 

Department  of  Operations  Research,  Code  55EI 
Naval  Postgraduate  School 
Monterey,  California  93940 


34.  CPT  Charles  J.  McKenzie  III  1 

OS  Readiness  Coaaand 
HacDill  AF  Base,  Plorida  33621 


35.  CPT  Thoaas  J.  Pawlowski  III  1 

Department  of  Mathematical  Sciences  (DFHS) 

OS  Air  Force  Acadeay,  Colorado  30840 


36. 


Nr.  Ernest  Boehner 


CACDA/CASAA 
Ft.  Leavenworth,. 


Kansas  66027 


37. 


Nr.  David  Goodwin 

1301  Stillhouse  Creek  Bd 

Chesterfield,  Missouri  63017 


1 


38.  CPT  Jaaes  V.  Nudd  1 

35  Heussy  Ave. 

Buffalo,  New  Tork  14220 


39.  CPT  Bruce  Hock  1 

288  Alexander  Drive 
Manchester,  New  Hampshire  03103 


40.  CPT  Nilliam  T.  Farmer 

TRADOC  Research  Element  Monterey 
Naval  Postgraduate  School 
Monterey,  California  93940 


178 


DATE 

FILMED 


DTIC 


